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CONVEKSION  FACTORS,  U.S.  CUSTOMARY  TO  Ml-TRIC  (SI) 
UNITS  OF  MIiASUFU:MI;NT 


U.S.  customary  units  of  measurement  used  in  this  report  caji  be  converted 
to  metric  (SI)  units  as  follows: 


Multiply 

by 

To  obtain 

indies 

25.4 

mi  llimeters 

2.54 

centimeters 

square  inches 

6.452 

square  centimeters 

cubic  inches 

16.  39 

cubic  centimeters 

feet 

30.48 

centimeters 

0.3048 

meters 

square  feet 

0.0929 

square  meters 

cubic  feet 

0.0283 

cubic  meters 

yards 

0.9144 

meters 

square  yards 

0.836 

square  meters 

cubic  yards 

0.7646 

cubic  meters 

mi  les 

1.6093 

kilometers 

square  miles 

259.0 

hectares 

knots 

1.8532 

kilometers  per  hour 

acres 

0.4047 

hectares 

foot-pounds 

1.3558 

newton  meters 

mi  1 libars 

1.0197  X 10' 

^ kilogriuns  per  square  centimeter 

ounces 

28.35 

grams 

pounds 

453.6 

gr.ams 

0.4536 

kilograms 

ton,  long 

1.0160 

metric  tons 

ton,  short 

0.9072 

metric  tons 

degrees  (angle) 

0.1745 

radians 

Fahrenheit  degrees 

5/9 

Celsius  degrees  or  Kelvins^ 

^To  obtain  Celsius 

(C)  temperature 

readings  from  Falirenheit  (F)  readings. 

use  formula:  C = 

(5/9)  (F  -32). 

To  obtain  Kelvin 

(K)  readings,  use 

formula:  K = (5/9)  (F  -32)  + 273.15. 
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SYMBOLS  AND  DEFINITIONS 


A cross-sectional  area  of  a channel 

Aj,  effective  surface  area  of  a block. 

Ac  cross-sectional  area  of  a channel 

As  surface  area  of  an  estuary  at  MSL 

ap  amplitude  of  input  tide  to  an  estuary 

B 8/3it  m(AsCo)2  Rq,  a parameter  which  determines  the  phase  lag  of 

tidal  response  in  an  estuary 

BN  right-hand  side  of  equation  (48) 

BP  right-hand  side  of  equation  (46) 

b (8A/9s)j.]  const.!  a characteristic  of  a channel 

C(j  dimensionless  discharge  coefficient  characterizing  a constricted 

opening  between  bay  and  sea 

Cg  admittance  coefficient  (with  dimensions  of  velocity);  nominally 

represents  the  wave  speed  in  the  sea 

Cq  dimensionless  overflow  coefficient  (generally  less  than  0.5 

for  a broad-crested  barrier) 

Cs  dimensionless  discharge  coefficient  for  a submerged  barrier 

(generally  less  than  ^ ) 

D total  depth  of  water  at  position  x,y  at  time  t 

D a mean  depth  for  the  effective  fetch  across  a block;  also  mean 

depth  for  a channel  (Dm  + Dp)/2 

D},  depth  of  water  over  the  crest  of  a barrier 

Dp  effective  depth  of  a channel  A^/w 

maximum  depth  to  be  expected  anywhere  in  the  system  during  a 
storm  surge 

h'l  contribution  to  the  forcing  term  in  equation  (17)  due  to  lateral 

transfer  of  mass  and  momentum 

f dimensionless  bed  resistance  coefficient  for  blocks 

fj.  channel  bed  friction  coefficient 


II 


SYMBOLS  AND  DEF INITIONS--Continued 
G damping  factor  for  channels,  see  equation  (44) 

Gj  damping  factor  for  x-transport  on  blocks,  see  equation  (35) 

G2  damping  factor  for  y-transport  on  blocks,  see  equation  (36) 

g acceleration  due  to  gravity 

H water  level  elevation  relative  to  local  MSL  datum 

HB  water  elevation  on  the  water-connected  block  of  a channel 

lie  common  water  elevation  for  a channel  junction 

HM  mean  water  level  anomaly  of  connected  channel  and  blocks 

HX  water  level  at  the  lower  end  of  an  x-channel 

HY  water  level  at  the  left  end  of  a y-channel 

H at  point  B in  a channel 

water  level  on  the  high  side  of  a barrier 

Hg  input  tide  level  at  time  t outside  a bay  entrance 

ll(i,j)  water  level  anomaly  H for  block  identified  by  x and  y 

indexes  i,j 

H*  tentative  predicted  H for  a ponding  block  in  the  absence  of 

any  contribution  by  longitudinal  discharge  to  or  from  the 
channel  which  terminates  adjacent  to  that  block 

H'  value  of  H at  new  time  level 

new  H value  at  point  1 in  channel 

Hj  6 water  levels  on  the  two  sides  of  a barrier  (both  of  which 

exceed  Zb),  equation  (10) 

i x-index  for  grid  blocks 

j y-index  for  grid  blocks 


dimensionless  wind-stress  coefficient,  equation  (6) 
effective  fetch  length 


SYMBOLS  AND  DEFINITIONS- -Continued 


Lf  net  time  rate  of  gain  of  water  volume  per  uiiit  distance  along 

the  channel  by  lateral  transfer  and  rainfall 

1^  net  time  rate  of  gain  of  momentum  (divided  by  water  density) 

per  unit  distance  along  channel 

m fL/gDcA^  or  l/gCCdAd)^ 

N denotes  negative  characteristic 

n time  index 

P wind  "push"  term  XAt  or  YAt;  also  denotes  positive 

characteristic 

Q volume  transport  through  cross-sectional  area  of  a channel 

Q mean  Q value  for  channel,  equation  (45) 

QCXP|^  flow  at  the  upper  end  of  an  x-channel  for  channel  block  K 

QCYN|^  flow  at  the  left  end  of  a y-channel  for  channel  block  K 

QCYP|^  flow  at  the  right  end  of  a y-channel  for  channel  block  K 

QCXN|^  flow  at  the  lower  end  of  an  x-channel  for  channel  block  K 

Q at  point  A of  positive  characteristic 

Qg  Q at  point  B of  negative  characteristic 

Q^j  discharge  from  channel  to  ponding  block 

qf  the  flow  (per  unit  length  of  channel)  from  the  channel  to 

the  adjacent  block 

q^  lateral  volume  flux  per  unit  length  into  the  channel 

qn  outward  component  of  volume  flux  at  a boundary 

qo  lateral  volume  flux  per  unit  length  out  of  the  channel 

qj,  flow  (per  unit  length  of  channel)  from  the  channel  block  to 

the  channel  (across  the  interior  side  of  the  channel) 


Q' 

new  Q 

value 

new  Q 

at  point  N 
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SYMBOLS  AND  DEFINITIONS--Continued 


Q{>  new  Q at  point  P 

QJ.  specified  river  discharge 

R rainfall  rate 

R(i,j)  rainfall  rate  for  block  i,j 

r relative  amplitude  response 

s distance  along  the  axis  of  a channel 

T tidal  period 

Tg  longitudinal  component  of  wind  stress  (divided  by  water  density) 


appropriate  wind-stress  component  (X  or  Y)  corresponding  to 
time  level  t for  the  associated  channel  block 


U vertically  integrated  x-component  of  volume  transport  per 

unit  width 

UCF(K)  lateral  transport,  per  unit  width  per  unit  time,  nominally  from 
an  x-channel  of  block  K to  an  adjacent  block;  also  denoted 
UCFk 

UCT(K)  lateral  transport,  per  unit  width  per  unit  time,  nominally  to 
an  x-channel  from  the  interior  of  block  I;  also  denoted  UCT^ 

UN  U value  on  left  side  of  block 

U(i,j)  value  of  U at  the  left  side  of  block  i,j 

UCi+l,;)  value  of  U at  the  right  side  of  block  i,j 

u typical  fluid  speed  in  the  bay 

U'  value  of  U at  new  time  level 

V vertically  integrated  y component  of  volume  transport  per 

unit  width 


VCF(K)  lateral  transport  per  unit  width  per  unit  time,  nominally  from 
an  y-channel  of  block  K to  an  adjacent  block;  also  denoted 
VCFk 


SYMBOLS  AND  DEF INITIONS--Continued 


VCT(K) 


VN 


V(i,j) 
V(i,j  + 1) 
V 
W 


y 


z 

Z(iJ) 

Zb 

Zc 

a 

r 

AH 

Aq 


lateral  transport  per  unit  width  per  unit  time,  nominally  to 
an  y-channel  from  the  interior  of  block  K;  also  denoted  VCTj( 

value  of  V at  the  lower  side  of  a block 

value  of  V at  the  lower  side  of  block  i,j 

value  of  V at  the  upper  side  of  block  i,j 

value  of  V at  new  time  level 

windspeed  at  10-meter  elevation  over  the  water 

a critical  speed  taken  as  14  knots  (7  meters  per  second) 

surface  width  of  a channel  (conveyance  width) 

x-component  of  the  wind  stress  divided  by  the  density  of  the 
water 

value  of  X for  right  side  of  block  i,j 

horizontal  Cartesian  coordinate  nominally  alongshore,  positive 
to  the  right  when  facing  shore 

y-component  of  the  wind  stress  divided  by  the  density  of  the 
water 

value  of  Y for  top  side  of  block  i,j 

horizontal  Cartesian  coordinate  nominally  normal  to  shore, 
positive  landward 

elevation  of  the  seabed  relative  to  MSL  datum 
value  of  Z for  block  i,j 
barrier  crest  elevation 
channel  bed  elevation 

(gD)**  At/As  (Courant  number);  also  Lc/DcAc,  equation  (77) 
L(CbDb)^/^At 

a head  differential  dependent  upon  barrier  type 

net  lateral  flow  to  the  channel  per  unit  length  of  channel 
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SYMBOLS  AND  DEFINITIONS- -Continued 


As  grid  size  for  blocks  (distance  between  successive  H values 

in  both  the  x and  y directions);  also  written  AS  or  DELS 

At  time  step  (time  interval  between  successive  H values  at  given 

location);  also  written  DELT 

6 the  angle  between  the  wind  velocity  vector  and  the  x-axis 

X w (gn‘)VG 

If  3.14159  ... 

o wflQl/A^ 


Ip  latitude 

n absolute  angular  speed  of  the  earth 

u radian  frequency  2fr/T 


DEVELOPMENT  OF  SURGE  II  PROGRAM  WITH  APPLICATION  TO  THE 
SABINE-CALCASIEU  AREA  FOR  HURRICANE  CARU  AND  DESIGN  HURRICANES 


by 

Robept  0.  Reid,  Andrew  C.  Vaatano,  and  Thomas  J.  Reid 
I.  INTRODUCTION 


Numerical  techniques  for  the  solution  of  equations  representing  storm 
surges  in  coastal  areas  were  significantly  augmented  in  1966  by  the  de- 
velopment of  a two-dimensional  model  (referred  to  in  this  study  as  the 
SURGE  I program)  for  the  U.S.  Army  Engineer  District,  Galveston  (Reid  and 
Bodine,  1968).  At  about  the  same  time  a number  of  bay  models  emeiged. 
Notable  among  these  are  the  models  of  Leenderste  (1967)  and  Masch,  et  al. 
(1969),  which  have  been  applied  to  problems  of  both  surge  and  circulation 
in  bays.  These  models  include  the  Coriolis  force  which  is  neglected  in 
the  Reid-Bodine  model.  However,  the  Reid-Bodine  model  produced  the  first 
successful  inclusion  of  flooding,  recession,  barriers,  and  flow  over 
barriers  in  the  study  of  inundation  of  low-lying  coasts.  The  actual 
model  is  a nonlinear  system  of  equations  and  boundary  conditions  solved 
by  numerical  integration  of  time -dependent,  forced  motion.  Its  use  pro- 
duces the  water  response  to  stormwinds  over  the  region  for  a given  storm 
tide  at  the  seaward  boundary.  The  initial  application  was  a hindcast  of 
the  Hurricane  Carla  surge  generated  in  Galveston  Bay  during  9 to  12  Sep- 
tember 1961. 

During  Hurricane  Carla,  the  wetted  perimeter  of  Galveston  Bay  essen- 
tially doubled,  as  accurately  reproduced  in  the  hindcast  computations. 
Serial  observations  of  water  levels  for  the  storm  period  available  from 
stations  throughout  the  bay  were  compared  to  levels  computed  with  the 
numerical  algorithm.  These  records  produced  a standard  deviation  of  less 
than  4 inches,  overall.  The  maximum  deviation  of  the  water  level  predic- 
tion was  1.5  feet  and  occurred  at  the  grid  square  corresponding  to  the 
location  of  the  Pelican  Island  Bridge  which  spans  the  channel  between 
Galveston  and  the  Pelican  Islands.  Although  this  disparity  was  rela- 
tively large,  its  effect  on  the  computations  was  effectively  reduced  by 
the  smoothing  operation  of  the  numerical  integration.  However,  this 
difference  points  out  a basic  problem  confronting  any  model — the  minimum 
definition  of  topographic  features. 


The  basic  problem  of  indicating  subgrid  scale  effects  in  numerical 
modeling  is  normally  solved  by  parameterization  of  the  omitted  physical 
mechanism.  Often,  an  analytic  relationship  is  introduced  that  requires 
the  specification  of  empirically  derived  constants;  e.g. , the  wind-stress 
equation  for  the  transfer  of  momentum  from  wind  to  water.  Another  sinple 
and  pertinent  instance  is  the  a priori  rotation  of  wind  vectors  over 
certain  grid  squares  in  the  Hurricane  Carla  computations  for  Galveston 
Bay.  The  model  Galveston  entrance  channel  was  not  in  the  proper  orienta- 
tion on  the  Cartesian  numerical  grid  system  and,  as  a result,  did  not 
admit  a realistic  amount  of  water  to  the  bay.  A programed  shift  in  the 
wind  vectors  indicated  this  subgrid  scale  feature. 
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SURGE  I has  been  applied  to  the  study  of  Texas  coastline  surge  sus- 
ceptibility. The  topographic  features  of  this  region  are  characterized 
by  barrier  islands  and  shallow,  river-fed  bay  systems  surrounded  by  near 
sea  level  land  and  marshes.  The  specific  applications  of  the  program 
have  therefore  centered  interest  on  the  immediate  environs  of  a bay.  The 
requirement  for  surge  studies  of  appreciable  distances  inland  from  the  bay 
system  has  only  recently  been  placed  on  the  numerical  model.  The  propa- 
gation of  the  surge  to  higher  ground  throu^  necessary  subgrid  scale  topo- 
gr^hic  features  has  required  an  extension  of  the  basic  algorithm. 

The  new  algorithm  developed  for  the  study  of  the  Sabine-Calcasieu 
region  is  referred  to  as  the  SURGE  II  program.  TTiis  program  incorporates 
all  the  features  of  SURGE  I with  the  further  option  of  representing  vari- 
able depth  and  width  channels  along  the  sides  of  eadi  grid  square.  The 
flow  computations  for  the  channels  interact  with  the  normal  grid  square 
confutations  and  permit  a complete  suite  of  flooding  conditions  for  over- 
topping of  levees.  In  this  manner  SURGE  II  provides  a time-dependent, 
subgrid  scale  transport  of  water  through  the  model. 

II.  THEORETICAL  DEVELOPMENT  FOR  SURGE  II 


1.  Summary  of  Two-Dimensional  Theory. 

The  development  of  SURGE  II  was  based  on  the  SURGE  I concept  by  Reid 
and  Bodine  (1968).  A part  of  this  study  is  presented  here  to  provide  a 
complete  description  of  SURGE  II. 

The  advection  of  momentum  (or  field  acceleration)  is  considered  neg- 
ligible except  at  singular  regions  of  the  bay  (submerged  barriers  and 
narrow  channels)  where  the  effect  is  included  implicitly  throu^  the  use 
of  appropriate  nonlinear  discharge  relations.  TTie  effect  of  the  earth's 
rotation  is  also  neglected;  this  approximation  appears  justifiable  for 
systems  of  small  spatial  scale  and  shallow  depth  t^ere  frictional  forces 
are  more  dominant. 


Within  the  normal  domain  of  the  bay  and  immediate  adjoining  sea,  the 
vertically  integrated  equations  of  motion  and  of  continuity  appropriate 
to  the  problem  are  taken  as  follows: 


f +gDf  -X-fqUD-^ 


(1) 


f +gD|^-Y-fqVD“^ 


(2) 


iH  + 3V  _ 
at  3x  3y  ' 


(3) 
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where 


X and  y = horizontal  Cartesian  coordinates; 
t » time; 

U and  V = vertically  integrated  x and  y components,  respectively, 
of  transport  per  unit  width; 

g » gravity; 

H = water  level  elevation  relative  to  the  local  mean  sea 
level  (MSL)  datum; 

D = depth  of  water  at  position  x,  y at  time  t; 
q = magnitude  of  the  transport  per  unit  width; 
f = dimensionless  bed- resistance  coefficient; 

R = rainfall  rate; 

X and  Y = X and  y components  of  the  wind  stress  divided  by  the 
density  of  the  water  (the  density  assumed  constant) . 

Normal  values  of  f are  in  the  range  10"^  to  lO"^  for  typical  seabed 
conditions . 

The  value  of  q is  obtained  from  U and  V by 

q a (u2  + (4) 

which  is  a positive  quantity. 

The  kinematic  forms  of  the  wind-stress  coiq)onents  in  the  absence  of 
rainfall  are  taken  as 

X ■ K cos  0 

Y - K W2  sin  0 , (5) 

where  N is  the  windspeed  at  a 10-meter  elevation  over  the  water,  and 
0 is  the  angle  between  the  wind  velocity  vector  and  the  x-axis.  The 
dimensionless  coefficient,  K,  used  in  the  calculations  is  presumed  to 
be  a function  of  windspeed  as  implied  by  the  van  Dom  (1953)  relation 
for  wind  stress.  Specifically,  it  is  assumed  that 

K - Ki  for  W < W 

* — c 

W 2 

K - Ki  + K2  (l  - ^)  for  W > , (6) 
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where  the  constants  Kj  and  K2  are  taken  as  1.2  x lO"®  and  1.8  x 10"®, 
respectively,  and  Wj.  is  a critical  speed  which  is  taken  as  14  knots 
(7  meters  per  second).  For  large  windspeeds,  K approaches  the  limiting 
value  of  3;6  x lO"^  which  corresponds  to  a resistance  coefficient  of  about 
3.0  X 10”^  if  the  ratio  of  air  density  to  water  density  is  taken  as 
1.2  X 10-3. 

In  the  presence  of  rainfall  an  added  flux  of  momentum  proportional  to 
RW  occurs  (van  Dom,  1953).  The  effect  can  be  included  by  augmenting  K 
by  R/W.  For  heavy  rainfall,  the  resulting  K is  increased  about  10  per- 
cent. 

The  variables  H and  D are  related  by  the  single  expression, 

D = H - Z , (7) 

where  Z is  the  elevation  of  the  seabed  relative  to  the  MSL  datum. 
Presumably,  Z is  a function  of  x and  y only;  i.e.,  the  time-depen- 
dent  scour  of  the  seabed  is  ignored. 

The  above  equations  ignore  the  direct  effect  of  variable  atmospheric 
pressure  which  is  relatively  minor  in  a small,  shallow  bay.  The  effect 
over  the  sea  is  included  implicitly  through  the  specification  of  an 
appropriate  surge  height  versus  time  in  the  adjoining  sea  where  the  com- 
bined effects  of  winds  and  differential  atmospheric  pressure  give  rise 
to  a coastal  storm  surge.  This  is  presumed  to  be  determined  independently 
of  the  detailed  calculations  for  the  bay  and  enters  as  a boundary  condi- 
tion. 


a.  Boundary  Conditions.  Four  different  types  of  boundary  conditions 
are  used  in  this  system  of  computations.  Two  of  these  conditions  apply 
to  the  water- land  boundary,  one  condition  applies  to  the  artificial  bound- 
ary representing  the  seaward  end  of  the  bay  system,  and  one  applies  at 
partial  barriers  internal  to  the  system.  (Additional  internal  conditions 
are  needed  in  the  presence  of  imbedded  channels  as  discussed  later  in 
Section  III, 2.)  All  four  conditions  relate  the  normal  component  of  flow 
at  the  boundary  to  the  state. of  the  water  level  at  the  boundary. 


In  general,  the  boundary  between  bay  water  and  Icind  depends  on  the 
water  elevation  and  the  land  topography.  The  shoreline  for  different 
uniform  elevations  of  the  surface  of  the  bay  is  readily  established  from 
a knowledge  of  the  topography.  For  a bay  with  low-lying  terrain,  the 
rate  of  increase  of  surface  area  of  water  per  unit  increase  of  water 
level  can  be  considerable.  In  the  actual  rising  stage  of  storm  tide  the 
amouit  of  inundation  is  controlled  by  the  rate  at  which  the  water  can 
flow  into  the  potential  ponding  areas.  In  the  present  scheme,  which  uses 
a representation  of  the  bay  in  terms  of  a discrete  grid,  the  elevation  of 
the  seabed  or  land  is  regarded  as  uniform  over  each  grid  square,  thus 
forming  a two-dimensional,  stairstep-type  approximation  of  the  actual 
topography.  The  boundary  condition  on  the  normal  consonant  of  flow,  q^, 
at  the  juncture  of  a flooded  square  and  a dry  square  is  taken  as 
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(8) 


\ “ O' 

if  the  elevation,  H,  of  the  water  is  less  than  that  of  the  adjacent 
dryland.  However,  if  the  water  level  is  greater  than  that  of  the  dryland, 
then  the  rate  of  flooding,  per  unit  length  of  land  barrier,  is 

given  by 

<1„  - ± 

where  Dj,  is  water  depth  over  the  crest  of  the  barrier,  and  Cq  is  an 
appropriate  dimensionless  overflow  coefficient,  generally  less  than  0.5 
for  a broad-crested  barrier.  The  choice  of  sign  depends  on  whether  the 
flooding  is  from  bay  to  land  or  from  flooded  land  back  to  the  bay  during 
the  recession  stage. 

Equation  (9)  is  considered  valid  for  any  barrier  within  or  at  the 
boundary  of  the  system  for  which  the  water  level  on  one  side  of  the 
barrier  is  greater  than  the  barrier  crest  elevation,  Z^,,  and  for  which 
the  water  level  on  the  other  side  is  less  them  Zb.  Moreover,  Db  is 
simply  Hb  - Zb,  where  Hb  is  the  water  level  on  the  high  side. 

In  the  case  where  the  water  level  on  both  sides  of  an  internal  barrier 
exceeds  the  barrier-crest  elevation,  the  discharge  is  taken  as  that  for  a 
submerged  wier, 

‘In  ” ± - Hal)**  . (10) 

where  Db  is  the  water  depth  over  the  crest  of  the  barrier,  Hj  and  Ha 
are  the  water  levels  on  the  two  sides  of  the  barrier  (both  of  which  exceed 
Zjj) , and  Cg  is  an  appropriate  dimensionless  discharge  coefficient  for 
the  submerged  barrier  (generally  less  than  /2) . In  this  case,  is 

taken  as  (H^  ♦ H2)/2  - Zb.  Again,  the  sign  is  taken  such  that  the  flow 
is  directed  toward  the  low-head  side  of  the  barrier.  Both  equations  (9) 
and  (10)  presume  that  the  velocity  of  approach  to  the  barrier  is  much 
less  than  the  velocity  over  the  barrier. 

In  the  numerical  computational  scheme,  emphasis  is  placed  on  the  eval- 
uation of  flow  and  water  levels  within  a bay  which  is  connected  to  a sea 
of  essentially  unlimited  extent.  An  appropriate  boundary  condition  is 
required  either  at  the  mouth  of  the  bay  system  or  along  some  line  within 
the  sea  which  delineates  the  outer  limit  of  the  computational  grid.  The 
correct  approach  would  be  to  treat  the  development  of  the  surge  in  the 
sea  and  bay  as  a single  problem.  However,  the  difference  in  spatial  res- 
olution required  for  the  two  different  regions  of  the  system,  as  well 
as  coiqjuter  storage  limitations,  makes  this  inpractical.  The  assumption 
is  made  that  the  effect  of  the  conditions  in  the  bay  has  only  a minor 
influence  on  the  development  of  the  surge  in  the  sea  and  over  tlie  Conti- 
nental Shelf.  The  evaluation  of  the  latter  can  be  determined  independ- 
ently of  the  bay  problem  or  obtained  from  observation  and  used  as  an 
outer  boundary  condition  for  the  bay. 
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The  simplest  condition  at  the  seaward  boundary  is  of  the  form 

H = Hg  . (11) 

where  Hg  is  the  prescribed  water  level  which  would  exist  in  the  absence 
of  the  bay  at  time  t at  the  outer  boundary  of  the  bay  system.  SURGE  II 
presently  uses  this  condition  at  the  seaward  boundary  and  at  lateral 
boundaries  on  the  limited  shelf  part  of  the  system.  An  alternative  con- 
dition for  the  lateral  boundaries  on  the  shelf  is  to  prescribe  that 
3U/3x  = 0 at  these  boundaries  where  x is  taken  alongshore  (Jelesnianski , 
1966,  1967).  An  alternative  condition  for  the  seaward  boundary  is  one 
which  allows  for  radiation  of  energy  to  the  sea.  The  latter  condition 
is  of  the  form 


H = Hg  + q^/Cg  . (12) 

where  qjj  is  taken  positive  outwards  from  the  bay  to  the  sea,  and  Cg 
is  an  appropriate  admittance  coefficient  (with  dimensions  of  velocity). 
Nominally,  Cg  represents  the  wave  speed  in  the  sea.  The  generalized 
condition  (eq.  12)  is  nearly  equivalent  to  the  simplest  condition  (eq. 

11)  if  Cg  greatly  exceeds  the  wave  speed  for  the  bay. 

b.  Initial  Conditions.  Since  the  system  includes  allowance  for  fric- 
tional dissipation  as  well  as  radiation  of  energy,  the  solution  for  given 
fields  of  X and  Y and  given  boundary  function,  Hg,  should  be  reason- 
ably insensitive  to  the  nature  of  the  initial  conditions  after  a suitable 
lapse  of  time  from  the  initial  state.  Thus,  the  initial  conditions  can 
be  somewhat  arbitrary.  As  in  the  laboratory  model  experiments,  it  is 
reasonable  to  start  from  a state  of  equilibrium  in  which  U and  V are 
zero  and  H is  uniform  throughout  the  system,  in  order  to  minimize  the 

introduction  of  transient  oscillations  related  to  the  starting  conditions. 

Moreover,  a reasonable  period  (depending  on  the  characteristic  decay  time) 
can  be  allowed  for  the  system  to  reach  that  state  where  its  response  re- 
flects only  the  effect  of  the  forcing  functions. 

2.  Theory  of  Embedded  Channels. 

Let  s denote  distance  along  the  axis  of  a channel  whose  cross- 
sectional  area  is  A and  surface  width  is  w at  position  s and  time 
t.  Let  Q be  the  volume  transport  througji  A in  the  positive  sense  of 
s,  and  let  H be  the  water  elevation  above  MSL  datum  at  the  same  section. 

In  general,  A and  w are  known  functions  of  H for  a given  cross  sec- 

tion, as  determined  by  the  geometry  of  the  cross  section  (Fig.  1).  In 
particular,  3A/3H  = w for  given  s.  The  width  w is  to  be  the  "convey- 
ance" width,  as  used  by  Dronkers  (1964). 

The  channel  is  considered  an  "open  system"  in  the  sense  that  water 
and  momentum  may  enter  or  leave  the  channel  laterally;  i.e.,  exchange  of 
fluid  with  adjacent  bay  area  or  flooded  land  can  exist.  If  the  longi- 
tudinal velocity  in  the  channel  is  considered  uniform  for  evaluating  the 
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Figure  1.  Schematic  channel  cross  section  showing 
pertinent  parameters. 


longitudinal  transport  of  momentum,  then  the  equations  of  motion  and 
continuity  for  a given  channel  reach  are  (Stoker,  1957,  Ch.  11;  Dronkers, 
1964,  Ch.  9) 


3t  3s 


(q2/a) 


A 3H 
8^^ 


wT  - oQ  + L 
s m 


(13) 
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(14) 


where 

Tg  = longitudinal  component  of  wind  stress  (divided  by  water 
density) ; 

0 = wfjQl/A^  where  f is  a dimensionless  channel-friction 
coefficient ; 

Lf  “ net  time  rate  for  gain  of  water  volume  per  unit  distance 
along  the  channel  by  lateral  transfer  and  rainfall; 

* associated  net  time  rate  of  gain  of  momentum  (divided  by 
water  density)  per  unit  distance  along  channel. 

The  units  of  Lf  are  square  feet  per  second;  Lj^  has  the  units  cubic 
feet  per  second  squared. 

It  is  convenient  in  the  analysis  of  the  channel  dynamics  to  trans- 
form the  above  equations  into  a characteristic  form.  There  are  several 
different  possible  characteristic  forms.  The  approach  used  by  Stoker 
(1957)  is  to  work  with  u and  H (where  u = Q/A)  as  the  dependent 
variables.  Dronkers  (1964)  works  with  either  Q and  H directly  or 
with  Q and  total  head  (H  (Q/A)^/2g).  Each  method  has  certain 
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advantages  and  disadvantages.  In  the  present  analysis,  the  variables 
Q and  11  are  used  to  be  as  consistent  as  possible  with  the  computa- 
tions in  the  two-dimensional  regions  of  the  system. 


r 


In  transforming  equations  (13)  and  (14)  to  characteristic  form,  it 
is  noted  that 

aA  _ 3H 

at  ” at 


3A 

as 


b , 


(15) 


where 


b 


const  . 


(16) 


(For  a channel  of  uniform  cross  section  the  latter  quantity  would  be 
zero.)  It  can  be  shown,  following  Dronkers'  (1964)  analysis  and  con- 
sidering equation  (15),  that  a characteristic  form  of  equations  (13) 
and  (14)  is 


- oQ  + L + b 
m 


(17) 


along  the  path  s(t)  where 


^ = S + f&h.  . 

dt  A — w 


(18) 


The  path  line  where  the  plus  or  minus  sign  is  taken  in  equation  (18)  is 
referred  to  as  the  positive  P characteristic  or  the  negative  N 
characteristic  path,  respectively.  These  are  illustrated  in  Figure  2 
where  x corresponds  to  s,  the  two  paths  having  point  C in  common. 
Equation  (17)  with  the  upper  sign  applies  along  P and  equation  (17) 


) 


with  the  lower  sign  applies  on  path  N.  Thus,  information  with  regard 
to  Q and  H at  points  and  Xg  at  time  t and  along  the  two 

paths  can,  in  principle,  be  lised  to  predict  the  values  of  Q and  H 
at  point  C from  two  equations. 

For  a laterally  closed  channel  {Lf,  Ljjj  = 0)  of  a uniform  cross  section 
(b  = 0)  without  friction  (cr  = 0)  , in  the  i)sence  of  wind  stress  (T^  = 0) , 
then  the  quantity  in  braces  on  the  right-hand  side  of  equation  (17)  van- 
ishes. In  this  case  only  the  information  at  points  A and  B of  Figure 
2 is  needed  to  predict  values  of  H and  Q at  C.  To  show  that  equa- 
tion (17)  is  consistent  with  Stoker's  (1957)  analysis  for  this  special 
case,  let  u = Q/A  and  D = A/w.  For  a uniform  cross  section  at  given  H, 
dll/dt  = dD/dt,  so  equation  (17)  reduces  to 

i(DU)  + (_  u + ^ = 0 


along 


= u i 


(20) 


Equation  (19)  simplifies  further  to 


wD  ^ (u  + 2/^)  = 0 . (21) 

dt  — 

Thus,  for  this  special  case  (u  + 2»'gD)  is  conserved  along  P where 
dx/dt  = u + /gD,  while  (u  - 2/gD)  is  conserved  along  N where  dx/dt  = 
u - /gD.  Thus,  u and  D (hence,  Q and  H)  can  readily  be  evaluated 
at  C. 

In  the  more  general  case  the  time  integral  of  the  right-hand  side  of 
equation  (17)  must  be  estimated  in  a rational  way.  This  is  considered 
later  in  Section  III, 2.  Also,  in  the  general  case  it  is  usually  not 
possible  to  put  the  left-hand  side  of  equation  (17)  in  the  sinple  form 
shown  in  equation  (21) . 

a.  Lateral  Transfer  Terms.  In  the  absence  of  direct  rainfall,  Lf 
must  equal  the  net  gain  of  volume  per  unit  length  per  unit  time  due  to 
lateral  flow  into  the  channel  on  either  or  both  sides.  J.et  q^  and  qQ, 
respectively,  represent  the  volume  fluxes  per  unit  length  into  and  out  of 
the  channel.  Then,  Lf  = di  ■ 9o  absence  of  rainfall,  or 

Lf  = qi  - % ^ wR  (22) 

with  rainfall.  The  corresponding  lateral  transfer  of  momentum  (divided 
by  water  density)  is 


L 

m 


(23) 
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the  transfer  from  rainfall  being  included  in  the  wind-stress  term  as 
discussed  in  Section  11,1.  In  equation  (23)  the  quantity  Uq  is 
simply  Q/A  for  the  channel  while  is  the  channel -directed  com- 

ponent of  velocity  of  fiLiid  from  the  adjoining  block  water  area.  In 
equation  (17)  the  terms  and  Lf  contribute  to  the  right-hand  side 

the  quantity. 


_ S L + /fiA 


L - ^ L 
m A f 


(24) 


Using  equations  (22)  and  (23)  yields 


Trap 

■ “o^  " "^“o  - (wj  ■ *^0 


(25) 


The  lateral  flows  into  or  out  of  the  channel  can  be  evaluated  by 
relations  such  as.  equations  (8) , (9),  and  (10).  This  is  also  discussed  in 
Section  111,2. 


b.  Simplifications.  The  SURGE  II  program  uses  certain  simplifica- 
tions of  the  above  equations.  For  normal  conditions,  the  propagational 
speed  (gA/w)^  significantly  exceeds  the  speeds  ui  or  Uq;  i.e.,  Q/A. 
Accordingly,  F is  approximated  by 


F = + 
Li 


w 


* 


(26) 


Elsewhere  in  equations  (17)  and  (18),  Q/A  is  neglected  compared  with 
(gA/w)^.  Moreover,  each  channel  reach  within  a grid  block  is  considered 
of  uniform  width  and  bottom  elevation  Zc 5 however,  w and  Zc  vary 
from  one  reach  to  another.  Thus,  b = 0 for  each  reach  and 

A/w  = D = H - Zc  . (27) 

Under  these  conditions  equations  (17)  and  (18)  take  the  form, 

^ = {wT^  - f (q(q/^2„)  + " *10 

along 

+ /gD  (29) 


ds 

dt 


where  Tg  » X or  Y as  s ■ x or  y,  depending  on  channel  orientation.  Equa- 
tion (28)  can  also  be  expressed  in  the  form. 


^ (0  i-l-  wo/gD  ) - F 


(30) 
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for  a given  channel  reach  where  F is  the  right-hand  side  of  equation 
(28).  The  neglect  of  Q/A  relative  to  in  the  above  approximate 

channel  equations  is  tantamount  to  neglect  of  longitudinal  advection  of 
momentum  in  the  original  equation  (13),  an  approximation  already  made  in 
the  two-dimensional  equations  in  Section  II,  1. 

III.  SURGE  II  PROGRAM 

Numerical  algorithms  for  two-dimensional  blocks  and  subgrid  scale 
channels  are  given  in  this  section,  and  the  coupling  between  these  is 
discussed.  A conplete  listing  of  the  SURGE  II  program  is  in  Appendix 
A.  A description  of  the  program,  as  adapted  for  the  GE-400  computer, 
and  the  required  input  and  output  options  are  discussed  in  /^pendix  B. 
Appendix  C is  a user’s  guide  to  the  SURGE  II  program.  The  block  algo- 
rithm is  essentially  as  discussed  by  Reid  and  Bodine  (1968)  except  for 
a change  in  the  barrier  computation  and  incorporation  of  coupling  with 
the  subgrid  scale  channels. 

1.  Block  Algorithm. 

In  the  numerical  analog  of  the  prognostic  equations  (1),  (2),  and 
(3),  values  of  H are  evaluated  on  a uniform  Cartesian  mesh  at  spacing. 
As,  for  uniform  time  steps.  At.  The  values  of  H are  representative 
of  the  water  level  for  the  grid  square  i,  j which  is  centered  at 
X = (i  - 1/2)  As,  y = (j  - 1/2)  As,  at  time  nAt,  in  which  i,  j,  and 
n are  integers.  Values  of  Z are  specified  as  permanent  storage  for 
the  same  locations  as  H so  that  D can  be  evaluated  as  needed  at  these 
locations.  Values  of  U are  evaluated  at  even  half  steps  of  x,  odd 
half  steps  of  y,  and  odd  half  steps  of  t (Fig.  3).  This  staggered 
system  gives  the  least  storage  consistent  with  a given  spatial  resolution. 
It  corresponds  to  the  simplest  scheme  discussed  by  Platzman  (1958)  and 
requires  only  half  the  storage  compared  with  the  coupled  scheme  used  by 
Miyazaki  (1963). 

The  variables  X and  Y are  supplied  at  spatial  locations  consis- 
tent with  U and  V,  respectively,  but  at  even  half  steps  of  t.  Values 
of  Hg  are  supplied  for  positions  and  times  on  the  outer  boundary  of  the 
bay  consistent  with  the  locations  and  times  for  the  H values  on  that 
line.  Values  of  R are  supplied  at  locations  consistent  with  H but 
at  a one-half  time  step  out  of  phase  with  H.  Arrays  of  X,  Y,  and  R, 
for  a single  value  of  j and  n,  and  the  array  of  Hg  values  for  given 
n are  read  from  tape  as  required.  The  fields  of  X hnd  Y are  gene- 
rated from  a coarse  spatial  and  temporal  array  evaluated  from  the  basic 
meteorological  data  and  then  evaluated  for  the  detailed  mesh  by  linear 
interpolation. 
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Figure  3.  Exanple  of  grid  blocks  showing  staggered 
arrangement  of  variables  U,  V,  and  H. 

Information  pertinent  to  the  position,  elevations,  and  discharge 
coefficients  for  barriers  (those  not  resolved  by  the  limitations  of  the 
grid  system)  is  stored  as  permanent  storage  along  with  the  field  of  Z. 

The  numerical  analogs  of  equations  (1),  (2),  and  (3)  use  values  of 
U,  V,  H,  Z,  X,  Y,  and  R at  locations  shown  in  Figure  4 for  a typical 
calculation.  In  the  present  application  a common  valiie  of  R for  given 
time  is  used  for  the  whole  spatial  array.  The  following  notation  is  used 
in  the  recursion  equations:  H(i,j)  represents  H centered  in  block 
i,  j at  t * nAt;  U(i,j)  represents  U for  the  left  side  of  block  i,j  at 
t = (n  - 1/2)  At;  v(i,j)  represents  V for  the  lower  side  of  block  i,j 
at  t » (n  - 1/2)  At. 

Primed  symbols  are  used  to  denote  values  of  these  variables  at  time 
step  At  later.  Thus,  the  difference  U'  - U is  centered  in  time  at 
the  level  of  H,  and  the  difference  H'  - H is  centered  in  time  at  the 
level  of  U'  or  V.  The  notation  for  Z or  D is  consistent  with 
that  of  H. 
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Figure  4.  Basic  block  triad  showing  variables  used 
in  confutation  of  U,  V,  and  H for  block  1. 

The  frictional  terms  in  equations  (1)  and  (2)  are  represented  by 
fAU'D"  and  fQV'D"^,  respectively,  where  the  estimation  of  Q and  D 
is  centered  spatially  at  the  position  for  U'  or  V.  Since  U,  V,  and 
D are  not  available  at  common  locations,  this  requires  a suitable  spatial 
average  in  order  to  obtain  centered  values  of  Q and  D.  The  resulting 
recursion  equations  for  U,  V,  and  H,  using  centered  differences  for 
the  spatial  derivatives,  are  as  follows: 

+ ffl 

+ D(i,j)l  lH(l,j)  - H(i+l,j)]  + X(i+l,j)At}  (31) 

“ Gl(i;jy  + ff~  [DCi.J+D 

+ D(l,j)]  lH(i,J)  - H(1,J+1)1  + Y(i,j+l)At}  (32) 


^ [U'(i,j)  + V'(i,j) 

- U'(i+l,j)  - V'(l,j+1)]  + R(l,j)At,  (33) 

where 

D(i.j)  “ H(l,j)  - Z(i,j).  (34) 

and  Gi  and  G2  are  the  factors  which  incorporate  the  effect  of  the 
friction.  These  are  given  by: 

Gi(i,j)  = 1 + f At  {[4U(i+l,j)]2  + [V(i,j)  4 V(i+l,j) 

+ V(l.j+1)  + VCi+l.j+Dl^}**  [D(i,j)  + D(i+l,J)l"^  (35) 

and 


Gzd.j)  - 1 + fAt  {[4V(l,j+l)2  + [U(l,j)  + U(i+l,j) 

+ U(i,j+1)  + U(i+l,j+l)]2}*5  [D(i,j)  + D(i,j+1)]"^  . (36) 

The  latter  factors  are  always  somewhat  greater  than  unity  unless  the  flow 
or  friction  factor  vanishes. 

The  prediction  relation  for  H given  by  equation  (33)  does  not  con- 
sider any  possible  contribution  of  flow  to  or  from  the  block  due  to  the 
presence  of  a subgrid  scale  channel.  This  will  be  considered  in  a subse- 
quent section. 

It  should  also  be  emphasized  that  the  effect  of  Coriolis  force  is  not 
considered.  The  relative  importance  of  the  Coriolis  force  compared  with 
bottom  friction  can  be  estimated  in  terms  of  the  ratio,  r,  of  these  two 
forces  which  is  of  the  order, ' 

r * XD/fu  , (37) 


where 

X » Coriolis  parameter  (2£1  sin  ip,  (I  being  the  absolute  angular 
speed  of  the  earth  and  p the  latitude); 

D = mean  depth; 

'f  ■ bottom- friction  coefficient; 
u * typical  fluid  speed  in  the  bay. 
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For  30*  latitude  X = 7.3  x 10”^;  typical  D and  f for  gulf  coast  bays 
are  10  feet  and  2 x io“^,  respectively.  For  u = 3 feet  per  second,  which 
is  reasonable  for  storm  conditions,  r is  only  1/10.  However,  for 
normal  circulational  regimes  u may  be  only  a fraction  of  1 foot  per 
second  and  r is  of  order  unity.  Hence,  while  it  may  be  justifiable  to 
neglect  the  Coriolis  term  for  short-duration  storm  surge  studies  for 
shallow  bays  of  limited  horizontal  dimensions  it  cannot  be  neglected  in 
long-term  circulational  studies. 

Althougji  it  does  not  appear  difficult  to  add  the  effect  of  Coriolis 
force,  it  can  be  shown  (Platzman,  1958)  that  a different  scheme  for  the 
U,  V,  and  H arrays  is  necessary  for  numerically  stable  con^)utations  using 
an  explicit  time-marching  procedure  eis  used  here.  The  coupled  scheme 
required  for  stable  explicit  computations  at  least  doiijles  the  computing 
time.  The  present  scheme  could  be  used  with  an  implicit  time-marching 
procedure  to  maintain  stability  and  similar  accuracy,  but  this  too  can 
be  achieved  only  af  the  cost  of  an  increase  in  computing  time  by  a factor 
of  at  least  two.  In  the  presence  of  friction,  the  destabilizing  effect 
of  the  Coriolis  terms  in  an  explicit  scheme  such  as  that  used  by  Masch 
(1969)  is  suppressed;  however,  this  is  acconplished  only  at  the  sacrifice 
in  rendition  of  the  frictional  terms.  Thus,  the  omission  of  the  Coriolis 
force  from  a program  intended  primarily  for  gulf  coast  estuaries  is  moti- 
vated primarily  for  reasons  of  economy  of  operation,  in  respect  to  surge 
calculations . 

a.  Stability.  Numerical  stability  requires  that  At  be  taken  at 
less  than  the  value  AS/(2gDj,33j)^,  where  is  the  maximum  depth  to 

be  expected  anywhere  in  the  system  during  the  storm  surge  (Platzman, 

1958) . 


b.  Barrier  Algorithm.  Equations  (9)  and  (10)  are  assumed  to  apply 
for  values  of  q^^,  D^,,  and  AH  at  the  same  time  and  in  the  immediate 
vicinity  of  the  barrier.  In  the  grid  scheme  used,  however,  the  flow  and 
the  water  level  are  staggered  in  time;  moreover,  the  water  levels  like 
Hi  and  H2  represent  in  effect  the  spatial  average  for  blocks  1 and  2, 
respectively,  at  a given  time  rather  than  local  values  in  the  vicinity 
of  a given  barrier,  which  in  the  schematization  are  presumed  to  occur  on 
lines  separating  two  blocks.  As  a consequence  the  above  relations  can- 
not be  applied  directly.  Instead,  the  evaluation  of  U or  V across 
a barrier  (if  the  water  level  allows  such  flow)  is  carried  out  by  a 
modified  version  of  the  predictive  equations  (1)  and  (2),  or  their 
numerical  counterparts,  equations  (31)  and  (32),  where  f is  replaced 
by  an  effective  value  related  to  the  barrier  discharge  coefficient  so  as 
to  be  consistent  with  equations  (9)  or  (10).  The  effect  is  to  maintain 
proper  time  phasing  and  to  consider  possible  tilt  of  water  level  across 
the  block;  i.e.,  difference  of  H at  barrier  relative  to  the  mean  value 
for  the  block. 

^ecifically,  the  frictional  terms  in  equation  (1)  or  (2)  are  taken 
as  (®/LC^) |q^|q^/D^  where  is  the  barrier  discharge  coefficient 
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(Cq  or  Cg,  depending  on  type  of  barrier),  is  the  transport  per  unit 

width  normal  to  the  barrier  (either  U'  or  V* , depending  on_barrier 
orientation)  , Djj  is  the  water  depth  over  the  barrier,  and  D is  a 
mean  depth  for  the  effective  fetch  L across  the  blocks.  The  gravita- 
tional slope  term  involves  the  same  scale  length,  L,  and  mean  depth, 

D.  The  resulting  relation  for  prediction  of  at  a barrier,  given 

q^  at  the  previous  time  step,  is: 

kAkA  + rqA  - F . (38) 

where 


and 


r = 


DAt 


(39) 


F = g(CbDb)2AH  + r • (q^  + P)  , (40) 

P being  the  wind  "push"  term  (XAt  or  YAt) , and  AH  a head  differential 
dependent  on  barrier  type.  For  steady  state  (q^  = q^)  and  no  wind  (P  = 0) , 
the  above  reduces  to 


q;  = ± CfeDb  . (41) 

which  is  consistent  with  equation  (9)  or  (10)  with  and  AH  taken 

as  Cq  and  Dj,  or  Cg  and  (Hj  - H2) , respectively,  depending  on  the 
barrier.  The  more  general  relation  (eq.  38)  provides  an  added  effect  of 
the  wind  and  of  the  inertia  of  the  water  on  the  blocks.  For  a submerged 
barrier,  L is  taken  equal  to  AS;  i.e.,  from  the  center  of  block  1 to 
the  center  of  block  2.  For  an  overflow  barrier,  L is  taken  as  half 
this  distance  since  the  inertia  and  wind  setup  are  effective  only  on  the 
higher  of  the  two  blocks. 

Thus,  Cj,,  L,  H,  and  Dj,  are  taken  as  follows: 

Submerged  barrier  (Hj  > Zj,  and  H2  > Z^) 


Cb  = Cg 

L =»  AS  AH  = Hi  - H2 


Db  =[fHl  ^ Zb 

Overflow  barrier  (Hj  > Zb  or  H 
Cb  “ Cq 

L » AS/2  AH  » - 

Db  » |ah1 


2 > Zb) 

Hi  - Zb 

or 

Zb  - H2 


(a) 

(b) 


» 


(42) 


T 
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where  Zjj  is  the  elevation  of  the  barrier  crest,  relation  (a)  being  for 
Hi  3b  and  (b)  for  H2  > Zb-  If  Zb  exceeds  both  Hi  and  H2,  then 
= 0.  The  meaningful  solution  of  the  quadratic  equation  (38)  is 

qA  = ± nlF|  + (r/2)2]»i  - r/2},  (43) 

where  the  sign  is  taken  as  that  of  F,  as  verified  from  equation  (38). 

The  above  relations  for  barriers  differ  from  that  used  in  Reid  and 
Bodine  (1968)  and  in  the  original  SURGE  I program.  The  present  barrier 
relations  have  a more  realistic  response  when  applied  to  the  numerical 
simulation  of  a natural  oscillation  of  a bay  having  a submerged  barrier 
across  it. 

c.  Barrier  Specification.  Since  only  certain  blocks  contain  barriers, 
they  must  be  identified  by  I,  J location;  specifically,  the  program 
identifies  the  Kth  barrier  block  by  location  I = IB(K)  and  J = JB(K) , 

K = 1,2  ...  KM.  A given  barrier  block  potentially  has  a barrier  on  the 
right  and  upper  side  of  the  block  in  an  x,  y plot.  These  are  designated 
X and  y,  respectively;  i.e.,  an  x barrier  is  one  normal  to  the  x-axis 
(the  flow  over  it  being  in  the  x sense) . For  both  potential  barriers 
on  a barrier  block,  values  of  Zjj,  Cq,  and  Cg  must  be  prescribed.  A 
real  barrier  is  one  where  is  larger  than  the  Z value  for  either  of 

the  adjoining  blocks.  A null  barrier  is  one  where  Z5  equals  the  larger  of 
the  Z values  for  the  adjoining  blocks  (thus,  in  effect,  the  higher  block 
is  a potential  barrier) . The  program  requires  that  information  pertinent 
to  both  null  barriers  (Z^,  Cq,  and  Cg)  and  real  barriers  be  provided. 

d.  Volume  Check.  During  the  recession  stage  of  flooding  when  water 
is  draining  off  flooded  blocks  (via  the  barrier  overflow  relation) , it  is 
possible  for  the  volume  leaving  in  one  time  step  as  computed  from  qA^t 
to  exceed  the  available  volume.  Therefore,  a test  is  included  in  the 
program  such  that  if  this  occurs,  the  flow  is  adjusted  to  only  drain 

the  block  dry  (D  = 0) , and  the  flow  to  adjacent  blocks  adjusted  to  be 
consistent . 

e.  Depth  Check.  When  the  water  depth  is  very  shallow  the  effect  of 
the  wind  is  such  that  a given  block  could  become  partially  dry  unless  the 
fluid  is  flowing  fast  enough  for  the  bottom  stress  to  balance  the  wind 
stress.  To  avoid  anomalous  computations  for  very  small  D (e.g.,  in 
areas  where  rainfall  is  occurring  over  regions  above  the  surge  level) , 
the  wind  stress  is  arbitrarily  set  zero  when  D is  less  than  0.1  foot. 

2.  Channel  Algorithm. 

a.  Channel  Specification.  As  in  the  case  of  barriers,  those  blocks 
on  which  channels  occur  are  identified  by  the  I and  J values;  for 
channel  block  K these  are  denoted  by  ICG(K)  and  JCG(K),  respectively, 
where  K = 1,2  ...  KCM.  Also  each  "channel  block"  may  contain  two  channels, 
one  on  the  right  denoted  the  x channel  and  one  on  the  upper  side  denoted  • 
the  y channel.  Each  of  these  channel  reaches  is  characterized  by  a 
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channel  width  (v) , a channel-bed  elevation  (Zj.) , and  a channel- 
friction  coefficient  (f{;)«  Figure  5 shows  a schematic  of  a channel 
block  indicating  nomenclature  for  dimensions  as  used  in  the  SURGE  11 
program.  Figure  6 shows  the  dependent  variables  pertinent  to  the 
channels  as  used  in  the  program  and  stored  for  the  channel  block  K. 
These  include  the  channel  flows,  Q,  at  each  end  of  the  channel,  one 
end  designated  N,  the  other  P (corresponding  to  the  negative  and 
positive  characteristic  ends  of  the  channel,  respectively).  Also  in- 
cluded is  the  height,  H,  of  the  water  level  at  the  point  in  common  to 
the  two  channels  for  block  K (UC(K)).  The  lateral  transport  (per  unit 
width  per  unit  time)  nominally  to  the  channel  from  block  K and  from 
the  channel  is  also  indicated:  UCT(K)  and  UCF(K),  respectively,  for  the 
channel  normal  to  the  x-axis,  and  VCT(K)  and  VCF(K),  respectively,  for 
the  channel  normal  to  the  y-axis.  In  the  formulas  in  this  study,  these 
are  referred  to  as  and  qf,  respectively.  Note  that  UCF(K)  and 

VCF(K)  correspond  to  U and  V,  respectively,  on  the  right  and  upper 
sides  of  the  general  block  flow.  Also,  the  quantity  HP(K)  corresponds 
to  the  block  (pool)  height  for  the  channel  block.  Values  of  H at  the 
"negative"  ends  of  the  channels  for  channel  block  K are  stored  as  HC 
values  in  adjacent  channel  blocks  to  minimize  duplication  of  storage. 

b.  Computation  of  Channel  Variables.  The  time  phasing  of  block 
variables  versus  channel  variables  is  indicated  in  Table  1.  The  H 
values  occur  at  common  times  thus  facilitating  evaluation  of  head  dif- 
ferentials used  in  determining  lateral  flow  between  channel  and  adjacent 
blocks . 


Table  1.  Time  phasing  of  computations 
for  blocks  and  channels. 


Time 

Block 

Channel 

t ♦ At 

H 

H,Q 

t ♦ At/2 

Q 

t 

H 

H,Q 

t - At/2 

Q 

t - At 

H 

H,Q 

For  a given  channel  reach,  application  of  equations  (28)  and  (29) 
can  be  made  for  two  characteristic  paths,  as  shown  schematically  in 
Figure  7.  As  in  the  case  of  the  block  computations,  the  friction  term 
in  equation  (28)  is  taken  proportional  to  the  product  of  a new  Q and 
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BLOCK  K 
I=ICGk 
J=JCGk 


Figure  6.  Channel  block  K at  coordinates  I = ICG(K)  and 
J = JCG(K) , showing  associated  flows  and  water 
level  variables. 


N P 


5 > 


Figure  7,  Characteristic  paths  on  the  time-distance 
diagram  for  an  individual  channel  reach. 

the  absolute  value  of  the  old  Q.  Specifically,  for  the  positive  char- 
acteristic path  from  A to  P'  in  Figure  7,  equation  (28)  is  approxi- 
mated by 

■ V " V “ •^Aq]  At  , (44) 

where  D = (Dj^  + Dp)/2,  Tg  is  the  appropriate  wind-stress  component 
(X  or  Y)  corresponding  to  time  level  t for  the  associated  channel 
block,  Aq  is  the  net  lateral  flow  per  unit  width,  and  Q is  taken  as 

Q •=  KQS  + Qp)/2]‘®  . (45) 

The  subscripts  on  Q,  H,  and  D designate  the  points  at  which  these 
apply  (see  Fig.  7)  and  primes  denote  new  time  level. 

After  regrouping  terms,  equation  (44)  can  be  written  as 

+ (w  /Oir  = [(Q^  + w H^)  + (WT^  + At]/G  , (46) 

where 

G = 1 + f^At|Q|/(D)2w  . (47) 

Similarly,  for  the  negative  characteristic  from  B to  N' , 

- (w*^ /G)ll^  = ((Q^  - w.^ Hg)  + (WTg  - Y^Atl/G  , (48) 

where  D and  G are  as  defined  for  the  positive  characteristic. 


The  values  of  Q and  11  at  points  A and  B are  determined  by 
interpolation  from  values  at  N and  P at  time  t,  using  equation 
(29)  for  the  path.  The  distance  from  A to  P or  _B  to  N,  using 
the  mean  wave  speed  for  the  channel  at  time  I is  /gllAt.  The  interval 
N to  P is  equal  to  As.  Let 


a 


t^D  At/As  ; 


(49) 
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this  should  always  be  less  than  or  at  most  unity  for  stability  of  com- 
putation. The  linearly  interpolated  values  at  A and  B are  then 


Qa  “ “ Qn  Qp 

Qb  “ (1  -a)  Qn  + « Qp  . 

and  similarly  for  and  Hg  in  terms  of  Hf^  and  Hp. 


(50) 


The  evaluation  of  Aq  is  the  most  sensitive  part  of  the  computations 
and  is  discussed  in  a subsequent  section.  Presuming  Aq  is  known,  the 
problem  of  evaluating  the  new  Q and  H individually  at  the  channel- 
end  points  is  considered.  Note  that  equations  (46)  and  (48)  yield  pre- 
dictions for  linear  combinations  of  Q and  H at  two  different  points. 
Thus,  information  from  adjoining  channels,  or  other  information  in  the 
case  of  channel  end  points,  is  needed  to  solve  for  the  new  channel  Q 
and  H.  For  a simple  continuous  channel  without  branches  and  consisting 
of  a series  of  reaches  of  length  As  but  not  necessarily  of  equal  width 
or  depth,  then  Q and  11  are  readily  solved  at  a common  junction,  using 
the  information  from  the  positive  characteristic  from  one  channel  and  the 
negative  characteristic  from  the  adjoining  channel.  However,  branches  do 
occur  and  it  is  therefore  desirable  to  use  a sufficiently  general  pro- 
cedure which  will  accommodate  either  brcinching  channels  or  continuous 
channels. 


In  the  scheme  chosen  for  representing  channels  in  SURGE  II  it  is 
possible  to  have  four  channels  merging  at  a common  junction.  Figure  8 
shows  this  jujiction  with  four  different  volume  transports,  but  with  a 
common  H.  The  designation  of  the  different  Q shown  in  this  figure  is 
that  used  in  the  coded  program  (see  App.  B) ; QC  for  channel  transport, 

X or  Y denoting  the  channel  (not  the  direction  of  flow),  and  N or  P 
denoting  whether  the  flow  is  at  the  negative  or  positive  end  of  a given 
channel  reach.  Each  is  identified  by  a channel  block  index  K. 


For  any  given  channel  reach  equations  (46)  and  (48)  predict,  for  a 
given  point,  values  of  the  quantities 


BP  = Q'  + XH' 


BN  = Q'  - XH'  , 


where 


X = w v^/G  . 


(51) 


(52) 
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Figure  8.  General  channel  junction,  showing 
flows  and  channel  identification. 
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V 


For  simplicity  of  notation  let  1, 
with  1 being  the  lower  channel, 
channel,  and  4 the  right  channel 


2,  3,  and  4 denote  the  merging  channels 
2 the  left  channel,  3 the  upper 
(Fig.  8).  Then,  with  this  notation 


Ql'  + 

XI 

• H 

' = BPl 

Q2'  + 

X2 

• H 

' = BP2 

(53) 

Q3'  - 

X3 

• H 

' = BN3 

Q4'  - 

X4 

• H 

' = BN4  . 

Now,  continuity  requires  that 

Ql'  + Q2' 

- 

Q3' 

- Q4'  =0 

(54) 

at  a common  junction.  Thus, 

(XI  + X2  + X3  + X4) 

H' 

' = BPl  + BP2  - BN3  - BN4 

(55) 

from  which  11'  can  be  calculated  at  the  junction.  With  11'  known,  the 
values  of  Ql',  Q2',  Q3',  and  Q4'  are  readily  evaluated  from  equation  (40). 


For  those  cases  where  one  or  two  of  the  above  merging  channels  do  not 
exist  (i.e.,  null  channels),  then  their  width  and  X value  are  zero. 
Moreover,  the  program  yields  zero  for  the  BP  or  BN  values  of  any  null 
channel.  Thus,  equations  (55)  and  (53)  can  apply  for  a general  junction 
consisting  of  from  two  to  four  real  merging  channels. 


c.  Net  Lateral  Flow.  The  net  time  rate  of  water  accumulation  in  the 
channel  per  unit  length  due  to  lateral  exchange  with  blocks  and  by  rain- 
fall is 


= qt  ■ ^f  '''R  « (56) 

where  q^  corresponds  (if  positive)  to  the  flow  (per  unit  length  of 
channel)  from  the  channel  block  to  the  channel  (across  the  "interior" 
side  of  the  channel.  Fig.  6)  and  q^  (if  positive)  is  the  flow  (per 
unit  length  of  channel)  from  the  channel  to  the  adjacent  block.  These 
flows  can  be  positive,  negative,  or  zero.  To  allow  for  channels  which 
have  widths  w much  smaller  than  the  block  grid  size  As,  and  since  the 
above  q values  are  comjjarable  to  tliose  \ hich  exist  across  the  sides  of 
blocks,  tlie  change  in  channel  water  level  can  be  vc'y  sensitive  to  the 
difference  q^  - qf.  Hence,  special  care  must  be  taken  in  the  nwdel  to 
avoid  possible  instal)i lities  caused  by  improper  calculation  of  these 
transverse  flows.  However,  there  is  no  particular  difficulty  witli  the 
rainfall  term  in  equation  (Sb)  whicli  is  generally  at  least  one  order  of 
nuagnitude  smaller  than  tliat  of  tlie  "net"  lateral  flow.  In  a sense,  tlie 
potential  difficulty  with  the  transverse  flows,  q^  and  q^,  arises 
liecause  the  At  chosen  for  stal)le  calculation  on  the  blocks  is  usually 
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too  large  for  stable  calculation  for  narrow  channels,  unless  the  coupling 
with  blocks  exists  only  in  respect  to  longitudinal  flow  from  the  channels 
to  blocks  at  end  points  of  such  channels. 

On  a given  side  of  a channel,  basically  four  physically  distinct  sit- 
uations can  occur:  (a)  a barrier  (levee)  or  block  ground  level  of  suffi- 
cient height  exists  to  prevent  lateral  flow;  (b)  overflow  exists  from  an 
adjacent  flooded  block  into  a channel  where  the  water  level  is  less  than 
the  adjacent  barrier  or  ground  level;  (c)  overflow  of  adjacent  barrier 
(levee)  exists  from  the  channel  to  an  adjacent  dry  block  or  one  where 
the  water  level  is  lower  than  the  barrier  elevation;  or  ^d)  both  the 
channel  water  level  and  the  water  level  on  the  adjacent  block  exceed 
the  height  of  any  intervening  barrier  and  the  lateral  flow  depends  on 
the  difference  of  water  level.  These  four  situations  are  illustrated 
in  Figure  9.  In  the  fourth  situation,  the  water  level  could  also  be 
lower  on  the  diannel  side  with  the  associated  lateral  flow  reversed. 

For  situation  (a)  there  is  no  problem,  the  appropriate  lateral  flow 
(q^  or  q^-)  being  constrained  to  zero  value.  For  situation  (c)  , the 
predictiVe-type  barrier  relation  (eq.  55) , with  auxiliary  relations 
(eqs.  39  and  40),  could  be  used.  In  principle,  the  above  predictive 
barrier  relations  should  apply  for  situation  (b)  as  well,  provided  that 
L in  equation  (39)  is  taken  as  the  channel  width  w.  However,  since 
w can  be  much  less  than  As  for  many  applications,  r can  be  so  small 
that  the  relation  for  q^  reduces  virtually  to  a diagnostic- type  rela- 
tion of  equation  (40),  or  more  specifically  of  equation  (9)  for  barrier 
overflow.  Since  situation  (b)  might  occur  on  one  side  of  the  channel 
and  situation  (c)  on  the  other,  and  since  both  should  be  evaluated  by 
relations  compatible  with  a common  time  level,  the  simple  diagnostic 
relation  (eq.  9)  has  been  adopted  for  both  situations  in  the  SURGE  II 
program.  This,  however,  still  demands  special  checks  and  possible 
adjustments,  as  will  be  discussed  later.  Finally  for  situation  (d) , a 
submerged  barrier-type  calculation  might  seem  appropriate  if  the  depth 
over  such  a barrier  is  small  compared  with  that  of  the  channel  or  adja- 
cent block;  however,  use  of  such  relations  in  preliminary  versions  of 
the  program  proved  to  be  very  vulnerable  to  numerical  instability.  The 
reason  for  this  is  related  to  the  above  discussion  concerning  the  usual 
case  where  w/As  is  very  small.  As  a consequence,  for  situations  of 
type  (d) , a special  calculation  is  required  which  treats  the  channel  as 
essentially  an  integral  part  of  the  associated  channel  block  or  the 
adjacent  block. 

As  stated  above,  for  overflow  situations  (b)  or  (c),  i,e.,  to  or 
from  the  channel,  the  relation, 

qn  = ^ CoDb(gDb)'4  , (57) 

is  used  where  = H - Zjj,  11  being  the  water  level  on  the  high  side  of 
the  barrier.  While  this  relation  gives  a valid  value  of  qrj(<lt  qf^ 
the  time  t,  the  value  of  q^  may  change  significantly  over  the  predic- 
tion interval  At  if  (gHj,)**  > w/At. 
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Figure  9 


Different  situations  along  a given  side  of  a channel 


Under  such  circumstances,  an  approximate  prediction  based  on  the 
initial  values  of  q^,  could  lead  to  physically  impossible  changes  of 
channel  level.  Thus,  tests  are  included  in  the  program  to  constrain  the 
lateral  flow,  such  that  q^.  - q^  alone  will  not  cause  the  channel  level, 
H^,  to  fall  below  a minimum  possible  value  nor  rise  above  a maximum 
possible  value,  depending  on  the  situation.  Six  different  situations 
requiring  tests  are  illustrated  in  Figure  10  (the  "mirror"  version  of 
each  is  also  a possible  situation).  Situations  where  one  side  of  the 
channel  is  blocked  are  special  cases  of  those  indicated.  For  situations 
A,  C,  and  E,  outflow  exceeds  inflow  and  the  horizontal  dashline  repre- 
sents a minimum  level  based  on  the  sill  depth  of  the  channel.  On  the 
other  hand,  for  situations  B,  D,  and  F,  the  horizontal  dashline  repre- 
sents a maximum  possible  level.  In  each  case,  the  maximum  possible  change 
in  is  indicated  as  AH^.. 


For  any  of  the  situations  illustrated  in  Figure  10,  the  SURGE  II  pro- 
gram compares  |q^  - qf|  with  |wHj,/At|.  If  the  latter  is  exceeded  by  the 
trial  value  of  |q^  - q^j  then  an  adjustment  is  made  in  q^  or  qj  such 
that  |q^^  - qfl  equals  [wAHc/At| . For  cases  A,  B,  C,  and  D,  both  q^.  and 
q£  are  prorated  by  a common  factor  to  satisfy  the  above  constraint.  For 
cases  E and  F,  only  the  overflow  q is  adjusted  to  be  consistent  with 
the  above  constraint. 

For  situation  (d)  where  the  channel  and  block  are  connected  by  a 
continuous  water  surface  (Fig.  10),  the  net  lateral  flow  to  the  channel, 
Aq,  is  taken  to  be  that  which  would  be  required  to  bring  HC  to  a value 
equal  to  the  existing  mean  level,  HM,  of  the  connected  channel  and 
block.  For  a channel  connected  to  a block  on  one  side  only  then, 

HM  = HB  . L + HC  - W (58) 

(L+W) 

where  UB  is  the  water  elevation  on  the  water-connected  block,  L is 
its  width,  while  HC  and  W are  the  water  elevation  and  width  for  the 
channel.  The  block  width  L is  AS  - W if  the  connected  block  is  the 
channel  block  containing  the  channel,  or  is  AS  for  an  adjacent  water- 
connected  block.  If  the  channel  is  water  connected  on  both  sides,  then 
the  above  relation  is  replaced  by  an  appropriate  average  over  both  blocks 
plus  the  channel. 

The  Aq  for  either  of  these  situations  is  taken  as 

Aq  = (HM  - HC)w/At  . (59) 

To  determine  the  individual  q^  and  q£  on  either  side  of  the  channel, 
the  mean  of  these  is  taken  to  be  that  which  is  calculated  as  the  flow 


between  blocks,  ignoring  the  presence  of  the  channel  (but  considering 
barriers).  Letting  this  be  denoted  q^,  then 

and  (60) 

qf  = q„,  - Aq/2  . 

This  system  of  calculation  leads  to  stable  results. 

d.  Channel  End-Point  Computations.  At  the  end  point  of  a given 
channel  system,  special  con^iutations  are  required.  Two  types  of  end 
conditions  are  used:  an  "H-end  condition"  is  used  where  a channel  dis- 
charges into  a lake,  bay,  or  sea,  in  which  ceise  the  channel  H value  at 
the  end  point  is  taken  equal  to  the  H of  the  adjacent  channel  block 
into  which  the  channel  discharges  (or  vice  versa);  a "Q-end  condition" 
is  used  at  the  head  of  a channel  or  river  at  which  point  the  discharge 
is  specified. 

For  a Q-end  point 


Q'  = ± 

H'  = (Q'  - B)/X  , (61) 

where  is  the  specified  river  discharge  (taken  as  zero  if  not  speci- 

fied); B equals  BP  or  -BN,  as  defined  by  equation  (51),  for  end 
points  occurring  at  the  positive  or  negative  end  of  the  channel  reach, 
respectively,  and  X is  as  defined  in  equation  (52).  The  sign  of  Q' 
is  taken  such  that  Q'  is  directed  into  the  channel,  depending  on  the 
channel-end  orientation.  There  are  four  possible  orientations  (see  App. 
B,  Fig.  B-3). 

The  H-end  points  also  have  four  possible  configurations;  these  are 
depicted  along  with  the  associated  adjacent  "ponding"  areas  (i.e.,  a 
block  with  Z < 0)  in  Figure  11.  For  an  H-end  point  neither  the  longi- 
tudinal flow  to  or  from  the  channel  nor  the  H at  the  junction  with  the 
ponding  block  is  specified  a It  is  required  only  that  the  pre- 

dicted H at  the  channel-end  point  and  that  of  the  ponding  block  be  the 
same.  Let  H*  be  the  (tentative)  predicted  H for  the  ponding  block 
in  the  absence  of  any  contribution  by  longitudinal  discharge  to  or  from 
the  channel  which  terminates  adjacent  to  that  block.  Thus,  H*  corre- 
sponds to  the  H resulting  from  the  routine  block  calculation  using 
equation  (33)  with  appropriate  adjustments  for  contained  channels  as 
might  occur  for  situations  3 and  4 shown  in  Figure  11.  These  adjustments 
are  discussed  in  a subsequent  subsection.  The  correct  predicted  H for 
the  ponding  block  in  the  presence  of  longitudinal  discharge  from  a 
channel  is  given  by 


H'  = H*  + (Q^  ^ Qd)At/2Aj,  , 


(62) 


configurations  and  index 
to  4) . 


where  and  are  the  new  and  previous  values,  respectively,  of 

the  discharge  from  channel  to  ponding  block,  and  is  the  effective 

surface  area  of  the  block.  For  situations  1 and  2 in  Figure  11, 

Ab  = (^s)2,  but  for  situations  3 and  4 a channel  might  exist  on  the 
ponding  block  in  which  case  Aj^  = (As  - w)As. 

Equation  (62)  involves  two  unknowns  H'  and  . However,  for  the 
channel , 


+ XH'  = B , (63) 

where  B = - BN  for  end-point  type  1 or  2 and  B = BP  for  end-point 
type  3 or  4,  BN,  BP  and  X being  those  quantities  defined  by  equations 
(46)  to  (52).  Note  that  for  end-point  type  1 or  2,  is  the  negative 

of  the  QC  value  for  the  channel. 

The  resulting  H'  and  for  an  "H-end"  condition  are 

H'  = (F  + BAt/2A^)/(l  + XAt/2Ab)  (64) 

Q'  = (B  - XF)/(1  + XAt/2Ab)  . (65) 

where 

F = H*  + Q^At/2Aj^  . (66) 

e.  Calculation  of  H on  Channel  Blocks.  For  blocks  with  D > 0 
and  containing  one  or  two  channel  reaches,  the  prediction  relation  for 
H given  by  equation  (33)  is  not  valid.  The  correct  relation  for  a 
channel  block  k having  location  i,j  is 

H’d.j)  = H(i,j)  + [U'(i,j)  - UCT’(k)]At/(As  - wx) 

+ [V’(i,j)  - VCT’(k)]At/(A8-wy)  (67) 

where  UCT  and  VCT  are  as  shown  in  Figure  6 and  correspond  to  the  q^ 
discussed  previously.  If  only  one  channel  exists  (i.e.,  if  wx  or  wy 
is  zero) , then 


UCT' (k)  = U' (i  + l,j)  if  wx  = 0 
or 

VCT'(k)  = V'(i,j  + 1)  if  wy  = 0 . 

IV.  APPLICATION  TO  THE  SABINE-CALCASIEU  SYSTEM 
1.  Adopted  Grid  and  Simulated  Topography. 

The  Sabine-Calcasieu  system  geographically  bridges  the  Texas-Louisiana 
border  and  is  physically  linked  by  a system  of  manmade  channels  and  a low- 
lying  region  extending  25  miles  between  Sabine  Lake  and  Lake  Calcasieu. 
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A local  chart  of  the  region  is  shown  in  Figure  12.  The  rectangular 
border  indicates  the  region  included  in  the  numerical  analog.  The 
selection  of  the  size  of  this  rectangle  is  dictated  by  the  basic  hydro- 
dynamic features  required  to  adequately  represent  the  region  and  then 
the  logistical  and  economic  limitations  placed  on  the  computations  by 
the  availability  of  computer  storage.  The  region  selected  is  56  x 40 
nautical  miles.  The  grid  size  (DELX)  is  taken  as  2 nautical  miles,  so 
that  IM  =28  and  JM  = 20. 

Figure  13  is  a contoured  plot  of  the  schematized  topography  super- 
imposed on  the  selected  grid  system.  The  offshore  topography  is  regular 
with  the  exception  of  a shallow  region  adjacent  to  Sabine  Pass  and  a 
slight  embayment  lying  between  Sabine  Pass  and  the  outlet  from  Lake 
Calcasieu  at  Cameron.  Both  lakes  are  adequately  represented  by  the  grid 
interval  tf  2 nautical  miles.  Figure  14  clearly  delineates  three  high 
topographic  areas  in  the  numerical  model:  the  Beaumont  rise  in  the 
northwest,  the  Orange  rise,  and  a more  gradual  rise  northeastward  to  the 
Lake  Charles  area.  The  low-lying  region  between  the  lakes,  immediately 
behind  the  shoreline  barrier,  and  forward  of  the  rises,  forms  a large 
ponding  area  during  the  inundation  sequences.  Between  each  rise  a major 
channel  is  present,  the  Neches  River,  the  Sabine  River,  and  in  the  Lake 
Charles  region,  the  Calcasieu  River  runs  northeastward  from  Lake  Calcasieu. 

The  deepest  block  in  the  system  is  -24  feet  fMSL) . Assuming  a 10- 
foot  surge,  a value  of  DELT  equal  to  or  less  than  260  seconds  (Sec.  Ill, 

1 , b)  is  required.  The  value  chosen  for  DELT  is  240  seconds. 

2 . Channel  and  Barrier  Schematization. 

The  numerical  discretization  of  the  area  shown  in  Figure  12  is  given 
as  an  overlay  in  Figure  15.  In  this  illustration  the  channel  network 
(shown  by  full  lines)  shows  the  landward  interconnection  of  Sabine  and 
Calcasieu  as  well  as  the  link  with  the  Intracoastal  Waterway  as  the 
lower  left-  and  right-hand  channels.  Each  channel  segment  has  been 
provided  the  physical  characteristics  of  width  and  cross-sectional  area 
that  best  reproduce  the  pertinent  information  for  the  channel  reach  that 
was  provided  by  the  Corps  of  Engineers.  The  extent  of  the  channel  system 
was  chosen  on  the  basis  of  past  inundation  history  and  the  judgment  of 
the  authors. 

Tlie  barrier  system,  also  shown  in  Figure  15,  represents  the  major 
manmade  and  natural  obstructions  to  flow  above  M5L.  At  the  shoreline 
the  major  dune  line  is  continuous  with  the  exception  of  an  apparent  open 
area  east  of  Sabine  Pass.  The  block  elevation  of  that  area  equals  the 
adjacent  barrier  heights.  Jetties  are  included  at  each  of  the  openings 
to  the  Gulf  of  Mexico.  Within  the  region  the  majority  of  barriers  are 
manmade  levees  erected  for  protection.  The  heights  of  all  barriers  were 
chosen  on  the  basis  of  data  provided  by  the  Corps  of  Engineers. 

Appendix  D has  a listing  of  all  data  used  for  the  Sabine-Calcasieu 
region  in  the  simulation  of  the  Hurricane  Carla  surge.  The  topography. 
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Figure  12.  Ma^D  of  Sabine-Calcasieu  region  showing  grid  boundary 


Figure  13.  Topography  contours  at  5- foot  intervals  for 

Sabine-Calcasieu  region  (broad  uncontoured  area 
between  Lakes  Sabine  and  Calcasieu  has  elevations 
between  0 and  5 feet) . 


barrier  data,  and  channel  data  are  the  same  for  the  astrotide  simulations 
and  the  standard  project  storms.  There  are  91  barrier  blocks  and  121 
channel  blocks  of  which  53  are  common  to  barrier  blocks.  Hxamples  of  null 
channel  blocks  are  for  K = 4,  6,  18,  21,  24,  etc.,  a total  of  19. 

Appendix  E shows  a plot  of  the  block  topography  with  the  channel  and 
barriers  sujierimposed.  This  plot  is  given  on  two  pages;  x (or  1)  runs 
from  left  to  right  and  y (or  J)  runs  from  bottom  to  the  top  of  the 
page  (I  values  are  indicated  along  the  top  of  both  pages  and  J along 
the  left  side  of  the  first  page).  Also  in  Appendix  E is  a listing  of 
the  key  arrays  for  channels  as  generated  by  the  program.  Note  that  the 
final  array  size  for  channels  is  128  (KCMP) , there  being  6 channels  which 
terminate  on  the  boundary  of  the  grid. 

As  an  illustration  of  barrier  input  note  from  Appendix  E that  for 
block  (2,2)  a y barrier  exists,  but  not  an  x barrier.  The  bed  eleva- 
tion of  block  (2,2)  is  -10  feet  while  that  of  block  (3,2)  is  -13  feet. 

Thus,  a value  of  ZX  of  -10  feet  should  have  been  input  for  this  block. 

The  listing  of  the  barrier  input  data  in  Appendix  D gives  the  information 
for  block  (2,2)  at  K = 12  with  ZX  = -100  (tenths  of  feet)  which  checks. 
The  actual  barrier  on  the  upper  side  indicates  a positive  6 feet.  How- 
ever, barrier  block  K = 13  at  the  adjacent  block  (3,2)  shows  a ZX 
value  of  -12  feet.  Reference  to  the  topography  in  Appendix  E indicates 
that  this  is  the  elevation  of  adjoining  block  (4,2)  which  is  higher  than 
block  (3,2)  and  hence  is  the  correct  entry. 

For  an  illustration  of  the  sign  coding  concerning  barriers  along 
channels,  refer  to  the  channel  input  data  in  Appendix  D and  the  plot  in 
Appendix  E.  Channel  block  K = 1 located  at  (8,1)  shows  a negative 
IWCX  and  a negative  IZCX  which  is  the  coding  for  double  levees  of 
equal  height  with  the  channel  in  between.  This  is  the  location  of  the 
double  jetty  entrance  channel  for  the  Sabine  region.  Channel  block  5 at 
location  (7,4)  shows  a (+,-)  signature  for  the  x channel  and  a (••■,  + ) 
signature  for  the  y channel.  Hence,  the  barrier  for  the  x channel  is 
on  the  inner  lateral  boundary  while  that  for  the  y channel  is  on  the 
outer  lateral  boundary  (see  App.  C,6).  Reference  to  Appendix  E key 
array  listings  shows  KCB  = 37  for  channel  block  5.  Barrier  block  37 
has  the  same  location  (7,4)  and  indicates  valid  barriers  of  a 5- foot 
elevation  above  MSL  for  both  the  x and  y channels. 

3.  River  Input  and  Hydrograph  Gage  Locations. 

There  are  three  river  discharge  locations  provided  for  tlie  Sabine- 
Calcasieu  region.  These  locations,  as  given  in  block  9 of  the  input 
(App.  D) , are  (28,15),  (4,19),  and  (14,19)  which  are  respectively  for 
the  Calcasieu  River  near  Lake  Charles,  the  Neches  River  north  of  Beau- 
mont, and  the  Sabine  River  north  of  Orange. 

Nine  gage  locations  for  the  astrotide  calibration  and  Hurricane 
Carla  simulation  are  shown  as  small  circles  in  Figure  16.  All  of  these 
with  the  exception  of  the  North  Sabine  Lake  gage  are  located  on  channels. 
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Figure  16.  Plot  of  channels  and  barriers  showing  locations  of 
hydrograph  gages  (0) . 


V.  TIDAL  CALIBRATION 


1.  TiJe  Data. 

Tlic  tidal  calibration  is  a required  step  in  the  preparation  of  the 
numerical  storm  surge  model.  Tliese  coin)utations  permit  the  adjustment 
of  the  parameters  representing  the  frictional  effects  in  the  channels 
and  low-lying  regions  of  tidal  inundation.  The  calibration  adjusts  the 
tidal  flows  in  order  to  adequately  predict  proper  phasing  and  tidal 
excursions  in  the  model  region.  Comparisons  are  made  with  actual  tide 
records  from  geographical  locations  corresponding  to  blocks  or  channels 
in  the  grid. 

Calibration  of  the  Sabine-Calcas ieu  region  was  carried  out  for  the 
springtide  conditions  that  existed  from  0000  hours,  22  August  to  0000 
hours,  27  August  1973.  Tide  recordings  at  nine  locations  in  the  region 
were  furnished  by  two  IJ.S.  Army  Engineer  Districts  (Fig.  16):  Sabine 
Pass  (southwest  jetty),  Port  Arthur,  north  Sabine  Lake,  Brakes  Bayou 
(Beaumont),  and  Orange,  Texas,  provided  by  the  Galveston  District;  Cameron, 
Hackberry,  Calcasieu  Lock,  and  Lake  Charles,  Louisiana,  provided  by  the 
New  Orleans  District. 

The  tidal  calibration  must  be  accomplished  during  a period  when  the 
tide  is  effectively  the  only  forcing  function  operating  on  the  system. 

Tills  requires  no  abnormal  riverflows  into  the  region  and  winds  whicli  will 
not  substantially  alter  the  slope  of  the  water  surfaces.  Such  conditions 
existed  for  the  first  96  hours  of  the  120-hour  record  period  and  this 
interval  was  used  in  tlie  tidal  calibration, 

2.  Estimation  of  fj-  for  Entrance  Channels. 

Many  of  the  bays  or  lagoons  along  the  Texas  coast  are  of  such  dimen- 
sions tliat  their  largest  natural  period  is  small  compared  with  the  tide 
period.  Moreover,  virtually  all  have  narrow  connections  with  the  Gulf 
of  Mexico.  These  two  features  conspire  to  produce  a reduction  of  tidal 
range  and  a significant  lag  within  the  bay  compared  with  the  gulf  tide. 

In  addition,  the  tidal  range  is  nearly  uniform  throughout  tlie  bay  except 
possibly  in  some  of  the  upper  reaches  of  adjoining  rivers.  For  these 
systems,  the  approximate  response  can  be  calculated  in  terms  of  the 
channel-friction  coefficient,  f^.,  (or  discharge  coefficient)  plus 
appropriate  dimensions  of  the  bay  and  entrance  channel  (Love,  1959). 

These  relations  can  be  used  to  get  at  least  a preliminary  estimate  of 
fj,  from  tlie  observed  response. 

Gonsider  a bay  of  total  M.SL  surface  area.  As,  which  is  connected  to 
the  sea  by  a channel  of  cross-sectional  area  Ac,  surface  width  W, 
effective  depth  Dc  (defined  as  Ac/W) , length  Lc,  and  channel-bed 
friction  coefficient  fp. 


Let  H be  the  volumetric  response  in  the  bay  at  time  t (where 

H • Aj;  represents  the  impounded  tidal  volume  above  MSL  at  time  t)  ; let 

Q be  the  tidal  flux  from  the  sea  to  the  bay.  Then, 

A,  f ■ Q . (68) 

Neglecting  the  inertia  effects  in  the  channel  for  the  slow  tidal  variation, 

the  slope  force  in  the  channel  is  balanced  by  friction  at  any  time  t; 

thus , 

Ilg  - H = miQlQ  , (69) 

where  Hg  is  the  given  tide  level  at  time  t outside  the  bay  entrance 
and  m is  a dimensional  constant  for  the  system  given  by 


m 


fL 


8°  aJ 
c c 


(70) 


This  can  also  be  written  in  the  form. 


8(Cd-A^)2  ' 

where  C^j  is  the  discharge  coefficient  characterizing  the  constricted 
opening  between  bay  and  sea. 

Assuming  the  input  tide  Hg  is  simple  harmonic  with  period  T and 
an^litude  a©  then. 


Hg  = a^  cos  ut  , (71) 

where  u = 2tr/T.  Ignoring  the  second-order  compound  tide  due  to  non- 
linearity in  equation  (69),  the  response  will  be  roughly  of  the  form, 

H = rap  cos  ((t)t-()i) , (72) 

where  is  a phase  lag  and  r is  the  relative  amplitude  response.  If 
these  are  substituted  into  equations  (68)  and  (69)  and  the  quantity  1q|Q 
expanded  in  the  Fourier  series  form,  it  can  be  shown  that 


and 

r = cos  i{) 

(73) 

, , / 1 + - 1 
sin  4>  •=  g 

(74) 

where 

B - — m(A^(u)2a^ 

(75) 

56 


(Love,  1959).  A plot  of  r and  ((>  versus  the  dimensionless  parameter 
B is  shown  in  Figure  17.  The  timelag  of  the  high  tide  in  the  bay  rela- 
tive to  that  outside  the  bay  is  simply  T = ifi/Zir,  for  (ji  in  radians 
(or  T = ()i/360  for  ()>  in  degrees) . 

Thus,  if  r or  is  estimated  from  observations  it  is  possible  to 
get  an  estimate  of  B.  Generally,  the  value  obtained  from  the  observed 
r will  differ  from  that  obtained  from  the  observed  in  this  event 

an  average  of  the  values  of  B can  be  used  to  estimate  In  terms 

of  B,  fj,  is  given  by 

f gB 

C 8 „2a?-^2 

s 

where 

L 

a2  = 

D a2  ’ 
c c 

It  is  emphasized  that  the  above  analysis  pertains  to  a bay  system 
connected  to  the  sea  by  a single  channel  of  uniform  dimensions.  The 
results  can  be  generalized  for  the  case  of  a series  of  N channels  of 
different  dimensions  or  of  N channels  in  parallel  or  combination  of 
both  (as  in  the  Sabine-Calcasieu  system)  by  using  an  effective  value 
of  a^. 

Let  designate  the  value  of  for  an  individual  channel  as 

evaluated  by  equation  (77),  Then,  the  effective  value  of  a2  for  a 
series  of  N channels  is  simply 

“s  = I “n  • (78) 

n=l 


(76) 

(77) 


However,  for  N diannels  in  parallel  the  effective  is  given  by 


a 


2 

P 


N 

I 


(79) 


For  a series  containing  a parallel  subset,  the  effective  a^,  for 
the  latter  is  used  in  equation  (78).  If  two  or  more  complex  entrance 
channels  are  in  parallel  then  the  effective  values  of  a are  used  in 
place  of  Op  in  equation  (79). 

The  use  of  this  procedure  will  be  illustrated  for  the  Sabine-Calcasieu 
system.  In  the  numerical  simulation  scheme  there  is  a total  of  40  blocks 
of  2 X 2 nautical  miles  covered  with  water  and  in  communication  with  the 
sea.  This  represents  a surface  area  of  5.91  x 10^  square  feet.  In  addi- 
tion, the  channels  contribute  a total  of  0.64  x lo^  square  feet.  Thus, 
the  total  surface  area  for  the  combined  system  is  As  = 6.55  x lo^  square 
feet  (3.77  x 10®  square  feet  for  the  Sabine  part  and  2.78  x 10®  square 
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feet  for  the  Calcasieu  part).  The  two  parts  of  the  system  are  coupled 
via  the  Intracoastal  Waterway  and  their  responses  are  about  the  same, 
so  the  combined  system  is  treated  as  one, 

A summary  of  data  and  calculations  pertinent  to  the  entrance  channels 
for  the  Sabine-Calcasieu  system  is  given  in  Table  2 (see  also  Fig.  15  and 
App.  D) . The  simulated  Sabine  Pass  between  the  gulf  and  Lake  Sabine  con- 
sists of  two  sections  (1  and  2 in  Table  2)  of  different  dimensions  in 
series.  However,  Calcasieu  Pass  consists  of  a pair  of  parallel  channels 
(4  and  5 in  Table  2)  in  series  with  a simple  channel  (3  in  Table  2) . The 
individual  for  each  channel  is  also  shown  in  Table  2.  The  effective 

for  Sabine  Pass  is  the  first  partial  sum  shown  in  the  last  column.  The 
effective  value  of  for  the  parallel  part  of  Calcasieu  Pass  is  shown 

in  the  last  column,  opposite  entries  4 and  5.  The  effective  value  for 
Calcasieu  Pass  is  the  partial  sum  indicated  in  the  last  column.  The 
effective  value  for  the  entire  pass  system  is  evaluated  from  the  Sabine 
Pass  and  Calcasieu  Pass  values,  using  equation  (79)  for  parallel  systems: 

=0.32  X 10”^  (square  feet)"^  . 


Table  2.  Data  on  simulated  Sabine  Pass  and  Calcasieu  Pass. 


n 

Wc 

Dc 

^c 

Lc 

X 10® 

a2  X 10^ 

(ft) 

(ft) 

(ft2) 

(ft) 

(ft- 2) 

(ft- 2) 

Sabine  Pass 

1 

2,330 

20 

46,600 

24 , 360 

0.561 

0.561 

2 

2,860 

21 

60,060 

36,480 

0.482 

0.482 

Subtotal 

1.043 

Calcasieu  Pass 

3 

800 

32 

25,600 

24,360 

1.162 

1. 162 

4 

500 

40 

20,000 

12,160 

0.760 

O.455I 

5 

1,000 

16 

16,000 

34,480 

8.960 

Subtotal 

1,617 

^Evaluated  by  parallel  channel  relation. 


The  observed  ranges  and  times  of  minimum  tide  for  25  August  1973  for 
the  Sabine-Calcasieu  system  are  given  in  Table  3.  Gage  1 is  used  as  the 
input  gulf  tide.  The  average  of  all  other  gages  is  used  as  the  response. 
The  indicated  amplitude  response  is 


r 


1.50 

2.59 


0.58  . 


Using  a tidal  period  of  25  hours  the  indicated  phase  lag  is 


♦ 


(20.8  - 17.5) 


59 


360 

25 


47*  . 


» 


Table  3.  Ranges  and  times  (c.d.tj  of  available  observed 
tides  in  the  Sabine-Calcasieu  system  for  25 
August  1973. 


Gage  No. 


7 

8 
9 


Place 

Range 

(ft) 

Time 

(hr) 

Sabine  Pass,  southwest  jetty 

2.59 

17.5 

Port  Arthur 

1.53 

North  Sabine  Lake 

1.40 

Beaumont 

1.52 

Orange 

1.40 

23.0 

Cameron 

2.05 

17.5 

Hackberry 

1.06 

22.0 

Calcasieu  Lock,  west 

1.45 

20.5 

Lake  Charles 

1.60 

21.5 

2 to  9,  inclusive 

1.50 

20.8 

From  Figure  17  the  corresponding  values  of  B are  4.7  and  3.6,  respec- 
tively, with  an  average  of  4.1.  The  tidal  frequency  is 


(i)  = 25~x^^J~eob~  ^ ^ 10~^  radians  per  second 

and  Uq  = 2.57/2  or  1.3  feet.  Consequently,  the  estimated  f^,  for  the 
entrance  channels  is  from  equation  (76):  f^,  = 0.0018. 

The  final  selected  value  of  fc  for  the  entrance  channels  is  0.0015 
as  determined  by  trial  rtins.  This  is  somewhat  less  than  the  above  esti- 
mate. The  difference  might  be  accounted  for  by  the  fact  that  the  tidal 
hydrograph  is  not  really  simple  harmonic  but  contains  compound  tides  (of 
higher  frequency)  giving  the  sharp  minimum  and  broad  or  double-peaked 
maxima.  The  effective  frequency  is  consequently  somewhat  greater  than 
the  w given  above,  thus  yielding  a smaller  closer  to  0.0015. 

3.  Final  Calibration  for  Tide. 

The  major  control  on  the  response  of  the  bay  to  the  tides  are  the 
dimensions  and  friction  factor  for  the  entrance  channels  as  discussed 
above.  In  this  connection,  it  should  be  pointed  out  that  channel  dimen- 
sions (width  and  depth)  were  taken  such  that  the  average  cross-sectional 
area  (under  MSL  conditions)  for  a given  reach  is  represented  by  the  pro- 
duct of  these  dimensions.  Thus,  if  the  depth  is  taken  as  the  mean  for 
the  reach,  then  the  width  will  be  somewhere  between  the  width  of  the 
dredged  channel  and  the  surface  width  of  the  natural  channel. 

The  values  of  channel  friction  for  the  remaining  channels  and  of  the 
block  friction  were  selected  by  a tri al-and-error  procedure,  starting 
with  a uniform  value  throughout.  The  final  values  of  channel  friction 
for  the  upper  reaches  of  the  Neches  and  Sabine  Rivers  were  taken  as 
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0.0025  to  give  a reasonable  agreement  for  the  Beaumont  and  Orange  tide 
response;  it  was  necessary  to  use  a low  value  (0.0005)  for  the  upper 
reach  of  the  Calcasieu  River  to  reproduce  the  Lake  Charles  tidal  hydro- 
graph. The  latter  three  gages  (Beaumont,  Orange,  and  Lake  Charles)  have 
connections  to  the  inner  bay  areas  only  via  channels,  hence  their  re- 
sponses are  fairly  sensitive  to  the  channel  friction.  The  low  value  for 
the  Calcasieu  River  may  be  due  to  underestimates  of  the  effective  channel 
widths,  which  would  demand  a less  than  normal  friction  factor. 

The  block  friction  for  the  tide  calculations  was  taken  as  0.0015  to 
got  a reasonable  agreement  for  the  north  Sabine  Lake  gage.  However, 
later  calculations  for  the  Hurricane  Carla  simulation  (which  is  more 
sensitive  to  block  friction  than  the  astrotide)  indicated  that  0.0025 
(as  used  in  the  Galveston  Bay  simulations)  was  more  appropriate. 

The  results  of  the  final  astronomical  tide  simulations  for  a 96-hour 
period  starting  0000  hours  c.d.t.,  22  August  1973,  are  given  in  Figures 
18  to  26,  and  Appendix  F.  The  input  tide  (Fig.  18)  corresponds  to  the 
observed  tide  for  the  period  at  Sabine  Pass  (southwest  jetty).  In  the 
subsequent  eight  figures  the  computed  (full  line)  and  observed  (line 
with  circles)  are  compared  for  the  eight  different  gages  within  the 
system;  the  gages  are  identified  in  the  figures.  Note  that  the  observed 
values  for  each  gage  have  been  adjusted  with  respect  to  a local  datura, 
taken  as  the  gage  mean  for  a 120-hour  period  starting  0000  hours  c.d.t., 

22  August  1973.  In  all  cases,  the  computed  ranges  are  in  fairly  good 
agreement  with  the  observed;  however,  there  seems  to  be  a consistent 
tendency  for  the  computed  to  lag  the  observed.  This  might  be  due  to  a 
possible  time-shift  error  for  the  input  gage.  Although  a lowering  of  the 
frictional  coefficient  for  the  entrance  channel  would  decrease  the  lag 
within  the  system,  it  would  also  increase  the  range  of  the  tide  every- 
where in  the  system.  It  was  felt  that  it  was  more  important  to  reproduce 
the  range  than  the  times  of  high  and  low  water,  and  hence  the  value  of 
f(.  = 0.0015  for  the  entrance  channels  was  retained. 

For  the  upper  Calcasieu  River  (Figs.  25  and  26)  the  computed  water 
level  (which  refers  to  a common  MSL  dat  un  for  the  system)  and  the  observed 
water  level  display  an  apparent  vertical  shift.  This  could  be  related  to 
possible  wind  effects  in  the  second  part  of  the  record,  which  have  been 
ignored  in  the  computations. 

The  steady  river  discharges  adopted  in  the  astrotide  runs  were  800 
cubic  feet  per  second  for  Calcasieu  River,  1,100  cubic  feet  per  second 
for  Neches  River,  and  1,500  cubic  feet  per  second  for  Sabine  River. 

Serial  listings  of  the  computed  water  levels  at  the  gages  discussed 
above  are  given  in  Appendix  F,  along  with  listings  of  volume  transport 
at  six  channel  positions.  Flow  at  points  1 and  2 correspond  to  input 
(if  positive)  to  the  system  through  Sabine  Pass  and  Calcasieu  Pass, 
respectively.  Since  the  tide  amplitude  is  less  than  the  seaward  barriers, 
the  two  passes  represent  the  only  source  of  water  for  normal  conditions. 
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Figure  18.  Astronomical  tidal  hydrograph  for  Sabine  Pass,  southwest  jetty 
(input  for  tide  calibration). 
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Figure  19.  Astronomical  tide  for  Port  Arthur  corresponding  to  input  of 
Figure  18. 
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Figure  20.  Astronomical  tide  for  north  Sabine  Lake. 
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Figure  21.  Astronomical  tide  for  Beaumont,  Neches  River,  and  Brakes  Bayou. 
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Figure  22.  Astronomical  tide  for  Orange  Naval  Station,  Sabine  River. 
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Figure  23.  Astronomical  tide  for  Cameron,  Calcasieu  Pass. 


Figure  24.  Astronomical  tide  for  Hackberry,  Calcasieu  River  and  Pass. 
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Figure  25.  Astronomical  tide  for  Intracoastal  Waterway  at  Calcasieu  Lock,  west 
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ronomical  tide  for  Lake  Charles,  Calcasieu  River. 


Reproductions  of  channel  output  at  three  different  times  (30,  60,  and 
90  hours  from  start)  are  shown  in  Appendix  F.  The  output  shows  flows 
(in  cubic  feet  per  second),  direction  of  flow,  and  water  level  along  the 
various  channel  reaches  at  the  specified  times. 

VI.  HURRICANE  CARLA  VERIFICATION 


1 . Forcing  Function  Input. 


a.  Wind-Stress  Fields.  The  x and  y components  of  the  wind  stress 
for  each  3 hours  in  a 72-hour  period  for  an  8 by  6 coarse  grid  for  Hurri- 
cane Carla  are  given  in  the  input  listings  in  Appendix  D.  For  convenience 
in  spotting  possible  errors  in  input,  the  wind-stress  vectors  were  plotted, 
based  on  the  above  input,  by  a special  subprogram.  Samples  of  these  plots 
for  each  12  hours,  are  shown  in  Figures  27  to  32.  The  plots  showed  sus- 
pect entries,  which  were  subsequently  corrected  before  any  runs  were 
attempted,  and  have  I increasing  upward  and  J increasing  to  the  left; 

1. e.,  the  seaward  boundary  is  on  the  right. 

b.  River  Discharge  Input.  The  river  discharges  for  the  Calcasieu 
River,  Neches  River,  and  Sabine  River  for  each  3 hours  are  listed  as 
block  (IDENT)  12  in  Appendix  D. 

c.  Gulf  Hydrograph  Input.  The  final  input  for  HG,  the  water  level 
input  along  the  seaward  boundary,  was  taken  as  interpolated  values  between 
Sabine  Pass  and  Calcasieu  Pass  with  input  sequences  at  those  passes  adjust- 
ed to  match  the  observed  values  at  the  Sabine  Pass  U.S.  Coast  Guard  Station 
and  Cameron  after  some  modification  due  to  flow  through  these  passes.  The 
input  is  given  sequentially  at  5-hour  intervals  along  with  the  wind-field 
input  in  Appendix  D. 

2.  Further  Adjustments  and  Results. 

a.  Adjustments.  In  the  series  of  runs  for  the  Hurricane  Carla  simula- 
tion, it  was  necessary  to  make  some  adjustments  in  the  block  topography, 
particularly  in  the  upper  reaches  of  the  Neches  River,  in  order  to  provide 
more  ponding  area  at  the  levels  of  flooding  encountered.  These  changes, 
which  are  reflected  in  the  final  topography  (App.  D) , do  not  change 
the  results  of  the  astronomical  tide  calibration  because  the  changes  were 
at  levels  well  above  those  encountered  with  the  astrotide  runs. 

A further  modification  was  the  reduction  of  the  wind-stress  values  to 
80  percent  of  those  shown  in  the  listings  and  in  the  vector  plots  for  the 
u?iper  left-hand  region  of  the  grid.  Specifically  for  I.LE.3  and  J.GE.4, 
the  wind-stress  components  were  so  reduced  in  the  final  runs  for  Hurricane 
Carla.  Tliis  reduction  was  also  used  in  the  later  application  for  Standard 
Project  Hurricane  (SPH)  simulations.  The  rationale  for  this  adjustment  is 
based  on  the  greater  sheltering  in  this  region  due  to  both  topography  and 
vegetation.  The  initial  H for  all  locations  in  the  bay  was  taken  as 
3.2  feet. 


71 


Figure  27.  Wind-stre.ss  vectors  for  Hurricane  Carla,  over 
Sabine-Calcasieu  region  on  an  8-nautical  mile 
grid;  time  = 12  hours. 


Figure  29.  Wind-stress  vectors  for  Hurricane  Carla,  over 
Sabine-Calcasieu  region  on  an  8-nautical  mile 
grid;  time  = 36  hours. 
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Figure  32.  Wind-stress  vectors  for  Hurricane  Carla,  over 
Sabine-Calcasieu  region  on  an  8-nautical  mile 
grid;  time  = 60  hours. 


b.  Results ■ The  results  of  the  Hurricane  Carla  simulation  are  given 
in  Figures  33  to  41,  and  Appendix  G.  The  input  (observed)  hydrograph  for 
Sabine  Pass  is  shown  in  Figure  33  for  a 72-hour  period  starting  at  0000 
hours  c.s.t.,  10  September  1961.  These  results  are  based  on  a block 
friction  factor  of  0.0010. 

The  computed  and  observed  values  (where  available)  at  gages  2 to  9 
are  shown  in  Figures  34  to  41.  The  principal  discrepancy  occurs  at  Beau- 
mont where  the  computed  peak  surge  exceeds  the  peak  observed  value  by 
about  0.8  foot.  It  was  found  later  that  by  increasing  the  block  friction 
to  0.0025,  this  difference  was  reduced  to  0.4  foot  without  materially 
changing  the  results  at  other  key  locations  in  the  system. 

The  auxiliary  sample  output  for  the  simulated  Hurricane  Carla  run 
(App.  G)  gives,  in  addition  to  the  serial  listings  of  the  above  hydro- 
graphs and  flow  at  the  two  main  passes,  sample  listings  of  channel  out- 
put at  elapsed  times  of  30  and  60  hours. 


VII.  STANDARD  PRaJECT  HURRICANE  (SPH) 

1 . LR-ST  Storm  Data. 

The  large  radius,  slow  translation  (LR-ST)  storm  was  utilized  as  an 
atmospheric  forcing  function  for  the  verified  model  of  the  Sabine-Calcasieu 
system.  The  storm  parameters  were  extracted  from  the  pertinent  gulf  coast 
section  of  the  National  Hurricane  Research  Project  Report  No.  33  (Graham 
and  Nunn,  1959).  Table  4 lists  these  values  which  were  also  used  in  con- 
junction with  the  analytic  storm  representation  given  by  Jelesnianski 
(1965) . 


Table  4.  Atmospheric  parameters  for  the  large  radius,  slow 

translation  (ST)  and  medium  translation  (MT)  storms. 


Parameters 

ST  storm 

MT  storm 

Radius  to  maximum  winds 

27  nmi 

27  nmi 

Maximum  windspeed 

100  mi/h 

100  mi/h 

Central  pressure 

27.55  in 

27.55  in 

Translation  speed 

4 kn 

11  kn 

Wind-stress  vector  plots  have  been  prepared  beginning  at  t = 30  hours 
and  at  10-hour  increments  to  t = 80  hours  (Figs.  42  to  47).  The  storm 
track,  which  is  taken  normal  to  the  general  shoreline,  has  the  Sabine- 
Calcasieu  system  on  the  right-hand  side  of  the  storm  approaching  the 
coastline.  Landfall  of  the  storm  center  is  close  to  grid  block  1,1. 

The  orientation  of  these  plots  relative  to  the  topography  is  similar  to 
the  wind  fields  shown  for  the  Hurricane  Carla  verification.  The  gulf 
hydrographic  input,  provided  by  the  Galveston  District,  was  developed 
by  an  application  of  a one-dimensional  bathystrophic  mode]  (Marinos  and 
Woodward,  1968;  Bodine,  1971).  A tidal  component  has  been  added  to  this 
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Figure  33.  Hydrograph  at  Sabine  Pass,  southwest  jetty  for 
Hurricane  Carla. 
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Figure  34.  Ilydrographs  at  Sabine  Pass,  U.S.  Coast  Guard  Station  for 
Hurricane  Carla  (FK  = 0.0010). 
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Figure  36.  Hydrograph  at  north  Sabine  Lake  for 
Hurricane  Carla  (FK  = 0.0010). 
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Figure  37.  Ilydrographs  at  Beaumont,  Neches  River,  and  Brakes  Bayou  for 
Hurricane  Carla  (FK  = 0.0010). 
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Figure  38.  Ilydrographs  at  Orange  Naval  Station,  Sabine  River  for 
Hurricane  Carla  (FK  » 0.0010). 
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rigure  39.  llydrograph  at  west  end  of  Intracoastal  Waterway  for 
Hurricane  Carla  (FK  = O.OOIQ), 
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Figure  40.  Hydrographs  at  Cameron,  Calcasieu  Pass  for  Hurricane  Carla 
(FK  = 0.0010) . 
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Figure  44.  liiiid-stress  vectors  for  SPII  large  radius,  slow 
translation  (LR-STl  on  :in  8-nautical  mile  grid; 
t ' me  = 50  hours . 
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Figure  45.  Wind-stress  vectors  for  SPH  large  radius,  slow 
translation  (LR-ST)  on  an  8-nautical  mile  grid; 
time  = 60  hours. 
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Figure  47.  Wind-stress  vectors  for  SPH  large  radius,  slow 
translation  (LR-ST)  on  an  8-nautical  mile  grid; 
time  = 80  hours. 
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open-coast  surge.  A total  rainfall  of  16  inches  in  24  hours  is  included 
as  input.  The  results  are  based  on  a block  friction  factor  of  0.0010 
and  arc  therefore  tentative. 

2.  LR-ST  Storm  Results. 

Nine  simulated  hydrograplis  are  presented  in  Figures  48  to  56.  The 
hydrograph  locations  are  selected  to  coincide  with  the  locations  of  the 
Hurricane  Carla  graphs  and  are  shown  on  the  base  map  in  Figure  16.  The 
maximum  water  level  excursion  in  these  graphs  is  nearly  19  feet  and 
occurs  at  the  Beaumont  gage  (Fig.  52)  at  approximately  t = 80  hours. 

The  highest  elevation  at  Sabine  Pass,  southwest  jetty  (Fig.  48)  is  13 
feet  and  occurs  at  t = 77  hours.  The  Port  Arthur  surge  crests  (Fig.  52) 
at  slightly  greater  than  14  feet  shortly  after  high  water  is  reached  on 
the  open  coast.  The  gage  at  north  Sabine  Lake  (Fig.  51)  reaches  a maxi- 
mum of  15.3  feet  which  coincides  with  Port  Arthur.  The  surge  at  Beau- 
mont develops  continuously  from  the  7- foot  level  at  t = 66  hours  to  a 
maximum  at  80  hours  in  direct  correlation  with  the  surge  character  at 
the  coastline.  However,  the  7-foot  level  is  reached  at  the  coast  approx- 
imately 6 hours  before  Beaumont.  At  the  Orange  Naval  Station  (Fig.  53) 
the  surge  development  appears  more  monotonic  with  a steady  climb  from 
t = 54  hours  to  the  peak  surge  of  nearly  16  feet  shortly  before  the 
Beaumont  peak.  The  recession  stage  at  Sabine  Pass  occurs  quickly  with 
passage  of  the  storm,  and  the  open-coast  water  level  has  returned  to 
normal  level  by  t = 90  hours. 

Drainage  inland  slowly  reduces  water  levels  and,  at  t = 90  hours. 

Port  Arthur  and  north  Sabine  Lake  have  water  elevations  of  approximately 
10  feet.  Farther  inland,  the  water  elevation  at  Beaumont  and  Orange 
stands  at  13  feet.  The  system  continues  to  drain  over  the  next  10  hours 
of  prototype  computation  but  slows  considerably  since  the  runoff  roaches 
peak  rate  at  the  end  of  the  run. 

Recalculations  at  the  Galveston  District,  using  a block  friction  of 
0.0025  which  improved  the  Hurricane  Carla  simulation,  indicated  signifi- 
cant reductions  in  the  peak  surges  for  the  LR-ST  storm  (Table  5).  The 
greatest  reduction  is  for  Beaumont,  as  in  the  case  of  Hurricane  Carla; 
however,  the  amount  of  reduction  is  disproportionately  greater  than 
that  seen  for  Hurricane  Carla. 

Table  5.  Comparison  of  peak  surges  for  the  LR-ST  storm, 
using  two  different  block  friction  factors. 

Location  Surges  (ft  above  MSL) 

„ FK  = 0.0010  FK  = 0.0025 

Port  Arthur  14.3  14.1 

Beaumont 

Orange  Naval  Station 
Lake  Charles 


18.7 

15.9 

14.1 
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Figure  50.  Hydrograph  for  SPH,  LR-ST  at  Port  Arthur  (FK  = 
0.0010). 
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3.  LR-NTT  Storm  Data. 

The  large  radius,  mcdiuin  translation  (LR-MT)  storm  has  identical 
characteristics  to  the  LR-ST  storm  with  the  exception  of  a higher  trans- 
lation speed  of  11  knots.  Wind  vector  plots  from  t = 15  hours  to  t = 40 
hours  are  shown  at  5-hour  increments  in  Figures  57  to  62.  The  storm 
track  is  identical  to  that  of  the  LR-ST  storm.  The  gulf  hydrographic 
input  was  derived  by  one-dimensional,  bathystrophic  analysis  and  pro- 
vided by  the  Galveston  District.  Runs  were  made  both  with  and  without 
rainfall.  Again,  the  results  given  graphically  below  are  the  tentative 
results  based  on  FK  = 0.0010. 

4.  LR-NfT  Storm  Results. 


The  more  rapid  movement  of  the  storm  center  across  the  Sabine- 
Calcasieu  system  yielded  generally  smaller  water  level  excursions  in- 
side the  bay  system  in  comparison  with  the  LR-ST  storm.  Hydrographs  at 
the  established  prototype  locations  are  shown  in  Figures  63  to  71  for 
the  computer  run  with  rainfall  (16  inches)  and  without  rainfall.  Note 
that  direct  comparison  between  the  LR-ST  results  and  LR-MT  results 
should  be  made  on  the  basis  of  Figures  48  to  56  and  63  to  71,  respec- 
tively. All  of  the  SPH  runs  use  an  initial  water  level  of  about  2.5 
feet  in  the  bay  system. 

A summary  of  the  peak  values  and  relative  times  of  water  level  at 
seven  locations  for  the  three  different  SPH  runs  is  given  in  Table  6. 
Although  the  absolute  values  of  the  water  levels  depend  on  the  value 
of  FK  (as  discussed  in  previous  sections),  all  results  in  Table  6 are 
based  on  the  same  FK  and  hence  the  difference  between  values  is  not 
too  sensitive  to  FK. 


Table  6.  Coicparison  of  peak  surge  and  time  of  peak  surge,  showing  effects  of  translational 
speed  of  storm  and  rainfall  (FE  > 0.0010  for  all  three  cases). 


Location 

Slow  speed 

1 Medium  speed 

With  rainfall 

1 With  rainfall  I 

1 Without  rainfall 

(ft  above  MSL) 

(ft  above  MSL) 

(time) 

(ft  above  MSL) 

(time) 

Sabine  Pass  entrance 

13.0 

0 

14.9 

0 

14.9 

0 

Port  Arthur 

14.3 

2 

13.2 

2 

12. S 

2 

North  Sabine  Lake 

15.3 

1 

IS. 3 

1 

14.7 

1 

Beaumont 

IS.  7 

4 

IS.l 

S 

11. S 

6 

Orange  Naval  Station 

IS. 9 

4 

14. S 

S 

11.7 

6 

Cameron 

11.3 

1 

11.0 

1 

10.8 

1 

Lale  Charles 

14.1 

6 

14.2 

6 

13.2 

6 

Lsearest  hour  after  that  of  Sabine  Pass  entrance. 


Comparison  of  the  first  and  second  sets  of  peak  levels  in  Table  6 
indicates  a reduced  response  at  nearly  all  stations  within  the  Sabine- 
Calcasieu  system  with  an  increase  in  the  translational  speed  of  the 
storm,  in  spite  of  the  increased  surge  at  the  shoreline  (Sabine  Pass 
entrance).  A reduction  in  volume  response  within  the  system  is  expected 
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Figure  57.  Wind-stress  vectors  for  SPII  large  radius,  medium 
translation  (LIl-MT)  on  an  8-nautical  mile  grid; 
time  = 15  hours. 


Figure  58.  Wind-stress  vectors  for  SPH  large  radius,  medium 


translation  (LR-OT)  on  an  8-nautical  mile  grid; 
time  = 20  hours. 


Figure  60.  Wind-stress  vectors  for  SPll  large  radius,  medium 
translation  (LR-MF)  on  an  8-nautical  mile  grid; 
time  = 30  hours. 
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Figure  61.  Wind-stress  vectors  for  SPH  large  radius,  medium 
translation  (LR-MT)  on  an  8-nautical  mile  grid; 
time  = 35  hours. 
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Figure  62.  Wind-stress  vectors  for  SPH  large  radius,  medium 
translation  (LR-f-fT)  on  an  8-nautical  mile  grid; 
time  = 40  hours. 
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Figure  63.  Hydrographs  for  SPH,  LR-^T^  (with  and  without 
rainfall)  at  Sabine  Pass,  southwest  jetty. 


Figure  64.  Hydrographs  for  SPH,  LR-HT  (with  and  without 
rainfall)  at  Sabine  Pass,  U.S.  Coast  Guard 
Station  (FK  = 0.0010). 
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Figure  65. 


Hydrographs  for  SI’H,  LR-fTF  (with  and  without 
rainfall]  at  Port  Arthur  (FK  = 0.0010). 


Figure  66. 


Hydrographs  for  SPH,  LR-ffT  (with  and  without 
rainfall)  at  north  Sabine  Lake  (FK  = 0.0010). 
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Figure  67.  Hydrographs  for  SPH,  LR-MT  (with  and  without 

rainfall)  at  Beaumont,  Neches  River,  and  Brakes 
Bayou  (FK  = 0.0010). 
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Figure  68.  Ilydrographs  for  SPH,  LR-MT  (with  and  without 

rainfall)  at  Orange  Naval  Station,  Sabine  River 
(FK  » 0.0010). 
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Figure  71.  Ilydrographs  for  SPH,  LR-MT  (with  and  without 


rainfall)  at  Lake  Charles,  Calcasieu  River 
(FK  = 0.0010). 
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for  the  greater  speed  (shorter  duration)  storm  because  of  the  constricted 
connection  to  the  sea.  Port  Arthur  shows  a reduction  of  1.1  feet  for  the 
NfT  storm  relative  to  the  ST  storm;  north  Sabine  Lake  appears  to  show  no 
change.  An  examination  of  the  wind  fields  close  to  the  time  of  the  peak 
surges  (Figs.  46  and  61)  indicates  that  a greater  wind-induced  setup 
within  the  lake  occurs  between  Port  Arthur  and  the  north  Sabine  Lake 
station  for  the  medium  speed  storm,  due  to  the  favorable  orientation  of 
the  winds  near  the  time  of  peak  surge  at  the  lake  entrance. 

The  response  at  Beaumont  and  Orange,  both  of  which  are  well  inland 
of  the  main  lake  area,  shows  a significant  reduction  (3.6  and  1.4  feet, 
respectively)  as  well  as  a greater  timelag  for  the  faster  storm.  More- 
over, the  peak  elevations  for  both  of  these  stations  are  somewhat  less 
than  that  at  north  Sabine  Lake  for  the  MT  storm  in  contrast  to  the  situa- 
tion for  the  ST  storm.  The  limited  access  of  water  to  these  regions  is 
apparently  responsible  for  this  sensitivity  to  storm  duration. 

The  influence  of  rainfall  and  associated  runoff  from  drainage  areas 
well  inland  is  shown  very  dramatically  from  a comparison  of  the  second 
and  third  sets  of  peak  levels  in  Table  6,  particularly  for  Beaumont  and 
Orcinge  Naval  Station,  where  runoff  produces  a differential  flooding  of 
3.6  and  2.8  feet,  respectively.  A differential  of  about  0.6  foot  due  to 
runoff  and  rainfall  occurs  even  within  Lake  Sabine.  The  effects  within 
Lake  Calcasieu  and  upstream  to  the  northeast  are  less  pronounced  due  to 
the  smaller  runoff. 


VIII.  CONCLUSION 

The  use  of  a modified  program  for  inclusion  of  subgrid  scale  channels 
has  been  demonstrated  to  be  essential  for  simulation  of  tides  in  the 
upper  reaches  of  a system  like  the  Sabine-Calcasieu  region,  where  the 
primary  connection  to  locations  such  as  Beaumont,  Orange,  and  Lake  Charles 
is  via  river  channels  which  would  not  otherwise  be  resolved  by  a grid 
scheme  of  the  order  of  a 1-nautical  mile  scale.  Even  for  conditions  of 
extreme  flooding,  as  occur  during  hurricanes,  the  incorporation  of  the 
subgrid  scale  channels  provides  a degree  of  freedom  for  return  flow  in 
the  presence  of  water  level  gradient,  which  would  otherwise  not  exist  in 
models  which  exclude  subgrid  scale  channels.  The  simulation  of  Hurricane 
Carla  in  particular  is  improved  over  that  attainable  with  the  SURGE  I 
program  which  did  not  allow  for  the  subgrid  scale  channel  subroutine. 

While  programs  such  as  SURGE  I can,  in  principle,  simulate  the  effects 
of  channels,  provided  the  grid  scale  is  of  the  order  of  the  channel  width, 
the  required  computer  time  is  usually  prohibitive  at  least  for  explicit 
numerical  models.  Some  advantage  can  be  gained  in  respect  to  economy  by 
the  use  of  implicit  numerical  models  such  as  that  of  Leendertse  (1967); 
however,  the  accuracy  of  such  schemes  when  used  on  a competitive  basis, 
from  the  standpoint  of  economy  (large  time  steps)  can  suffer  relative  to 
that  which  can  be  achieved  with  the  subgrid  scale  channel  routine.  How- 
ever, the  best  procedure  for  such  numerical  simulation  remains  to  be 
determined. 
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APPENDIX  A 


SURGE  II  PROGRAM 

This  appendix  includes  a comprlete  listing  of  the  SURGE  II  program. 
Except  for  SUBROUTINE  CHANL,  the  program  is  much  the  same  as  that  used 
in  Reid  and  Bodine  (1968).  I-t' should  be  emphasized  that  the  coding  of 
calculations  of  flow  and  water  level  for  blocks  does  not  include  the 
effect  of  Coriolis  force.  Moreover,  no  attempt  has  been  made  to  opti- 
mize the  coding  since  the  original  version.  The  actual  new  part  of  the 
program  is  embodied  in  SUBROUTINE  CHANL  and  the  way  in  which  the  channel 
computations  mesh  with  the  block  calculations.  Thus,  while  many  users 
may  prefer  their  own  version  for  calculations  over  the  main  grid,  it 
should  be  possible  to  incorporate  SUBROUTINE  CHANL  with  their  own  pro- 
gram when  applied  to  systems  like  the  Sabine-Calcasieu  region  in  which 
allowance  for  channels  is  essential. 
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APPENDIX  B 


DESCRIPTION  OF  THE  SURGE  II  CODED  PROGRAM 


The  general  strategy  of  the  program  is  discussed  and  certain  special 
features  are  pointed  out  which  may  not  be  apparent  without  detailed  study 
of  the  program.  Operational  aspects  of  the  program  are  discussed  in  some 
detail  in  Appendix  C. 


The  version  of  the  program  adapted  for  use  on  the  GE  400  computer 
system  by  the  Corps  of  Engineers  consists  of  the  following  parts  or 
subroutines : 


MAIN 

PART  2 

aiANL(l) 

aiANL(2) 

aiANL(3) 

CHANL(4) 

LIST(l) 

LIST(2) 

LIST(3) 

SAVE(l) 

SAVE(2) 

C0NTIN{1) 

rONTIN(2) 


whose  primary  job  is  to  read  and  check  the  sequencing  of  the 
basic  data  for  the  block  computations; 

which  controls  the  basic  computational  sequencing,  initializa- 
tion, and  updating  of  storage,  interpolation  of  coarse  wind 
fields  for  the  actual  grid,  and  routine  computation  of  U,  V, 
and  II  for  all  blocks,  considering  barriers  (basically,  the 
SURGE  I program) ; 

which  is  called  only  once  to  read  channel  data  and  to  estab- 
lish certain  key  arrays  for  routine  calculation; 

which  is  called  routinely  to  compute  flow  and  water  levels  in 
channels  and  at  channel  end  points; 

whose  task  is  the  routine  calculation  of  H on  blocks  con- 
taining channels; 

which  is  called  for  listing  of  channel  computations; 

which  is  called  only  once  to  read  control  data  for  block 
listings  and  to  list  the  topographic  Z field; 

which  lists  the  H field  for  blocks  if  called; 

which  lists  the  U,  V,  and  H fields  for  blocks  if  called  in 
place  of  LIST(2); 

which  is  called  only  once  to  read  the  positions  of  certain 
gage  locations  for  water  level  or  flow; 

which  is  called  routinely  at  preselected  time  intervals  to 
save  water  levels  and  flow  for  gage  locations  defined  by 
SAVE(l) ; 

which  is  called  only  once  to  read  basic  storage  in  COMMON 
BLOCKS  1 to  10  in  the  case  of  a continuation  of  a given 
problem; 

which  is  called  at  the  termination  of  a run  to  output  the 
continuation  data  called  for  by  CONTlN(l). 
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The  version  of  the  program  used  in  the  testing  and  calibration  work, 
using  an  IBM  3b0/65  computer  system,  has  an  additional  assembler  language 
subroutine  for  plotting  positions  of  barriers  and  channels  (see  Fig.  15), 
This  is  useful  in  checking  input  data  for  channels  and  barriers  to  spot 
possible  errors  in  coding  the  positions  of  channel  blocks  and  barrier 
blocks.  Unfortunately,  this  subroutine  is  not  compatible  with  the  GE  400 
system.  Subroutine  PLOT  in  Appendix  A however  can  be  used  for  this  pur- 
pose. Subroutine  LIST  is  not  used  in  the  version  of  the  program  in 
Appendix  A. 

1 . Flow  Diagram. 

A schematic  flow  diagram  for  the  SURGE  II  program  is  given  in  Figure 
B-1.  If  a new  problem  is  being  run  then  the  first  phase  is  reading  in 
the  basic  data  and  checking  the  data  sequencing  to  make  sure  it  is  in 
order  ;uid  complete.  This  is  carried  out  in  MAIN  and  the  beginning  of 
PART  2 which  calls  suliroutines  CUANL(l),  SAVE(l),  and  LIST(l). 

Initialization  of  block  arrays  is  carried  out  in  PART  2;  initializa- 
tion of  channel  arrays  and  establishing  of  key  arrays  are  carried  out  by 
aL\.NL(l).  These  key  arrays  are  discussed  in  a subsequent  subsection. 

Step  4 of  the  flow  diagram  is  the  beginning  (or  reentry  point)  of  the 
routine  computations  for  each  time.  After  generating,  the  detailed  inter- 
polated fields  of  X and  y components  of  wind  stress  for  the  blocks 
(step  4)  and  all  blocks  (i.e.,  all  I ,,J)  are  swept  to  compute  the  flow 
components,  U and  V,  ignoring  at  first  the  presence  (if  any)  of  sub- 
grid scale  channels,  but  considering  barriers  for  any  barrier  blocks 
(step  5) . 

In  step  6 aiANL(2)  is  called  to  sweep  through  all  channel  blocks  to 
evaluate  all  channels  Q and  11  except  those  for  11-end  points  and  all 
lateral  flows  to  and  from  channels.  In  the  latter  operation,  the  flows 
U and  V computed  in  step  5 are  replaced  by  corrected  U or  V be- 
tween blocks,  considering  the  presence  of  the  chajinels. 

Step  7,  which  is  carried  out  in  PART  2,  sweeps  all  I ,J  to  compute 
water  levels  on  blocks  ignoring  for  the  present,  the  presence  of  any 
subgrid  scale  channels. 

In  step  8,  CHANL(3)  is  called  to  correct  the  block  11  values  on 

those  blocks  containing  channels  and  to  compute  the  H and  Q values 

at  11-end  points  of  channels.  This  also  provides  corrected  H values 
for  those  blocks  into  which  the  channels  discharge. 

Steps  10  and  11  are  output  operations  for  block  and  channel  computa- 
tions carried  out  in  PART  2 and  aiANL(4).  This  is  followed  by  a time 
updating  and  test  for  end,  dependent  upon  a prescribed  maximum  number  of 

time  steps.  Before  termination  of  a run,  the  contents  of  all  data  in 

COMMON  arc  saved  for  possible  continuation  of  the  problem,  if  desired. 

2.  Identification  of  Adjacent  Channel  Blocks. 

To  provide  rapid  access  to  values  of  11  and  Q in  channels  adjoin- 
ing  a given  channel  reach,  special  arrays  are  generated  in  subroutine 
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Figure  B-1.  Generalized  flow  diagram  for  SURGE  II. 


□ lANL(l).  Tliere  are  four  sucli  arrays;  KCX(K) , KCY(K),  KCXP(K) , and 
KCYP(K) . These  give  the  channel  block  identification  index  for  those 
channel  blocks  which  are  adjacent  to  the  Kth  channel  block  as  indicated 
in  Figure  B-2.  llius,  KCX(K3  is  the  identification  of  the  channel  block 
which  has  an  x-side  channel  adjoining  channel  block  K on  the  negative 
characteristic  side  (i.e.,  on  a preceding  row),  while  KCXP(K)  is  the 
identification  of  the  channel  block  which  has  an  x-side  channel  adjoining 
channel  block  K on  the  positive  side  (i.e.,  on  a following  row).  KCY(K) 
and  KCYP(K)  have  analogous  meanings  for  blocks  with  y-side  channels 
adjoining  that  of  block  K.  These  arrays  are  generated  by  an  appropriate 
series  of  tests  in  which  the  I ,J  values  of  blocks  adjacent  to  that  of 
channel  block  K are  compared  with  the  ICG  and  JCG  values  of  all 
other  channel  blocks.  This  is  carried  out  only  once  during  any  run,  and 
is  not  particularly  time  consuming;  moreover,  it  avoids  any  human  error 
which  may  easily  occur  if  such  arrays  were  required  as  input. 


CHANNEL 

BLOCK 

KCXPk 

CHANNEL 

BLCJOC 

KCYk 

CHANNEL 

BLOa< 

K 

CHANNEL 

BLOCK 

KCYPg 

CHANNEL 

BLOCK 

KCX^ 

Figure  B-2.  Channel  block  identification  for 

channels  adjacent  to  those  of  block  K. 

The  arrays  KCX  and  KCXP  have  the  properties  KCXP(KCX(K))  = K and 
KCX(KCXP(K))  * K with  similar  relations  for  KCY  and  KCYP. 

As  an  exan|}le  of  the  use  of  such  arrays,  suppose  the  value  of  HC  in 
an  X channel  adjoining  that  of  channel  block  K is  needed.  This  could 
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be  addressed  as  HC(KX)  where  KX  = KXC(K) . Using  Figure  8 as  an 
example,  the  values  of  channel  flow  entering  the  junction  from  channels 
1 and  2 would  be  addressed  by  QCXP(Kl)  and  QCYP(Kl),  respectively, 
where  K1  designates  the  channel  block  containing  channels  1 and  2. 
However,  the  flow  leaving  the  junction  would  be  addressed  by  QCYN(K2) 
where  K2  = KaP(Kl)  and  QCXN(K3)  where  K3  = KCXP(Kl).  tVhile  redun- 
dant storage  of  such  H and  Q values  would  also  satisfy  the  require- 
ment of  rapid  access  to  such  values  adjoining  a given  channel  block,  the 
use  of  the  integral  arrays  KCX,  KCY,  KCXP,  and  KCYP  saves  storage  for 
most  com[)uter  systems. 

An  examination  of  the  listings  of  the  values  of  the  arrays  KCX,  KCY, 
KCXP,  and  KCYP,  as  output  by  the  program,  indicates  that  the  maximum 
value  of  any  of  these  can  and  usually  does  exceed  the  number  of  input 
channel  blocks  (KCM) . The  reason  for  this  is  that  dummy  storage  posi- 
tions are  created  for  blocks  adjoining  channel  end  points.  This  is  an 
artifice  of  the  program  which  allows  routine  computation  for  all  channel 
reaches  before  special  comjiutation  for  channel  end  points. 

3.  Barrier  Identification. 

The  position  of  the  Kth  barrier  block  is  given  by  the  array  pair, 
IB(K)  and  JB(K),  which  is  input  to  the  program.  It  is  convenient  to 
have  rapid  access  to  barrier  information  for  those  barriers  which  happen 
to  fall  on  a given  channel  block.  The  array  KCB(K)  gives  the  identifi- 
cation of  tl.f'  barrier  block  which  coincides  with  channel  block  K.  Thus, 
ICG(Kj  = IB(KCB(K))  and  JCG(K)  = JB(KCB(K)).  If  no  barriers  exist  in  a 
given  channel  block  then  the  corresponding  value  of  KCB  is  zero.  Thus, 
in  the  routine  program,  a test  for  zero  value  KCB  is  made;  if  nonzero, 
then  a call  can  be  made  for  barrier  data  such  as  elevation  and  barrier 
coefficients  via  the  barrier  index  KB  = KCB(KC)  where  KC  is  the 
channel  block  concerned. 

llie  array  KCB(K)  is  generated  in  CHAML(l),  via  a scan  of  all  IB 
and  ./B  values  for  given  ICG  and  JCG  for  channel  block  K. 

An  array  KLB(K)  is  also  generated  which  identifies  those  barrier 
blocks  not  common  to  channel  blocks.  This  is  ased  only  in  the  IBM  360/65 
assembler  language  plotting  routine,  not  in  routine  calculations. 

4.  Channel  F.nd-Point  Identification. 


As  a signal  that  at  least  one  diannel  end  point  occurs  in  a channel 
block  K,  the  value  of  ICG(K)  is  negative.  If  two  end  points  occur, 
the  value  of  JCG(K)  is  also  negative;  otherwise,  it  is  positive.  If 
no  channel  end  point  occurs,  then  both  ICG  and  .ICG  for  the  block  are 
positive.  This  positive-negative  coding  is  generated  automatically  in 
OlANL(l)  by  appropriate  testing;  namely,  to  check  if  a valid  channel 
connects  at  each  end  of  a valid  channel  in  the  block  concerned. 

In  addition,  the  arrays  KEN(1,K)  and  KF.N(2,K)  are  generated  in 
CIIAM.(l)  to  identify  the  type  of  end  point  for,  at  most,  two  potential 
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channel  terminations  in  given  channel  block  K,  If  there  is  no  channel 
termination  both  KKN(1,K)  and  K1:N(2,K)  are  zero;  if  one  termination 
occurs  for  block  K,  KEN(1,K)  will  have  an  integral  value  from  1 to  8 
and  KEN(2,K)  will  be  zero;  if  two  terminations  occur,  both  KI:N  arrays 
will  have  nonzero  value.  In  ase,  KI:N(2,K)  is  called  only  jf  JCG(K) 
is  negative. 

The  coding  for  the  type  of  end  point  is  indicated  schematically  in 
Figure  B-3.  Values  of  KEN  from  1 to  4 represent  "H-end"  type  termina- 
tions where  a ponding  block  immediately  adjoins  the  channel  end.  Values 
of  KEN  from  5 to  8 are  those  for  which  Q is  specified;  e.g.,  river 
discharge.  Values  within  either  group  indicate  the  relative  orientation 
of  the  channel  end  point  in  question  to  assure  calling  the  correct  data 
and  using  the  right  signs  in  the  routine  calculations. 


3.7 


4,8 


TYPE  OF  END 

KENk=  1,2. 3. 4 HC=H  OF  ADJACENT  BLOCK 
KENk=  5, 6, 7,8  Q SPECIFIED 

Figure  B-3.  Identification  of  type  and  orientation  of  a 

channel  end  point  by  the  coded  identifier  KEN(K) . 
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APPENDIX  C 


USER'S  GUIDE  TO  SURGE  II 

The  coded  progr.'un  SURGE  II  is  intended  for  use  in  the  numerical  simu- 
lation of  storm  surges  or  astronomical  tides  in  bays  and  estuaries  for 
specified  time  sequences  of  water  level  at  the  seaward  boundary  of  the 
bay  or  estuary  and  specified  wind  stress  and  other  storm  data  over  the 
bay  or  estuar>'.  The  user  may  use  one  of  two  distinct  modes  of  operation: 

(a)  the  storm  mode,  in  which  all  storm  data  are  required  as  well  as  sea- 
ward hydrograph  data;  or  (b)  the  tide  mode,  in  which  no  storm  data  are 
required,  the  only  forcing  being  the  input  water  level  variation  at  the 
seaward  boundary.  Moreover,  in  both  modes  the  user  has  the  option  of 
initiating  a new  simulation  or  continuing  a previous  simulation,  the 
input  requirements  being  different  for  each. 

In  general,  the  input  consists  of  the  following  types  of  information: 

(a)  Control  Data--For  input-output  operations,  initializa- 
tion, array  size,  time  stepping,  and  run  duration. 

(b)  Bay  Schematization  Data--including  block  topography, 
barrier  data,  and  channel  data. 

(c)  Forcing  Data--including  sequences  of  water  level  at 
seaward  boundary,  wind-stress  comjjonents  over  bay,  rainfall 
data  over  bay,  and  river  discharge  data. 

(d)  Problem  Specification  Information. 

Certain  checks  are  made  as  the  data  are  read  in,  with  regard  to  proper 
order  of  input,  proper  amount  of  sequential  data,  and  proper  size  arrays. 
All  stops  resulting  from  these  editing  checks  of  input  are  identified. 

In  the  subsequent  subsections,  the  individual  input  parameters  are 
identified  (with  appropriate  units),  the  sequence  of  data  input  for  the 
different  modes  of  operation  is  given  in  some  detail,  and  special  require- 
ments concerning  data  input  for  barriers  and  channels  arc  discussed, 
followed  by  a summary  of  output  information  and  output  options. 

1 . Definition  of  Input  Variables. 

The  following  variables  are  listed  in  the  order  in  which  they  are 
input  (asterisks  separate  data  blocks) ; 

ICARD  Control  index:  0 for  starting,  1 for  continuation. 

*****  Block  0 

I DENT  Data  block  identification; 
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IBl  starting  column  (I  value)  for  listing  of  block  H output 

(normal ly  taken  as  1) ; 

KCM  total  number  of  blocks  with  channels  (including  null  channels, 

see  subsec.  6 of  this  app.); 

NOK'IND  control  for  storm  data  input:  0 for  normal  input  operation  for 
wind  stress,  rainfall,  and  runoff;  -1  for  omitting  such  input 
for  tide  computations; 

INTER  interval  in  SAVE  operation  (time  interval  is  INTER*DELT) ; 

NGAGE  numlier  of  H gage  location.;  saved; 

NFLOW  number  of  Q gage  locations  saved; 

IMIN  minimum  expected  II  (feet) ; 

INIAX  maximum  expected  H (feet). 

NOTE IMIN  and  IMAX  are  used  only  in  subroutine  GRAF,  applicable  to 

IBM  360  or  370. 

*****  Block  1 

NTIME  Initial  time  level  (normally  0,  unless  a continuation  run  is 

being  carried  out,  in  which  case  NTIMi;  should  equal  the  final 
value  of  the  previous  run) ; 

NM  maximum  number  of  time  steps  for  the  problem; 

MMIN  minimum  "map  time"  for  wind-stress  input; 

MMAX  maximum  map  time  for  wind-stress  input; 

NFU  number  of  iterations  per  map  time  interval; 

lOirr  interval  for  routine  output  from  blocks  and  channels  equals 

I OUT  •*.  1; 

INFLI)  special  output  flag:  0 for  standard  output,  1 for  extra 

listing  of  channel  output  for  one  iteration  preceding  normal 
listing. 

*****  Block  2 

IM  Total  nunl)er  of  x-grid  intervals; 

.JM  total  nun*)er  of  y-grid  intervals; 

KM  total  nuralier  of  blocks  having  barriers; 
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KM\X 


total  nunil)er  of  coarse  x-griil  points  for  wind-stress  input; 
total  numlier  of  coarse  y-grid  points  for  wind-stress  input. 
Bloak  Z 


LMAX 


DfLX 

Spatial  grid  interval  or  block  size  (nautical  miles); 

Dl-LT 

time  interval  between  block  11  and  flow  computations 

(seconds) ; 

CDO 

overflow  coefficient  for  natural  low-lying  ground  such 
barrier  islands; 

as 

FK 

bed- resistance  coefficient  for  blocks; 

FC 

bed- resistance  coefficient  for  channels  (used  only  if 
for  individual  channels  are  not  entered); 

values 

IIGI 

initial  water  level  above  MSI.  in  the  bay  (feet). 

* ★ ★ ★ * 

Bloak  4 

KI 

Numljer  of  interpolation  subdivisions  of  each  coarse  x- 
interval  KI*(KMAX-1)  = IM; 

grid 

LI 

number  of  interpolation  subdivisions  of  each  coarse  y- 
interval  LJ*(LMAX-1)  = .JM; 

grid 

KIl 

number  of  coarse  x-grid  intervals; 

LJ.I 

number  of  coarse  y-grid  intervals; 

JBL,  JBR 

number  of  "open  boundary"  J-intervals  on  left  and  right 
of  system  (not  used  in  version  in  App.  A). 

Bloak  5 

IBdO 

I location  index  for  barrier  block  K; 

.IB(K) 

.1  location  index  for  barrier  block  K; 

IZX(K) 

elevation  of  x-barricr  (right  side)  on  barrier  block 
of  feet) ; 

K (tenths 

IZY(K) 

elevation  of  y-barricr  (upper  side)  on  barrier  block 
of  feet) ; 

K (tentlis 

I{:iK)X(K) 

overflow  coefficient  for  x-barrier  (value  x 1,000)  on 
barrier  block; 

Kth 

164 


I i 

i 

■ 

« 


ICDOY(K) 

ICDSX(K) 

ICDSY(K) 


IMRO 

•JMRO 

KR 

ISTR 

IND 

NOW 

KIM 


NORT 

***** 

RI 

CONST 

S 

***** 

LROI(K) 

TRO.I(K) 


overflow  coefficient  for  y-barricr  (value  x 1,000)  on  Kth 
barrier  block; 

subraerged  wier  coefficient  for  x-barrier  (value  x 1,000)  on 
Kth  barrier  block; 

subraerged  wier  coefficient  for  y-barrier  (value  x 1,000)  on 
Kth  barrier  block. 

Block  6 

elevation  of  ground  or  seabed  (feet)  relative  to  MSL  datum 
for  block  location  I,.J. 

Block  7 

Number  of  river  input  (runoff)  locations; 

number  of  map  times  with  runoff  values; 

number  of  channel-stress  values  (normally  same  as  JMRO); 

start  of  rain  (map  time); 

end  of  rain  (map  time); 

number  of  iterations  between  river  input  values  (normally 
same  as  NFU) ; 

number  of  iterations  between  channel-stress  values  (normally 
same  as  NFU) ; 

number  of  iterations  per  hour  for  rain  (normally  same  as  INTHR) . 
Block  8 

Total  rainfall  (inches); 

fraction  of  rainfall  not  absorl>ed  by  ground; 

conversion  factor  for  wind  stress  (5 ,280/3,()()()) ^ x l,l/io. 

Block  9 

I location  index  for  Kth  river  input  block; 

J location  index  for  Kth  river  input  block. 
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***** 


Block  10 


DIST(M)  Percent  of  total  rainfall  per  hour  for  24  hours. 

*****  Block  11 

CHST(M)  Qiannel-stress  values  at  map  time  M (entries  are  used  only 
if  KCM  = 0) . 

*****  Bloak  12 

R0(K,M)  Discharge  (cubic  feet  per  second)  from  Kth  river  input  block 
at  map  time  M. 

*****  Bloak  12 

MTIMT;  Map  time  for  given  block  of  wind-stress  input  and  seaward 

water  level. 


*****  Bloak  14 

HGR(K)  Seaward  water  level  above  MSL  (feet)  at  Ml'IME  for  coarse 
grid  position  K. 

*****  Bloak  IS 


HBR(J)  Water  level  on  right  open  boundary  above  MSL  (feet)  at 

NfTIME  for  grid  position  J (not  used  in  version  in  App.  A). 

*****  Bloak  16 


XR(K,L)  Wind-stress  component  in  the  x direction  (units  of  (miles 

per  hour)^/10)  for  coarse  grid  position  K,L  at  time  MTUE. 

*****  Bloak  17 

YR(K,L)  Wind-stress  component  in  the  y direction  (units  of  (miles 

per  hour) 10)  for  coarse  grid  position  K,L  at  time  MTIME. 

*****  Block  18 


ICG(K)  I location  index  for  channel  block  K; 

JCG(K)  J location  index  for  channel  block  K; 

IWCX(K)  width  of  x channel  (right  side)  on  channel  block  K (feet), 

with  sign  (see  subsec.  6 of  this  app.); 

IZCX(K)  depth  of  x channel  bed  on  channel  block  K (feet),  with  sign 
(see  subsec.  6 of  this  app.); 
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iiva(K) 


(feet) , 


width  of  y channel  (upper  side)  on  channel  block  K 
with  sign  (see  subsec.  6 of  this  app.); 

IZCY(K)  depth  of  y channel  bed  on  channel  block  K (feet),  with  sign 
(see  subsec,  6 of  this  app.); 

IFC(K)  bed-resistance  coefficient  for  channels  on  block  K (value  x 
10,000),  if  entry  is  zero  (blank)  then  IFC  is  taken  as  FC 
(entered  in  Block  J)  x 10,000. 

*****  Blodk  19 

I(1AGE(K)  Location  index  for  the  Kth  hydrograph,  if  JGAGE(K)  ^ 0 then 
IGAGE(K)  is  the  I location  of  a block  H;  if  JGAGE(K)  = 0 
then  IGAGE(K)  is  the  channel  block  index  for  a channel  H; 

JGAGE(K)  if  not  zero,  this  is  the  J location  of  a block  H;  if  zero, 
a channel  II  is  indicated; 


KFLOW(K)  channel  block  index  for  the  Kth  flow  gage,  the  flow  being 
that  of  the  lower  end  of  the  x channel,  or  the  left  end  of  a 
y channel  if  an  x channel  does  not  exist,  or  a channel  end 
point  if  one  exists  in  the  identified  ch;mnel  block., 

*****  Block  20 


lEND 

NF 

I BEGIN 


N.1 


NCARD 

ALPIIA(J) 


Maximum  I in  listing  of  block  arrays  of  II,  U,  and  V; 
numljer  of  iterations  between  listings; 
first  I in  listing  of  block  arrays; 
maximum  .1  in  listing  of  block  arrays; 

total  number  of  alphanumeric  problem  identification  cards; 

alphanumeric  character  data  which  identify  the  problem  and 
gage  locations  by  name. 


2.  Input  for  Initiating  Storm  Surge  Simulation. 

The  sequence  of  input  for  starting  a problem  in  the  storm  surge  mode 
is  given  below  in  the  form  of  a summary  of  the  REAP  statements  active  in 
this  mode,  together  with  a summary  of  the  appropriate  FORMATS  for  data 
input  in  different  blocks.  For  all  data  blocks  requiring  an  entry  of 
the  identification  integer  IDENT,  only  the  units  digit  of  the  data  block 
number  is  entered  in  column  1 of  the  data  input  card. 


Control  Card 


RI;AD  1 , I CARD  (0  for  starting) 

Block  0 (1  card) 

RIIAU  1 , IDENT,  IBL,  KCM,  NOWIND,  INTER,  NGAGE,  NFLOW,  IMIN,  IMAX 
NOTE IMIN  and  IMAX  are  left  blank  unless  subroutine  GRAF  is  used. 

Block  1 (1  card) 

READ  100  , IDEiVr,  NTIME,  NM,  MMIN,  MMAX,  NFU,  lOUT,  INFLD 
Block  2 (1  card) 

READ  100  , IDENT,  IM,  JM,  KM,  KMAX,  LMAX 
Block  2 (1  card) 

READ  250  , IDENT.  DELX,  DELT,  CDO,  FK,  FC,  HGI 
Block  4 (1  card) 

READ  100  , IDENT,  KI , LJ , KII,  LI J , JBL,  JBR 
Block  5 (total  of  KM  cards  of  barrier  data) 


DO 

500 

K = 1,  KM 

READ 

100  , 

IDENT,  IB(J),  JB(K).  IZX(K) , IZY(K) . ICDOX(K),  ICDOY(K) , 
ICDSX(K),  ICDSY(K) 

500 

CONTINUE 

Block 

6 (total  of  2*IM  cards  of  block  topography) 

DO 

550 

I = 1,  IM 

READ 

100  , 

IDENT,  (IZ(I,J),  J = 1,10) 

READ 

100  , 

IDENT,  (IZ(I,J) , J = 11,  JM) 

550 

CONTINUE 
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Bloak  7 (1  card) 

RJiAD  100  , ini;.\T,  IMRl),  JMRO,  KR.  ISTR,  IND,  NOW,  KIM,  NORT 
Bloak  fl  (1  card) 

RliAD  250  , IDHNT,  RF,  CONST,  S 

Bloak  9 (1  or  2 cards,  dependent  on  IMRO) 


RliAD 

100  , I DENT, 

(LROI(K) , 

LRaJ(K) , K = 1, 

IF  (IMRO.LT.b)  GO  TO 

575 

READ 

100  , I DENT, 

(LROI(K)  . 

LRai(K) , K = 6, 

575  CONTINUE 

Block 

10  (3  cards) 

ri;ad 

250  , I DENT, 

(DIST(M), 

M = 1,10) 

READ 

250  , IDENT, 

(DIST(M) , 

M = 11,20) 

READ 

250  , IDENT, 

(DIST(M) , 

M = 21,24) 

Bloak 

22  (L  + 1 card 

1 where  L • 

= KR/10.  If  KR 

omitted. ) 

READ 

250  , IDENT, 

(CIIST(K) , 

K = 1,11) 

1 RliAD 

1 

250  , IDENT, 

(aiST(K) , 

K = 11,20) 

1 READ 

250  , IDENT, 

(aiST(K) , 

K = KL,  KR  (KL 

1 

1 Bloak 

22  (JMRO  cards 

of  river 

discharge  data) 

i DO 

700  M = 1, 

JMRO 

\ RliAD 

250  , IDENT, 

(R0(K,M), 

K = 1,  IMRO) 

700  CONTI  NUF. 

Vind  Streas  and  Water  Level  Forcing 

(MIL  sets  of  blocks  13  to  17  where  MIL  = MMAX  - MMIN  + 1) 

710  CONTINUE 
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Bloak 

12  (1 

card) 

RliAD 

100  , 

IDENT,  m’lME 

Bloak 

14  (1 

card) 

READ 

250  , 

IDENT,  (IIGR(K),  K = 1,  KMAX) 

Bloak 

IS  (1 

card) 

R1;AD 

250  , 

IDENT,  (IIBR(.I),  .1  = 2,8) 

Bloak 

16  (KMAX  cards) 

DO 

790 

K = 1,  KMAX 

READ 

250  , 

IDENT,  (XR(K,L),  L = 1,  LMy\X) 

790 

CONTINUE 

Bloak 

17  (KMAX  cards) 

DO 

800 

K = 1,  KMAX 

ri-:ad 

250  , 

IDENT,  (YR(K,L),  L = 1,  I.MAX) 

800  CONTINUE 

IF  (WIMli  - MMAX)  710,  1,015,  1,015  (710  returns  to  read  block  13) 

1,015  (CONTINUE) 

Bloak  18  (KCM  cards  with  channel  data.  If  KCM  = 0,  the  READ  statement 
is  bypassed  and  block  18  should  be  omitted.) 

IF  (KCM. or. 0)  CALI,  CUANL(l) 

DO  50  K = 1,  KCM 

R1;AD  100  , IDENT,  ICG(K)  , .ICC(K),  IWCX(K),  IZCX(K),  IWCY(K),  IZCY(K), 

IFC(K) 

' I 

50  CONTINUE  i 

I ) 
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Block  19  (2  cards) 
CALL  SAVL(l) 


RFAD 

350  , 

(IGAGF(K),  JGAGn(K).  K = 

1,  NGAGE) 

RLAD 

350  , 

(KFLOW(K),  K = 1,  NFLOW) 

Block 

20  (NCARD  + 1 card) 

CALL 

LIST(l) 

RFAD 

1 , 

IDFNT,  IFND,  NF,  I BEGIN, 

NJ,  NCARD 

DO 

250 

J = 1,  NCARD 

RFAD 

450  , 

(ALPMA(J),  J = 1,40) 

250 

CONTINUE 

Format 

Statements  for  Input.  The  following  formats  were  used  in 

all  the  testing  operations.  It  is  recommended,  however,  that  for  routine 
operations  those  READ  statements  using  FORMAT  1 be  replaced  by  FORMAT  100 
to  make  all  basic  numerical  input  consistent  in  card  column  range. 


1 FORMAT  (II,  13,  19,  14) 

100  FORMAT  (II,  2X,  15,  9(3X,  IS) 

250  FORMAT  (II,  F7.0,  9F  8.0) 

350  FORMAT  (20  I 4) 

450  FORMAT  (15A2,  15A2,  10A2) 

3.  Input  for  Tide  Mode. 

For  calibration  of  a given  bay  system,  under  virtually  no  wind  con- 
ditions, for  its  response  to  forcing  by  astronomical  tide  at  the  seaward 
boundary  and  a steady-state  river  discharge,  allowance  is  made  in  the 
coded  program  to  bypass  the  detailed  input  of  wind-stress  con^ionents, 
and  rainfall  and  channel-stress  data;  moreover,  since  a steady  river 
discharge  is  assumed  only  a single  card  is  required  to  define  this  input. 
In  essence,  the  data  blocks  10  to  17  are  replaced  by  a shortened  version 
of  block  12  plus  a modified  version  of  block  14  in  which  tide  data  at  the 
seaward  boundary  are  prescribed  at  hourly  intervals  as  the  map  time  inter- 
vals. The  input  is  summarized  as  follows: 
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Control  Card:  0 in  column  1 


Bloak  0:  see  Section  IV, 1,  NOWIND  = -1 

Blocks  1 to  y see  Section  IV, 2 

Bloak  12  (1  card  for  steady  river  discharges) 

RI-;AD  250  , IDHNT,  (R0(K,M),  K = 1,  IMRO) 

Astrotide  Bloak  (1  card  for  each  12  hours) 

905  RliAD  910,  IGA,  MTIME,  (H(1,J),  J » 1,12) 

MU  = NTTIMli  + 12 
IF  (MU.LT.MMAX)  GO  TO  905 
910  FORMAT  (12,  14,  12F  6.2) 

(IGA  = 1) 

Blocks  Id  to  20:  see  Section  IV, 2 

Comments  on  Tide  Mode.  'Ilie  map  time  interval  for  the  tide  mode  is 
1 hour.  The  MTIME  entry  for  the  astrotide  block  is  the  time  (hour)  of 
the  first  of  12-hourly  values  of  HG  (entered  as  U(1,.I)).  The  tide  is 
assumed  uniform  along  the  seaward  boundary  of  the  bay  system,  hence  one 
HG  value  per  hour  is  sufficient. 

In  starting  the  tide  mode  from  rest  state  (U  = V = 0 and  II  = IIGI  = 0) , 
usually  one  or  two  diurnal  tide  cycles  are  required  for  the  numerical 
model  to  reach  a nearly  periodic  response  to  an  almost  periodic  input. 
Thus,  if  the  final  diurnal  cycle  is  to  be  free  of  initial  transients,  at 
least  72  hours  of  HG  data  should  be  provided,  lliis  may  require  an 
adjustment  in  the  dimensions  given  in  C0MM0N/BLK6/  which  appears  in  sub- 
programs M/MN,  PART  2,  and  CONTIN,  if  the  full  data  set  is  to  be  stored 
for  one  run.  An  alternative  is  to  make  use  of  the  continuation  option, 
using  less  data  input  per  run  (e.g.,  24  hours). 

4.  Input  for  Continuation  of  a Run. 

Since  the  main  purpose  of  the  tide  mode  is  for  calibration  of  the  bed 
friction  coefficients  for  blocks  and  channels,  it  is  expected  that  many 
trial  runs  will  be  made  for  a given  bay  system.  In  order  to  keep  the 
machine  time  to  a minimum  for  each  successive  run,  it  is  desirable  to 
ase  an  initial  field  of  IJ,  V,  or  H which  is  close  to  the  true  re- 
sponse at  the  starting  time.  This  can  be  accomplished  by  using  the  re- 
sulting U,  V,  and  H arrays  from  a previous  tide  run  for  the  bay  system 
as  the  initial  values.  (lliis  should  be  done  even  if  the  previous  run  has 
different  values  of  the  bed  friction  coefficients.)  The  mechanism  for 
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accomplishing  this  is  the  use  of  the  continuation  mode  option,  as  con- 
trolled by  ICARD.  In  this  mode,  the  contents  of  common  from  a previous 
run  arc  input  along  with  any  additional  forcing  function  data. 

To  make  tlie  program  as  flexible  as  possible,  the  continuation  option 
can  be  used  for  eitlier  storm  surge  problems  or  astrotide  problems,  the 
only  difference  in  input  being  in  the  type  of  forcing  function  input. 
Such  forcing  function  data  should  be  consistent  with  the  continuation 
time.  Moreover,  the  value  of  NTIMI;  input  in  data  block  1 should  be 
equal  to  the  final  NTIMh  in  the  previous  run  which  is  continued. 


The  sequence  of  input  for  continuation  of  a problem  is  as  follows: 
Contvol  Card:  1 in  column  1 

Contents  of  COMMON  output  from  a previous  run 
see  subsection  2 of  this  app.  (4  cards) 


Contin  Deck: 
Blocks  0 to  3: 
Forcing  Deck: 


For  storm  surge  mode,  blocks  13  to  17,  inclusive. 
For  tide  mode-astrotide  deck. 


A flow  diagram  summarizing  the  RFAD  operations  as  controlled  by  ICARD 
and  NOWIND  is  given  in  Figure  C-1. 

5 . Comments  on  Barrier  Input. 

a.  Possible  Barrier  Locations.  All  barriers  in  the  schemati zation 
occur  parallel  to  the  sides  of  a given  barrier  block.  Barrier  data 
qualified  by  an  X in  the  coded  name  (e.g.,  IZX,  ICDOX,  ICDXS)  refer 

to  barriers  normal  to  the  x-axis  on  the  right  side  of  the  barrier  block; 
those  qualified  by  a Y in  the  coded  name  (e.g.,  IXY,  etc.)  refer  to 
barriers  normal  to  the  y-axis  on  the  upper  side  of  the  barrier  block. 

If  a channel  exists  parallel  to  either  barrier,  then  such  a barrier  may 
occur  on  either  or  both  sides  of  the  parallel  channel,  depending  upon 
the  coding  of  the  associated  channel  input  data  (as  discussed  in  a sub- 
sequent subsection).  Barriers  which  might  exist  along  the  left  or  lower 
side  of  a given  block  arc  represented  by  appropriate  data  coding  of  a 
barrier  block  in  a previous  row  or  column. 

b.  Precaution.  It  should  be  emphasized  that  for  any  barrier  block 
it  is  up  to  the  user  to  supply  appropriate  barrier  elevations  ZB  for 
both  the  right  and  upper  sides  of  the  barrier  block  even  if  a real 
barrier  occurs  only  on  one  side  of  the  block.  The  important  point  to 
observe  is  that  the  specified  ZB  values  should  aUjcqjS  equal  or  exceed 
the  larger  of  the  block  elevations  at  or  adjacent  to  the  side  of  the 
barrier  block  in  question.  Otherwise,  errors  can  occur  in  the  corajiuta- 
tions . 


c.  Array  Size.  The  number  of  barrier  blocks  KM  is  normally 
limited  to  less  than  100. 
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Figure  C-1.  Flow  diagram  for  read  statements. 
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6.  Comments  on  Channel  Input. 

a.  Possible  Channel  Locations.  All  channels  in  the  schematization 

occur  along  the  right  side  or  the  upper  side  of  a given  channel  bloc' . 
Channel  data  qualified  by  an  X in  the  coded  name  (e.g.,  IWCX,  IZCX> 
refer  to  channels  normal  to  the  X-axis  on  the  right  side  of  the  channel 

block;  those  qualified  by  a Y in  the  coded  name  (e.g.,  IWCY,  IZCY) 

refer  to  channels  normal  to  the  Y-axis  on  the  upper  side  of  the  channel 

block.  If  a block  has  both  an  X and  Y channel,  one  data  card  speci- 

fies both. 

b.  Channel  Junctions.  In  the  schematization  of  a channel  system 
junctions  can  occur  with  adjoining  channel  reaches  parallel  to  each 
other  or  perfiendicular.  Moreover,  one-,  two-,  or  three-way  branches  are 
possible. 

Four  possible  right-angle  channel  junctions  are  illustrated  in  Figure 
C-2.  The  simplest  junction  is  that  shown  in  the  upper  right  panel  of  the 
figure  where  the  joining  channel  reaches  are  in  the  same  channel  block 
Kl.  Right-angle  junctions  involving  two  adjacent  channel  blocks  are 
illustrated  in  the  upper  left  and  lower  right  panels  of  Figure  C-2. 


CHANNEL 

BLOCK 

K3 

^ 

CHANNEL  1 
BLOCK  1 

Kl  I 

1 

NULL  CHRNNEL  i 
BLOCK 

(JUcx-iycY-oy 

c . 

CHANNEL 

BLOCK 

K2 

1 

CHANNEL 

BLOCK 

Kl 

CHANNEL 

BLOCK 

Kl 


Figure  C-2.  Four  possible  simple  bends  for  a channel  reach. 
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The  final  possible  right  turn  is  illustrated  in  the  lower  left  panel 
of  the  figure.  In  this  case,  the  program  requires  that  a channel  block 
(Kl)  join  the  connecting  channels  of  the  nonadjoining  channel  blocks 
(K2  and  K3)  even  though  no  channels  exist  on  the  joining  block  Kl.  In 
such  circumstances,  the  required  "null"  channel  block  would  have  zero 
width  for  both  the  X and  Y channels  (IWCS  = IWCY  = 0)  <as  input.  The 
H value  at  the  junction  of  the  connecting  channel  reaches  for  this  case 
is  stored  as  HC(Kl) ; i.e.,  in  association  with  the  null  channel  block. 

Colinear  adjoining  chiuinels  always  involve  two  adjacent  channel 
blocks.  Four  possible  junctions  of  this  type  are  illustrated  in  Figure 
B-2  in  relation  to  central  channel  block  K. 


c.  Channels  with  Levees.  The  program  allows  for  the  following  possi- 
ble situations  with  respect  to  barriers  parallel  to  channels: 

(a)  Single  barrier  on  the  "inner"  lateral  boundary  of  a 
cliannel ; 


(b)  single  barrier  on  the  "outei"  lateral  boundary  of  a 
channel ; 

(c)  barriers  of  equal  elevation  on  both  sides  of  a channel. 


.NOTE. --The  term  inner  or  outer  side  of  a channel  refers  respec- 
tively to  the  side  common  to  the  channel  block  contain- 
ing the  channel  or  the  side  common  to  an  adjacent  block. 


Tlie  barrier  elevation  information  is  input  separately  from  the 
channel  block  data  and  allows  only  one  elevation  for  the  rigfit  side  and 
one  for  the  top  side  of  a block  (hence,  the  restriction  of  equal  barrier 
lieights  for  the  double  levee  situation  c above).  The  specification  for 
situations  a,  b,  or  c is  accomplished  by  a sign  coding  in  the  channel 
block  data  as  follows; 


(a)  Channel  width  (ICW)  positive,  channel-bed  elevation  (IZC) 

(b)  channel  width  positive,  channel-bed  elevation  positive; 

(c)  channel  width  negative,  channel-bed  elevation  negative. 

It  is  understood  that  only  the  magnitude  of  IWC.  and  IZC  for  a given 
channel  is  used  in  calculations. 

d.  Channel  Terminations.  A channel  system  can  terminate  at  (a)  a 
larger  body  of  water  representing  a lake,  bay,  or  sea;  or  (b)  at  a bound- 
ary or  in  a landlocked  block  within  the  system.  In  the  second  case,  the 
program  assumes  that  the  flow  at  the  channel  end  is  zero  unless  a river 
discharge  to  the  channel  is  specified  (see  input)  and  that  the  channel  end 
block  is  one  block  inside  the  boundary  block. 
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e.  Restriction.  Only  channels  with  the  channel  bed  below  the  mean 
water  level  (MWL)  reference  are  allowed.  Tlie  actual  elevation  used  in 
calculations  is  - |lZC|,  regardless  of  the  sign  on  the  input  of  IZC  for 
a given  channel. 

f.  Array  Size.  The  numlier  of  channel  blocks  (including  null  channel 
blocks)  is  KCM.  However,  (aiANL(l))  creates  arrays  of  length  KCMI'  > KCM. 
The  value  of  KCMF’  exceeds  KCM  by  one  plus  the  number  of  channels  which 
terminate  on  the  exterior  boundary  of  the  grid  including  the  seaward  bound- 
ary. Since  KCMP  is  limited  to  130,  KCM  should  bo  less  by  the  amount 
described  above. 

7.  Output . 

a.  Listings  of  Input  and  Key  Arrays.  Ail  input  data  are  listed  in 
easily  identifiable  form  in  the  order  in  which  the  data  are  entered 
through  block  18.  Immediately  following  the  basic  channel  input  is  a 
listing  of  the  key  arrays  for  channels,  as  discussed  in  Appendix  A, 
including  the  assignment  of  sign  coding  for  ICG  and  JCG. 

Also  printed  out,  in  the  same  block  format  as  the  routine  listings  of 
II,  are  the  block  elevations. 

b.  Sequential  Output.  Normally,  the  routine  output  of  computed  values 
includes  block  H arr.ays  and  listings  of  all  channel  variables  at  pre- 
determined intervals  oi  time  (as  determined  by  lOUT) . It  is  possible  to 
list  the  U,  V,  and  H arrays  for  blocks  by  changing  the  CALL  LIST(2) 
statement  following  statement  2,100  in  PART  2 to  CALL  LIST(3). 

For  diannel  listings,  refer  to  Figure  6 for  notation;  the  listings 
are  ordered  by  channel  block  number  K.  The  block  location  I , .1  is 
repeated  (negative  signs  indicating  end  points).  This  is  followed  by 
IIX,  the  water  level  (feet)  and  QXN,  the  volume  transport  (cubic  feet 
per  second)  at  the  lower  end  of  the  x channel,  then  QCP , the  transport 
at  the  upper  end  of  the  x channel.  These  are  followed  by  IlY,  QYN,  and 
QYP  representing,  respectively,  the  water  level  and  flow  at  the  left  end 
and  flow  at  the  right  end  of  the  y channel.  Next  is  HC,  the  water  level 
at  the  junction  of  the  x and  y channels.  The  last  four  entries  in  the 
channel  listings  are  the  transports  (in  cubic  feet  per  second)  to  the 
channel  from  the  channel  block  and  from  the  channel  to  an  adjacent  block 
fo.'  the  X and  y charmels.  Tlie  IIC  value  is  meaningful  for  null  channels 
only. 

c.  Saved  Time  Sequences.  Subroutine  SAVE,  if  used,  saves  sequences 
of  water  level  and  flow  at  preselected  locations  (as  identified  in  block 
19  of  the  input).  In  the  original  version  of  the  subroutine  used  with  an 
IBM  360-65  computer  the  saved  information  was  punched  on  cards  to  facili- 
tate later  graphing  of  the  sequences. 
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APPENDIX  D 


COMPLETE  DATA  LISTING  OF  INPUT  FOR 
SABINE-CALCASIEU  REGION  WITH 
FORCING  DATA  FOR  HURRICANE  CARLA 
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-.00®® 

- . 0076 

-.00*5 

- .0056 

6 

-.020 

-.0225 

-.oqo® 

- .0076 

-.0085 

- . OOS6 

6 

-.02b 

- . 0225 

.."0®® 

-.0076 

-,"085 

-.0056 

7 

-.011 

-.  0037 

-.0055 

-.002® 

-.0025 

• . 002o 

7 

-.oow 

- , 0055 

-.0050 

-.oo?® 

-.0020 

-.0018 

7 

-.oo<* 

-.0057 

-.0050 

- . 0 0 2 .1 

-.0017 

-.0015 

7 

-.007 

- , Ou  36 

-.0051 

-.00?" 

-,0017 

-.0015 

7 

-.00e> 

-.0030 

-.002® 

-.0017 

-.0015 

-.0013 

7 

-.005 

-.oo«« 

-.0018 

- . 0 V 1 6 

-.0015 

-.0011 

7 

-.005 

-.00«<» 

-.0018 

-,"016 

-,0015 

-.0011 

7 

-.005 

-.O0“« 

-.0018 

-.0016 

-.0015 

-.0011 

MG®  FOLLH.f.O  BY 

Mb®  tHB*Y  (JM  FEFT)  *T 

M7 I“t« 

1 

5.500 

5. 5000 

5,5000 

5.5,100 

5,5000 

5 , S 0 ii  0 

5,5o00 

5 

<1.500 

U.5000 

6.5000 

6.5000 

6, 5c 00 

6,5000 

6.50"0 

X9  VlLHtS{  ID''T«6)  t\ri  YH  YIUI'FS 

0101x7*7)  *7  “iPTlYf® 

6 

-.O'i 

-.0112 

-.0112 

- . 0060 

-.0070 

- • 0056 

-.051 

-.0200 

-.0102 

- . 0 0 0 1 

-.0071 

- .0056 

6 

-.052 

-.0101 

-.0156 

-.Pu®l 

-.007] 

-.0055 

6 

-.052 

-.0172 

-.0162 

-.00®2 

-.0072 

-.0055 

6 

-.o«2 

- . 0 1 uu 

-.0117 

-.<1034 

-.0073 

-.0055 

6 

-.050 

- . 02Ui. 

-.0118 

-.0o*5 

-.0073 

-.0086 

6 

-.050 

• . n^Ub 

-.0118 

-.00*3 

- , 0 7 3 

- . 0066 

6 

-.050 

- . n2«b 

-.0118 

- . Ou«  J 

-.0073 

-.0086 

7 

-.015 

-.00x5 

-.0065 

- . 0,155 

-."027 

-.0021 

7 

-.011 

-.008® 

-.0355 

- . Oci5  1 

-.0025 

-.101® 

7 

-.010 

-.OyJi 

-.0058 

-.0050 

-.0025 

- , 0 0 1 8 

7 

-.00® 

-.0066 

-.0061 

-.0027 

-,002p 

-.0017 

7 

-.008 

-.  0035 

-."02® 

- . 0 ,1 2 3 

-.0020 

-.0015 

7 

-.008 

-.0068 

- . 0025 

-.001 » 

-.0017 

-.0015 

7 

- , OOo 

-.006* 

-.0025 

- . 0 0 1 8 

-.0017 

• . Oo  1 5 

7 

-.006 

-.0068 

-."025 

-.001  * 

-."017 

-.0015 

BtSl-AVAllABlt  COPV 
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^OLLOk-FO  by  «««  *nR*y  ( I v.  r^FTi  4T  "ir't»  3 


u 

b.ano 

6,?000 

5, 2000 

5 , 2 1)  0 0 

5.20O0 

8,2000 

5.2000 

5 

5.200 

^.?000 

5.2000 

5.2OO0 

5,2(100 

5. 2000 

5,2b00 

XB  y*uueS(  I0M«6)  *no 

yH  V*LU£S 

f lOMT,?) 

»T  ►'aPTI-BB 

b 

-.03T 

-,0l3b 

-.0123 

-."1 1 1 

-.00*8 

-.0053 

6 

• • 03b 

- . 02 1 b 

-.0112 

•.01"1 

-,0090 

- , 0080’ 

6 

• # 0 3b 

-.0159 

-.0130 

-.0102 

-.OOO] 

-.0081 

6 

• * 0 3b 

-.0187 

-.0155 

-.0103 

-.0091 

-.0081 

6 

-.03« 

-.0103 

-.0102 

- . 0 1 0(j 

-.0002 

-.00*2 

6 

-.055 

-.0202 

-,010B 

• a 0 09u 

-.00*3 

-."073 

6 

-.035 

- . 02*2 

-.0100 

- . 009u 

-.00*3 

-.0073 

b 

-.035 

-.0282 

-.0100 

• . OOOii 

-.00*5 

-.0073 

7 

-.015 

- , 0 0 5b 

-.0  052 

-.0007 

-.003* 

-.0022 

7 

-.013 

-,  0079 

-.0003 

-.00X0 

-."oxs 

-.0031 

7 

-.012 

•,0055 

-.0005 

-.  0035 

-.0031 

-,0028 

7 

-.011 

- , 0 0 5 0 

-.0008 

- . 0 0 3 2 

-.0050 

• , 0 028 

7 

-.009 

- ,003b 

-.0038 

-.0020 

-.0027 

• .oo2b 

7 

-.007 

-.0050 

-.  00  33 

-.  0022 

-."021 

-.0018 

7 

-.007 

- , 0050 

-.0033 

-.0022 

-,O02l 

-.0018 

7 

-.007 

- a OObO 

-.0033 

-.0022 

-.0021 

-.0016 

HRB  follo-fo  py 

MgP  »WB4y  (Jli  FfrT)  4T 

3 

a 

b . 0 00 

5,0000 

5. 0000 

b.OOOO 

0 . op  On 

b , no  00 

b.OOOO 

5 

5.200 

5.2000 

5.2000 

5 . 2 0 0 0 

5.2pnc 

5.2000 

5.2u00 

XB  v*LUE5(  ID''‘TBb)  *5C 

VIL'IF* 

f 1011797  ) 

47  “APOIffE* 

b 

- .000 

-.0138 

-.0125 

-.0112 

. , n 1 op 

-.0089 

b 

-.038 

-.0201 

-.0125 

-.0102 

-.0091 

-.00*1 

b 

-.039 

-,"20  3 

-."130 

- . 0 1 0 1 

-."092 

-.00*1 

b 

-.037 

• . 0205 

-.0102 

• . 0 1 oa 

-.0092 

-.0082 

b 

-.035 

-.0170 

-.0158 

-.0105 

-.00*3 

- ,0o73 

b 

-.033 

-.0305 

-.0105 

-.OCBj 

-.00*3 

-,0u7- 

b 

-.033 

-.0305 

-.0105 

- . "OOo 

- , 0 fi » j 

- , no7<. 

b 

-.053 

• , 0 305 

- , 0 4 5 

- . OiiO(4 

-.00*3 

• , no  7« 

7 

-.015 

-.0053 

-.0008 

- . 0 0 0 J 

-.OOBO 

-,0038 

7 

-.013 

-,"0bb 

-.0000 

-.0035 

-.0031 

-.0028 

7 

-.011 

-.0059 

-.0  052 

-.0032 

- . 0 0 2 • 

-.0025 

7 

-.009 

-.0051 

-.0038 

• . 0 u 28 

-.0027 

-.002b 

7 

-.007 

-.0037 

-.0037 

- , 0 u 2(1 

-.0021 

-.0018 

7 

-.005 

-.0050 

-.0025 

-.001 0 

-.0018 

-.0015 

7 

-.005 

-.0050 

- . 0 0 2 b 

-.0018 

- , on  1 * 

-.001b 

7 

- . 005 

-.0050 

• . 002b 

-.0018 

-.0018 

-.OOlb 

6.?(300 

z 


6.00^0 

3 


HCS  FOLLOwri)  BV  HHH  *B»AV  (IS  FEfTI.  IT  nTI^ts  u 


II 

5.eno 

5,0000 

5.aoou 

5.8uOn 

5.«000 

5,6u00 

5.8000 

5 

5.200 

5.2000 

5.2000 

5 . 2 0 0 0 

5,2ooo 

5,2000 

5.2000 

XP  V4LUt8( 10MT»6)  ASO  xr 

values 

( lOMT* 

7)  aT  -aPTIbE* 

6 

• .052 

-.0205 

-.0172 

-.0196 

-.019? 

-.0126 

6 

-.053 

• , 0202 

-.0173 

-.0158 

-.0130 

-.0117 

6 

-.050 

• , 0226 

-.0175 

-.Ol-s 

-.0131 

-.0116 

6 

-.(ipa 

-.029.7 

-.0170 

-.0196 

-.0118 

-.0095 

6 

-.0«5 

-.0229 

-.0211 

-.01X3 

-.0119 

• , 0 1 0 6 

h 

-.0113 

-.0376 

-.01911 

-.0133 

-.0107 

• , 0 u85 

6 

-.0113 

-.0376 

-.0194 

•.0133 

-.0107 

-.0085 

6 

-.0U3 

-.0376 

-.019- 

-.0133 

-.0107 

- . 0u85 

1 

-.0  12 

-.0051 

-.(1096 

-.0095 

-.00"  1 

-.  0039 

7 

-.010 

-.0060 

-.0090 

-.0037 

-.003? 

-.0029 

7 

-.00« 

-.0090 

..00«7 

-.0028 

-.0028 

-.0025 

7 

- .006 

• , 0 0 3 1 

- .0025 

-.0023 

-.0021 

-.0017 

7 

-.003 

• , 0 0 20 

• . 0 0 22 

-.001  ti 

-.001 7 

-.0017 

7 

-.002 

-.0020 

- ,00 1 9 

-.0012 

-.0011 

-.0010 

7 

-.0"2 

•,oo2o 

-.0019 

-.0012 

-.0011 

-,OU10 

7 

-.002, 

* , 0 0 20 

-.0019 

-.0012 

-.0011 

-.0010 

HG5  POLLOafLl  BY 

MBP  array  (i>i  Ffrn  AT 

mT 1 -t  • 

5 

(1 

6.200 

6.2000 

6,2000 

6.2000 

6.2000 

6,2000 

6.2000 

5 

5.500 

5,5000 

5.5000 

5.5000 

5.50O0 

5,5000 

5,5o00 

XR  VALUES!  I0‘'T»6)  asO  YH 

values 

(IomT.t)  aT  maotIme* 

6 

-.065 

-,02‘‘« 

-.0225 

-.ni-g 

-.0173 

-.0156 

6 

-.062 

-.0373 

-.0208 

-.0192 

-.0159 

-.0131 

6 

-.059 

-.0306 

-.0210 

-.0170 

-.0160 

-.0195 

6 

-.057 

-.0296 

-.0229 

-.0177 

-.0161 

-.0196 

6 

-.057 

- . 0266 

-.0296 

-.0162 

-.0197 

-.0133 

6 

-.051 

-.0959 

-.0212 

-.0162 

-."197 

• , 0 1 33 

6 

-.051 

- , 095u 

-.0212 

-.0162 

-.0197 

-.0133 

6 

-.051 

• , OuSli 

-.02' 2 

-.0162 

-.0197 

-.0133 

7 

-.012 

- ,0096 

-.noud 

.,0094 

- ,0090 

. ,0059 

7 

-.009 

- ,0060 

-.0037 

- . 0 0 3- 

-.0031 

• , 002e 

7 

• . 005 

• , no  33 

-.0093 

-.0028 

-.0026 

• , 0026 

7 

-.003 

-.0016 

- . 002o 

-.0019 

-.0020 

-.0021 

7 

-.001 

-.0010 

-.0013 

-.0011 

-.0013 

-.0019 

7 

.001 

0,0000 

-.0009 

• , 0006 

-.Ofl"8 

-.  0009 

7 

.001 

0,0000 

-.0009 

-.0006 

01)08 

-.0009 

7 

.001 

0,0000 

-.0009 

-.0006 

-.0008 

-.0009 

W'’'* 
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» 


»OLl.O*’ED  8V  (IN  FfFT)  4T  b 


a 

6.800 

6,8000 

6.5000 

6 , 5 C 0 6 

6,8000 

6,8000 

6,8UOO 

5 

8>600 

5,6000 

8.6000 

5,6000 

8. 6000 

5.600(( 

5.6000 

8R  VlLUtS(  tO'^Tab)  iNC 

VK  ViL'JEO 

( I0inT»7) 

4T  ►•iPTIxE* 

6 

-.C6<» 

-.6298 

-.021  1 

-.0194 

-.0161 

-.0196 

6 

-.063 

-.0332 

-.0211 

-.0177 

-.0162 

-,019b 

6 

-.060 

-.0266 

-.0220 

-.0162 

..0197 

-.0133 

6 

-.08“ 

-.0230 

-.0230 

-.0162 

-.0133 

-.OlCe 

6 

-.081 

-.0229 

-.0229 

-.0197 

-.0120 

- .0096 

6 

-■Q“8 

-.0382 

-.0169 

-.0133 

• . 0 I 2o 

-.oloa 

6 

- .0u8 

-.0352 

-.0169 

-.0133 

-."1 20 

-.0108 

6 

-.0«8 

-.0352 

- , 0 1 Ott 

-.0133 

-.0120 

-.0108 

7 

-.005 

-.0  022 

-.0022 

-.0029 

-.0023 

-.0023 

7 

-.001 

-.0012 

-.0011 

-.0016 

-.0017 

-.0016 

7 

.0O^ 

0.0000 

-.0006 

-.  0006 

-.0010 

-.Ou12 

7 

.00“ 

.OOOa 

,0009 

0 ,0000 

-.0005 

- . 0006 

7 

.005 

.0020 

.0012 

,0005 

.0002 

0 .0000 

7 

.006 

. 009  3 

."017 

. 3006 

. 0006 

, 0009 

7 

.006 

.0093 

.0017 

.0009 

.0006 

.0009 

7 

.006 

,0093 

.0017 

.0009 

,0006 

.0009 

6&9  fOLLO-EO  BV 

mhr  *ho»y  (16  rr^n  *t  - 

Tl-ta 

7 

a 

6.^00 

6.2000 

6,2000 

6,2000 

6.2OO0 

6,2000 

6.2OOO 

5 

5. “00 

5, 9000 

8,9000 

5.9000 

8,9000 

5 , uOOO 

5.9000 

«o  V4LUtS(  I(3‘‘Ti6)  anO 

79  V*Ul'P8 

(10^8,7) 

47  -iPTIffa 

6 

- . 065 

-.0297 

-.0211 

-.0199 

-.0178 

-.0162 

6 

-.089 

-.0307 

-.0211 

-.0177 

-.0162 

-.0197 

6 

-.086 

-.02*6 

-.0210 

-.0162 

-."197 

-.0133 

6 

-.080 

-.0207 

-.  0209 

-.0196 

-.0132 

-.0119 

6 

-.09“ 

- , 0205 

-.0207 

- . 0 1 8 1 

-.0119 

-.0107 

6 

-.059 

-.0276 

-.01*9 

-.0117 

-.0106 

-.Oo95 

6 

-.039 

-.0278 

-.01*9 

-.0117 

-.0106 

-.0995 

6 

-.039 

-.  0276 

-.01*9 

-.0117 

-.0106 

-.0095 

7 

.006 

.0026 

.0015 

.0010 

.0006 

.0003 

7 

.010 

.0093 

.0092 

.0015 

.001  1 

.0006 

7 

• Ct  0 

,0090 

.0030 

.001  9 

. 001  0 

.0007 

7 

.on 

,0099 

.0037 

.0023 

.0018 

.0015 

7 

.012 

.0051 

.0099 

.0025 

.0019 

.0015 

7 

.012 

,0079 

,0097 

.0027 

,"020 

.0017 

7 

.012 

.0079 

.0097 

.0027 

.(1020 

.0017 

7 

.012 

.0079 

.0097 

.0027 

.0020 

.0017 

6,8000 


6 


6 .?000 

7 


McB  roLUiioft)  dv  (i»,  rtfn  »t  « 

« (>.n''0  diOOdO  ft.^oco  dtPu^ft  ft,no''P  etOOOO  S.nwoo  b.PGOO 

5 s.jno  S.^000  5.^000  S.JUPO  ^,So"P  S.JOOO  S.iUOO 


KH  V»LUtS(  TO^-Tb#)'  Af'O  VH  V*LUeS  flflNT,?}  *T  “*PTImE«  6 


6 

• •OUb 

•,P1PS 

• • P 1 S b 

- . 0 1 '1 ) 

-.PlSO 

-.0110 

6 

• • OU^ 

•,P7tb 

-.OlSS 

-.0177 

- . P 1 1 b 

•.0105 

6 

-.CS« 

-.OlOb 

-.PlSO 

-.01  U 

-.010} 

-.009^ 

b 

• lOlB 

• , p 1 So 

•.Pl}b 

- . P 1 1 « 

-.0091 

• .P07«: 

6 

-.0^0 

•.PI lo 

-.01 }b 

-,00*9 

- . P070 

•.0077 

b 

• • 0 ?b 

•.OJ  75 

-.OlflB 

-.P07« 

- , pp  tq 

-.0080 

1 

b 

• 1 0 ?b 

• , P 1 75 

•.PIPS 

-.007« 

-.0079 

-.OObP 

b 

• • 0 ?b 

• .Pl7S 

-.OIOS 

-.0078 

-.0070 

-.00*0 

7 

.017 

.OObO 

.Pros 

, Oo  X7 

,PPAP 

.0075 

7 

• Olb 

,P07# 

.Ppso 

.00}7 

.POSl 

.0070 

7 

.01  7 

,P07b 

. P p Ml 

. 0 rt  A T 

.PPSl 

• 0u7e 

7 

.0  1 b 

, P ob  i 

.nPSi 

.PCuj 

.00*0 

. 0077 

7 

.PIS 

.0007 

.Pp«b 

.Pl'Ab 

,0070 

,0071 

7 

.0  1 B 

.00*5 

.onSl 

.POPS 

.P0S7 

.0079 

7 

.01B 

.POPS 

.0  1 

.0035 

.00}? 

.Ou?9 

7 

.010 

.POPS 

. . ) 1 

.POPS 

.00}? 

.0079 

“GB  fULOO^fO  H» 

khB  »b«*Y  (l>,  rt»T1  AT 

mT 1“E« 

0 

<1 

7.?P0 

7.?0"0 

7.?f/PU 

7 . 7 p p 0 

T.?noo 

7.7000 

7.7000  7,S000 

5 

b.^PO 

b.7t>00 

b.7PPli 

b.?OPP 

b,?OPO 

b.7000 

b.?000 

XU  VtOUtSdO'^Tbb)  ii.O 

7B  VALUES 

dlMTiT)  aT  -aPTIBE*  9 

b 

-.ObO 

-.P7b'» 

-.0710 

- , P7P  1 

-.0178 

-.0171 

b 

-.PS7 

•.PiOi 

-.071  } 

-.Pl*i 

-.Plb* 

-.0155 

b 

-.PSO 

•.07'»b 

-.P2«0 

-.PlhS 

-."1*7 

-.0119 

b 

• . 0 B 5 

-.0777 

-.o?po 

-.OIUS 

-.0117 

-.0170 

b 

-.057 

-.Pdu. 

-.0170 

-.01}} 

-.01?? 

-.0110 

t 

b 

-.OSb 

•.P7B« 

-.Pl«7 

-.0119 

-.0107 

• , 0 1 b 1 

b 

• . 0 ?b 

•.P7B« 

-.7107 

-.0119 

-.0108 

•.eibi 

b 

• .OXb 

• .P7Jb 

-.PI  “7 

-.0119 

-.0107 

-.Olbl 

7 

.077 

.0100 

.0P7B 

.0005 

.0057 

,00U9 

7 

.07b 

.PliS 

.00*0 

. 0 0 b 0 

,005« 

,0050 

7 

.077 

.PlSO 

.00®0 

.0007 

."0*7 

,0050 

7 

.0>S 

.011} 

.0007 

. OObb 

.0055 

.OOOb 

7 

.072 

.OOB# 

.Ofloe 

.0005 

.O05« 

,0089 

7 

.07} 

.Ploe 

.0079 

.OJbl 

.0051 

.0075 

7 

.0?} 

.0l«6 

.0079 

.0001 

,"051 

,0075 

7 

.071 

.01B« 

,0079 

.OObt 

.0051 

.0075 

1 

0grWMlABl.t  WP'( 
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HG®  ^OLLCfED 

My®  I®0*r  (1.)  r 

eCT)  4T 

IP 

a 

7.000 

7.0000 

7,0000 

7 , OOOo 

7.P0P0 

7.rooo 

7.0W00  7,o000 

5 

ft.ono 

6 . OUOO 

0 , OoOu 

o.Ooop 

6 , PflOP 

6 . Ou  PO 

0 . 00  00 

*R  V*LUtStinMT»e) 

''ii.Uf.  S 

f TOf  7, 

7)  *7  “iPTiME*  10 

6 

- .0UI> 

•.nj®6 

-.017® 

• . O'l  fcy 

-."1®" 

-.0115 

6 

-,"226 

.,0170 

-.OUT 

-."125 

-."111 

6 

••0  !7 

-.0107 

• , 0 1 50 

-.1)1  22 

-."1  12 

-.0101 

6 

• .012 

-.0170 

-.0112 

-.016® 

-.ppo® 

-,PO®0 

6 

• .07<) 

-.OHl 

-.0157 

• . 0 0 ® B 

-.PP*7 

• , 0u76 

6 

• .02U 

-.0165 

-.010® 

- .pp*5 

- . P07b 

• , u u b2 

6 

• . P ?M 

• , 0 1 65 

-,010® 

-.00*5 

• .P07h 

• , 0 U 6 fl 

t> 

• .02<i 

• .0 1 65 

-.01 ou 

-.00*5 

•."076 

• . 0 U B B 

7 

.P?® 

.00*6 

.0076 

. 0 0 6 u 

.Pp'ia 

.0015 

7 

.021 

.010® 

,Oo7® 

.00*5 

, POUP 

,00*1 

7 

.P21 

.00®® 

,00*0 

. 0 P*  U 

."065 

. 0 u 17 

7 

.020 

.O0®6 

, Op67 

.0051 

. P 0 11  > 

, 0 0 1* 

7 

.OI« 

.0C«I 

.00*1 

. Pa®  ® 

. " " ® 1 

. 0 u '5 

7 

.0  16 

.OlPi 

,00*0 

.P0«5 

, 0 0 7 J 

.0031 

7 

.016 

,0101 

,0(160 

.00®5 

.001® 

• 0 0 53 

7 

.P16 

.0101 

,0060 

.00®5 

.001® 

.('Oil 

MG®  fOLLO.^L  Bv 

•.yo  twRAY  (I*;  rtrri  *t 

BT I ®fc» 

11 

u 

S.POO 

5.OP00 

, OO^P 

5.BOPO 

5 , »0PO 

5 . B 0 G C 

5.8000  5, *000 

5 

5.000 

5 . ooOu 

5, *000 

5 .Boon 

<.*PPp 

5, *000 

5. *000 

*«  v*LUtS(  It)NT66)  »*.0  7® 

y‘Lt'ES 

f IP'.T,' 

T)  47  "tPTiBr*  11 

6 

- .0^“ 

-.021 1 

-.01°® 

- . p 1 0® 

-."15b 

- , 0 1 u 1 

6 

• . 066 

-.0265 

-.0106 

-.015J 

• ,Pl®ft 

-.0126 

6 

..pao 

-.0225 

-.0100 

-.01  M 

-.0126 

- . 0 1 1 * 

6 

-.017 

-.015® 

-.0ia7 

-.0121 

-.011? 

-.0102 

6 

-.013 

-.0157 

-.0l7i 

-.PUP 

- . "0°® 

• , 0*)*® 

<> 

-.020 

-.01®® 

-.0110 

• . 00®l< 

..P0»* 

-.007® 

6 

-.P20 

-.01®® 

-.0110 

• , P p 0 * 

..OP** 

-.007® 

6 

-.02B 

-.01®® 

-.01 1 ® 

• . Pc®? 

..P0*« 

-.007® 

7 

.026 

.00«e 

,op72 

.PP5S 

."075 

."0  76 

7 

.p>2 

.010® 

.007® 

. 0055 

.O0®5 

.0017 

7 

.021 

.0105 

.0075 

.0051 

. 0 0 ® 6 

.Ool® 

7 

.021 

,00®6 

.006® 

.00*2 

.oo®l 

.0015 

7 

.01® 

.0001 

.00*5 

.00®® 

.P0®2 

.OOJb 

7 

.010 

.0115 

,006® 

.00®5 

, " 0 1 e 

.0011 

7 

.010 

.0115 

,006® 

.00U5 

•.001® 

.foil 

7 

.010 

.0115 

.006® 

.00®5 

."PI® 

.0011 
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HG’’  F0LL0»FD  PV  tBB*Y  tiK  Ffm  «T  ►■TP'ta  l2 


*1 

fat?no 

fa.?0P0 

h.?0OU 

6 . ?0  0(1 

0.8000 

t^?000 

fa , 8000 

5 

5.100 

5,1000 

5.1000 

5.1000 

5.3000 

5 • 0 0 

5.1000 

XB  V»  LUt  S ( t O'*!**!!  »^0 

YB  ViLl'El 

( t 0 M T , 7 ) 

iT  “iPTlwr* 

6 

..cun 

-.0186 

..01  ue 

- . 0 1 IQ 

- . 0 1 ?B 

-.flllB 

6 

..0J5 

-.OlfaB 

.,01  ue 

-.018u 

-.0115 

-.0105 

6 

-,pJB 

-.0165 

-.0150 

-.omo 

-."C’l 

-.0081 

6 

-.0^5 

-.0156 

-.Olio 

- . 0 1) « 7 

-.nn«i 

-.0074 

6 

-.071 

-.nioi 

-.01  15 

- . 0 (,  7 7 

- . O070 

- , 00  0 1 

fa 

-.018 

-.01?! 

-.00B8 

-.005B 

-."001 

-.0055 

fa 

-.PI8 

-.Olii 

-.O0B8 

-.0C5Q 

-."001 

-.0055 

fa 

-.019 

-.ni?l 

-.noB8 

-.005Q 

-.0061 

-.0055 

7 

.011 

.0  115 

.OOBB 

. no  ‘u 

.0071 

.0001 

7 

.0>B 

.0187 

.0108 

.0080 

, OOOB 

.0058 

7 

.C?5 

.0151 

.00B5 

.0007 

.0057 

.Ou-e 

7 

.0?.i 

,01  lu 

.OOBl 

.OCfal 

. 0058 

,0uui 

7 

.e^o 

,0086 

.OOBi 

.0058 

, 0 0 UB 

.OOul 

7 

.018 

.OlOB 

. 007  7 

. noufa 

,"nuu 

,0017 

7 

.018 

.01  OB 

. 00  77 

.OilUb 

, "0  uu 

,0017 

7 

.0)  » 

.OlOB 

. 0077 

. OOUfa 

. nouQ 

.0017 

MCB  xOLLO.rO  bV 

wbb  »faB»Y  (I-Y  rtrn  *T  ►‘TI-E* 

11 

0 

fa.faOO 

fa , fa  0 0 1) 

• ftii  0 

0,0000 

fa  . 7 0 0 0 

6.7000 

5 

5.800 

5.8000 

5.8000 

5 , • 0 0 0 

5.8000 

5.8000 

XB  VlLUtS(  ID‘iT«6)  »vO  ' 

r iONi»7) 

17  ~»PTIME« 

fa 

-.OUB 

-.0855 

-."810 

-.0811 

-.080] 

-.01«1 

fa 

• . 0u2 

-.0851 

-.niBB 

-.0170 

-."lOB 

-.015B 

fa 

-.017 

-.08u7 

-.01*0 

-.OlTu 

-.0151 

-.0111 

fa 

-.011 

-.01B7 

-.017B 

-.01U5 

-.0180 

-.01  OB 

fa 

-.c?e 

-.Ol5u 

-.0150 

-.01"* 

-."111 

-.0108 

fa 

-.07u 

-.niB3 

-.0115 

- . 0 v;  B fa 

.."070 

- , OOoo 

fa 

-.0?u 

-.OlBi 

-.0115 

-.O.lBfa 

- . 0n7B 

-.0006 

fa 

• . O^u 

-.OlBl 

-;nil5 

- . 0 0 B fa 

-.0,170 

- . Oofae 

7 

.051 

, njue 

.01*B 

.0105 

. 0 1 U 6 

,0186 

7 

. 0u7 

.0851 

.OlTB 

. 0 1 uu 

.0187 

.0111 

7 

.cui 

.0855 

.0110 

.0151 

.0180 

.OOBi 

7 

.OlB 

.0811 

.0171 

.0110 

.0102 

.0076 

7 

.Oil 

.0171 

.0150 

.0101 

.0007 

. 0086 

7 

.P?B 

.081U 

.OlUO 

.0005 

.OOB] 

.0085 

7 

.0»B 

.081U 

.otuo 

.OOBl 

,0001 

.0085 

7 

.OJB 

.081U 

.OlUO 

.00B5 

.OOBl 

.0085 

fe.JOOO 


1? 


A. 7000 


11 
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HG»  »'OLLf’«fO  KB*’  iBKiV  ( P'  FfFT)  IT  MT;Mts  1 <* 


« 

5. BOO 

5,0000 

5.800U 

5 . 9 0 n f, 

6 .noon 

6,1000 

r.20oo 

5 

5tB00 

5.6000 

5.6000 

5.6CO0 

5.60O0 

5.6000 

5.6000 

XP  V»LUtS(  ID'*T»B)  7H 

v»LUt  S 

noM»7) 

at  maPTIkE* 

6 

-.0225 

-,Ol»B 

-.0179 

-.0171 

-.0165 

6 

-.OlB 

-.0255 

-.0157 

-.019) 

-.0195 

-.0196 

6 

-.052 

-.0201 

-.0150 

-.0162 

-.019? 

-.0121 

6 

-.0 1 06 

-.0195 

-.01X9 

-.6129 

-.0121 

6 

-.0?5 

-.0l«7 

-.0152 

-.0  110 

-.0115 

-.0007 

6 

«.0?2 

-.01«i 

-.0126 

-.0099 

-.010? 

-.0105 

6 

-.022 

-.0185 

-.0126 

-.0099 

-.0102 

-.0105 

6 

-.022 

-.0185 

-.0126 

• , 0099 

-.0102 

-.0105 

7 

.OOP 

.0215 

.0165 

. 0 1 9(1 

.0129 

.0106 

7 

.OuO 

.0259 

.0192 

.0122 

.0117 

.0102 

7 

.05e 

.0207 

,0190 

.0156 

.0107 

.00*5 

7 

.052 

.01P3 

.0199 

.0117 

.0100 

, 0085 

7 

.029 

.0152 

.0197 

. 0099 

. 009} 

.007] 

7 

.025 

.Ol’O 

.0126 

.00«9 

.00*6 

.0092 

7 

.025 

.01^6 

.0126 

.00*9 

.00*6 

.0  0, *2 

7 

.025 

.01<»6 

.0126 

,00*9 

.00*6 

.00*2 

KG«  xOLLO-FO  87 

Me9  iRRxv  fl9  xfrn  »t  ^ 

T I-t» 

15 

u 

9.500 

9. JOOO 

9.50O0 

9 , 5oon 

9 ,7o0n 

9,9000 

5.1000 

5 

5.200 

5.2900 

5.2000 

5.2000 

5.2000 

5. 2000 

5,2000 

XR  ViLUfcSt TONTiB)  ‘NO  7R 

V»U"ES 

f iri»(T,7) 

at  ••'aPTIkE* 

6 

-.090 

-.0222 

-.019- 

•.oi«s 

-.019| 

-.01*1 

6 

-.056 

-.0269 

-.0175 

-.0168. 

-.0] 7fc 

-.0169 

6 

-.051 

-.0255 

-.0170 

-.01*2 

- .0 1 60 

-.0159 

6 

-.050 

-.0205 

-.017 

-.0152 

-.0197 

-.0190 

6 

-.029 

•.0165 

-.0170 

• ,0 1 2* 

-.0152 

-.0115 

B 

-.025 

-.0218 

-.0195 

-.0115 

-.0110 

-.0102 

B 

-.025 

-.0218 

-.0195 

-.0115 

-.0119 

-.0  i 12 

B 

-.025 

-.0216 

-.0195 

-.0115 

-.0119 

-.0102 

7 

.099 

.0296 

.0l9u 

.0166 

.0160 

.0191 

7 

.09- 

.0509 

.0179 

.0157 

.0195 

.0127 

7 

.090 

.0267 

.01*0 

.0169 

.01X0 

.0112 

7 

.057 

."259 

.0179 

.0197 

.0127 

.01  10 

7 

.055 

.0187 

.01*2 

.0129 

.Olio 

.Ou95 

7 

.052 

.0250 

.0155 

.0115 

.0111 

.00*6 

7 

.052 

.0250 

.0155 

.0115 

.0111 

.0086 

7 

.052 

.0250 

.0155 

.0115 

.0111 

.0086 

B,?000 


lu 


^.1000 


15 
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UNCLASSZFXeO 


COASTAL  STUDIES  INC  COLLCOE  STaTiOn  TX  F/t  A/2 

OEVELOf»NENT  OF  SUROC  II  PR06RAM  xlTH  APPLICATION  To  THE  SABINE— ETC (U) 
NOV  77  R 0 REZO*  A C VASTaNO*  T J RCZO  OACW6A-7A-C-OOiS 


CERC-TP-77-1S 


NL 


A 


MC»  rOLl.'3"fO  BY  MB«  4MR*Y  (IW  FtrT)  *T  HTIMt*  |6 
« «>6nO  M.BOOO  U.BonO  «*T(i(>o  <>•‘*000  S»nOOO  5.1000  5.1000 

5 5.100  5.1000  5.1000  5.tuoo  R.lQOo  5.1000  5.1u00 

KH  V4l.U|3tI0NT«*J  iNO  »H  y*LIIE8  (lOxTiTJ  *T  'UPTIME*  t6 


ft 

■ . 0?7 

-.0156 

-.0199 

-.0162 

-."179 

-.0160 

ft 

• . 0 ?ft 

-.0^01 

-.0158 

-.0150 

-."150 

-.0160 

ft 

•.0?5 

-.Pl»« 

-.0190 

-.0150 

-.0158 

-.0125 

ft 

• .0?1 

-.0161 

-.0155 

-.015* 

-."178 

-.0115 

ft 

» . 0?  1 

-,oi50 

-.0192 

-.0104 

-.0117 

-.0125 

1 ‘ 

-.019 

-.0155 

-.0152 

-.0100 

-."110 

-.0099 

ft 

-.019 

-.0155 

-.0152 

-.0100 

-.0!  in 

-.0094 

ft 

-.019 

-.0155 

-.0152 

-.0100 

-.0110 

- , 0094 

7 

.057 

.0295 

.0297 

.0259 

."250 

.0190 

7 

.059 

.0578 

.0210 

.0222 

."198 

.0140 

7 

.051 

.0511 

.0290 

.0222 

.01*1 

.0199 

7 

.095 

,0290 

.0295 

.0206 

.0169 

.0116 

7 

.090 

.0271 

.0227 

.0161 

.0155 

.0199 

7 

.oia 

.0260 

.020 

.0195 

.0190 

.0110 

7 

.05B 

,0268 

.021  > 

.0198 

.0190 

.0110 

7 

.010 

.0266 

.0211 

.0198 

.0190 

.0110 

HCB  POLLn»tD  BY  MpQ  4MR*Y  (Jw  r£rTl  *T  MTlMt»  l7 
a 5.000  5. (1000  5.0000  5.60O0  O.Bcoo  O.CQOO  B.oooo  ft. 0000 

5 5.S00  5.J0O0  5.5000  5.1000  5.5000  5.1000  S.lUOO 

*»  V4LU£S(I0NT»ftJ  *«.0  YB  V41.UE9  flONTiT)  *T  MiPTlMC*  >7 


-.015 

-.0119 

-.01*9 

-.0158 

-."129 

-.0119 

-.015 

-.0125 

• ,00*9 

-.0111 

-.0121 

-.0151 

-.011 

-.0118 

-.0069 

-.0111 

-."118 

-.0122 

-.019 

-.0111 

-,0108 

-.0109 

-,0101 

-.0096 

-.015 

-.0097 

-.Ot"l 

-.00*6 

-,"099 

-.00«5 

-.015 

•.0101 

- . 0092 

-.0,0*6 

-,00*6 

-.0085 

-.01  5 

-.Old 

- , 0092 

-.00*6 

- , 00*6 

-.0085 

-.Oil 

-,0101 

-.0092 

- . Oi)*6 

- , 0 0 * 6 

-.00*5 

.071 

.0518 

,0*66 

.1178 

.05"* 

."291 

.Oft* 

.0*6* 

.027* 

.02*9 

.02*5 

.02*1 

.056 

.0*11 

.02*0 

.02*9 

.0269 

.0219 

.055 

.0162 

.011* 

.0269 

."227 

,01*9 

.0*9 

.0117 

."291 

.0211 

.0210 

,0190 

7 

.Oft* 

.0116 

.0252 

.0211 

.019^ 

.0175 

7 

.0*4 

.0116 

.0252 

.1211 

.019$ 

.0175 

7 

• Oft* 

.0116 

.0252 

.0215 

.0193 

.0175 

BBI  AVAIUBIE  COPY  ^ 
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Hc#  ^ouuo-ED  ay 

H6R  *rtR*V  (IM  FEFTl  *T 

nTI-e* 

18 

< 

«>7oo 

6.7000 

6, 7000 

6,9o0o 

9,1000 

5.3000 

5.6000  5,6000 

s 

5«000 

S.OOOO 

5,0000 

S.OOOn 

5,0000 

S.OoOo 

5.0VO0 

XR  VtLUlSC  IOMT*«)  tK'P  YN 

VXLUf S 

r jont* 

7)  *T  WARTIME*  IS 

6 

.002 

0,0000 

• , 0007 

•.00?0 

•,0020 

•,0060 

6 

.001 

• , 00  06 

•.0011 

•.0022 

•.00)0 

• .0037 

6 

O.OOO 

•.0011 

•.0026 

•.0026 

•.00)1 

•.oust 

« 

-.001 

• ,0016 

•.0016 

•.0026 

• , 0 0 2 0 

• . 0028 

k 

•.001 

•,0015 

•.0022 

•.002) 

•,002* 

•,002o 

6 

• . 002 

•,0016 

•.0020 

•. 002l 

•,0026 

•.0025 

6 

• . 002 

• ,0016 

•, . 0o20 

•.Oo?i 

•.0026 

•,0025 

b 

• . 002 

•,0016 

• . 0 0 2o 

•.0021 

•,0026 

• ,0025 

7 

.OUO 

,0269 

.0212 

.0226 

.0228 

.0226 

7 

.0<t5 

,0))2 

.0211 

.021  1 

,0210 

.0206 

7 

.0«0 

,0310 

.0210 

.019) 

.019? 

.0175 

7 

.010 

,0266 

.021  1 

.019} 

,0176 

,0160 

7 

.OIS 

,021  1 

.0211 

.0161 

,0160 

.0131 

7 

.0)1 

.022R 

.0193 

.0166 

.0165 

.0116 

7 

.0)1 

,0226 

.0193 

.016« 

,016) 

.0116 

7 

.0)1 

,0229 

.0193 

.0166 

,0169 

.0116 

NCR  YOLLOk{0  Sy 

hbR  *HR*y  (IN  FfFTl  *T 

MTr’Ei 

1» 

« 

}.200 

3,2000 

5,200U 

).5C00 

),*ooo 

6.O000 

6.1000  6.1000 

5 

X.ROO 

u.uq'ko 

6.600(/ 

6.6000 

6,6000 

6,6000 

6.6000 

XR  Y*LIJt8(10NT«*l  iNO  rH 

Y*UI'ER 

(I0nT,7}  *7  19 

.016 

,Oo75 

.0056 

.0090 

,00)9 

,0025 

.01<l 

,O0R6 

,0066 

.00)7 

,002* 

,0019 

.012 

,0072 

,0050 

.00)0 

.0.12? 

.0013 

.011 

,0061 

,0063 

.00)0 

.0017 

.0007 

.006 

,0066 

,00»5 

,00?o 

,0012 

,0u06 

.000 

,0066 

.00)2 

.0015 

.0011 

.0007 

.000 

,0066 

,00)2 

.0015 

."Oil 

,0007 

.006 

,0066 

,00)2 

.0015 

.0011 

,0uO7 

.0)7 

,0301 

.0262 

.0260 

.02*6 

,02*6 

.0)) 

,0617 

.0266 

.0266 

.0207 

.0266 

.0)2 

.0371 

.0200 

.0267 

.0268 

,0266 

.0)0 

,03«R 

.0)07 

.026? 

."26* 

.0212 

.0«7 

.0306 

.0266 

.0229 

."2»9 

' .0230 

7 

.06) 

,0326 

,0)06 

.021 1 

.021 1 

.0212 

7 

.06) 

.0326 

.0)06 

.0211 

.021  1 

.0212 

7 

.06) 

.0126 

,0306 

,0211 

,0211 

.0212 
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H6»  FOtLOKfD  BY  Mart  *Rrt*r  (IN  rEfT)  4T  wTIhfc*  go 


a 

2.70U 

2.7000 

2.7000 

2.9000 

1.2000 

1.900U 

3.6000 

9 

««0't0 

9.0000 

9.0000 

9.0000 

9.0000 

9.0000 

9.0900 

MU  Y4LUeS(I0^Ta«)  4N0  YR 

VFLoes 

(lONTaT) 

. 47  M4P7JHJB 

6 

lOIS 

,0086 

.0071 

• OOYfc 

.O0T5 

.007) 

* 

• 01$ 

.0111 

.0069 

.006} 

.0082 

.0019 

6 

• Oitt 

.0082 

.0060 

.0015 

,"059 

.0051 

« 

.012 

,0078 

.0081 

.0019 

.0096 

,00)9 

6 

• on 

.0060 

.0017 

• 0092 

."0)6 

,0026 

6 

• OOR 

.0062 

.0099 

.0012 

."")) 

• 00)9 

« 

• OOR 

.0062 

.0099 

.0012 

.00)) 

.00)9 

6 

.oo« 

,0062 

.0099 

.0012 

."0)) 

.0019 

T 

• Ol6 

,0191 

,01«0  . 

.0198 

.0217 

.02)6 

7 

.0$$ 

,0266 

.0186 

• 0119 

."2") 

.0222 

7 

.013 

.0211 

.0101 

.0169 

.0187 

.0201 

7 

.011 

.0216 

.0201 

.01*7 

."172 

.0118 

7 

.029 

.0181 

.01*6 

.0117 

."111 

.0191 

7 

.n?7 

.0201 

."171 

.0110 

."199 

.0119 

7 

.0*7 

.0201 

.0171 

• Olio 

."|99 

.0119 

7 

.fl|7 

.0291 

.0171 

.0110 

.0199 

.0119 

MG*  FOLtOaro  BY 

MBR  4rtR6Y  (IN  FEfT)  »T  MTI«e« 

21 

a 

1.700 

1.700U 

1.7000 

l.looo 

).9e"o 

9.0000 

9.1000 

5 

a.600 

9.6000 

9.60O0 

9.6000 

9.6000 

9.6000 

9.9000 

XR  VtLUESdONTa*)  4N0  YH 

V4l,l'£S 

f IOnTbT) 

4T  M4PTlrtE6 

6 

.016 

.0079 

.0062 

.OOYl 

,"070 

.00*1 

6 

.011 

,0091 

.0062 

.008} 

."067 

."070 

6 

• 012 

,0972 

.0060 

.0051 

.0055 

.0018 

6 

.011 

.0070 

,0069 

• 0061 

.0008 

.00)7 

6 

• 010 

,0011 

.0011 

.009) 

.30)9 

.Ou)l 

6 

.000 

.0061 

.0091 

.0017 

."Oil 

.00)0 

» 

.008 

.9061 

.0391 

.00)7 

.0015 

.00)0 

6 

.008 

.0061 

.0091 

.00)7 

."015 

.00)0 

7 

.029 

.0191 

.0119 

• 0161 

.0169 

.0191 

7 

.026 

.0190 

.0119 

.oi5n 

."169 

.0182 

7 

.026 

.0161 

,011(1 

.01)8 

.015) 

.0168 

7 

.02$ 

.0169 

.0166 

.0167 

."190 

.0119 

7 

.021 

.0118 

.0118 

."1.26 

."128 

• 0116 

7 

.022 

.0167 

.0126 

• 0119 

."111 

.0109 

7 

• 022 

.0167 

.0126 

.0119 

."HI 

.0109 

7 

• 022 

.0167 

.0126 

• 0119 

."111 

.0109 
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n(;<»  rOttOxfO  d»  *<8'»  firr)  at  2ir 
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APPENDIX  E 


LISTING  OF  KEY  ARRAYS  AND  CHANNEL  AND  BARRIER  PLOT 
FOR  SABINE-CALCASIEU  REGION 
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APPENDIX  F 


IDENTIFICATION  OF  GAGES 
FOR 

SABINE-CALCASIEU  ASTKOTIDE  CALIBRATION 

Gages  for  Sabine-Calcasieu  astrotide  calibration  and 
time  sequences  of  accepted  astrotide  simulation  at  those 
gages  for  72  hours  are  identified.  Also  included  are 
listings  of  the  channel  output  at  t = 30,  60,  and  90 
hours.  For  explanation  of  each  column  see  Appendix 
C,7,b. 
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*8TBU  TIOt  C^LIBHiTIOH  FOO  8 »B J NE-c AUC »S lEO  A^EA 
SABInE  bass  tides  used  as  riPMT 

PERIOD  OF  RECORD-  OPOO  *U5,?2  TD  ?<|00  AUC,£6.IR7l 
calcucatioss  allow  FOR  suB-ORio  scale  channels  and  Barriers 

TIHE  StOuENfLs  OF  -ateR  LEVEL  anq  FLO-  ARE  SAVED  FOR  THE  FOLLO»lNc  PLACES- 
GAGE  I SAblNF  Pass.  SOUTw.ejt  jetty 
gage  2 PORT  ARTHUR.  CE  area  OFFICE 
GAGE  3 NORTH  SABINE  LAKE 

GAGE  « bCAUHONT.  NECHtS  RIVER  AND  BRAKES  BatOu 

GAGE  S orange  naval  STATION.  SABINE  RIVER 

GAGE  6 CAMERON,  CALCASIEU  PASS 

GAGE  7 HAGKBEWRy.  CALCASIEU  RIVER  AND  PASS 

GAGE  S !.«,►,  at  CALCASIE.U  LOCR.  hEST 

GAGE  R LAKE  CHARLES.  CALCASIF.u  RTvFR 

FLO*  1 SABINE  Pass  InflOh 
FLO*  2 CALCASIEU  PASS  InFlO. 

FLOW  3 FLO-  TO  nfCheS  RIver  from  saBIne  L*ke  and  INTbaCO‘STal  wAteR-av 
FLU*  a IaST-aRn  flO*  via  IntRaCOASTaL  CANal  JuST  EAST  OF  SabInj  rivIp 
FLO-  b FLO*  TO  SAbl'.E  PlVfP  from  SAPiNf  LAKE  ANU  INTBaCOASVAL  *ateb-ay 
FLO*  0 FLO*  TO  CALCAsSItu  l-TvFR  fRQh  CaLCaSIEU  L*KE  ANO  INT»ACOaSTaL  *. 
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APPENDIX  G 


IDENTIFICATION  OF  GAGES 
FOR 

SABINE-CALCASIEU  HURRICANE  CARLA  VERIFICATION 

Gages  for  Sabine-Calcasieu  Hurricane  Carla  verification  and 
time  sequences  of  water  level  and  flow  at  the  identified  gage 
for  60  hours  are  identified.  Also  included  are  listings  of 
detailed  channel  output  at  30  and  60  hours.  For  explanation 
of  each  column  see  Appendix  C,7,b. 
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H'jfcRJCtkfc  C4»L*  C*LI«R»TU'M  fob  S*■I^t•C*LC^SIeo  4«E* 
ptHicrO  OP  "pcnpo.  oi'o*'  St®  10  ’o  eeoo  sep  ij«  ippi 
C*UCL'L* TIOnS  4LLC*  ''ts  SoB-GBTO  SC4LE  C"4^SELS  4N0  P4RHIE95 

Tt-E  SfCuE‘-CtS  OF  *tU*'  UrvFL  4^0  FLO*  4><E  S4VEO  FOP  The  FOLLO'^inq  PL4CE8- 
G4C>E  1 S43II-F  P*Sb.  south. est  JETtt 
C4&E  2 POST  ikThUP#  CE  40E4  OFFJCE 
G4GE  3 nOPTh  S4Hr«£  L4KE 

G4GC  tt  aE4U'*QNT,  nEChES  PIVE*  4NC  bP4hES  B4Y0u 
G4GE  S OPAkGE  N4V4L  STATIqh.  SABIhE  PlvED 
CAGE  b C4heGok«  CALCASIEU  PASS 
GAGE  T *fst  E'.o  of  InTBaCOASTal  hjiTERhA.T 

GAGE  s SAui <E  Pass,  coast  guapo  station 
CAGE  p LA<F  Charles.  CauCaSIFu  PIvEP 

Fuo»  j SABINE  Pass  ineuoh 
FLO*  I CALCASIEU  BASS  InFlO* 
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